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SUMMARY 


DA  Task:  3A99-27-005-08  Micrometeorology.  (USAEPG) 
Title:  Studies  of  the  Three-Dimensional  Structure  of  the 

Planetary  Boundary  Layer 

Contractor:  University  of  Wisconsin.  Madison,  Wisconsin 


The  broad  objective  of  this  research  is  to  characterize  the  three- 
dimensional  structure  of  the  planetary  boundary  layer.  A  series  of 
eight  investigations  is  reported,  which  fall  into  three  types:  (1)  a 
comparison  of  experimental  measurements  against  a  theoretical  model 
for  the  effects  of  solar  heating  at  the  earth' s  surface,  (2)  experimental 
measurements  of  the  effect  of  variations  in  surface  roughness  on  the 
dynamics  of  the  air  flow  near  the  ground,  and  (3)  a  theoretical  analy¬ 
sis  of  the  steady-state,  neutral  wind  profile  for  the  entire  planetary 
boundary  layer. 

The  results  may  be  summarized  as  follows: 

(1)  Data  are  presented  for  the  thermal  response  of  an  air-soil 
system  to  artificial  cycles  of  solar  radiation.  It  was  found  that  a 
theoretical  two-layer  (air-soil)  model  did  not  adequately  describe  the 
experimental  results  and  that  the  grass  cover  at  the  site  must  be 
treated  as  a  third  layer.  The  major  discrepancy  in  magnitude  of  pre¬ 
dicted  versus  observed  amplitude  ratios  of  soil  heat  flux  to  net  radi¬ 
ation  was  attributed  to  the  evaporation  factor.  From  this  it  is  evident 
that  thermal  response  experiments  represent  a  promising  method  for 
determining  evaporation.  For  the  investigation  of  two-layer  conduc¬ 
tion  processes,  it  is  recommended  that  thermal  gradient  experiments 
be  conducted  over  surfaces  which  are  free  of  both  moisture  and  vege¬ 
tation. 

(2)  Results  are  reported  from  other  studies  and  experiments  in  a.  _ 
broad  attack  on  the  evaluation  and  significance  of  the  aerodynamic 
roughness  of  the  earth' s  surface.  Wind  profile  data  from  the  litera¬ 
ture  were  re-analyzed  with  special  emphasis  on  measurements  above 
tall  vegetation,  and  a  regression  was  established  between  vegetation 
height  and  aerodynamic  roughness.  A  review  of  recent  literature  on 
drag  coefficients  for  vegetation  is  included.  Data  on  the  areal  dis¬ 
tribution  of  plant- cover  types  were  evaluated  for  the  state  of  Wisconsin 
and  for  latitudinal  zones.  From  these  data,  estimates  were  obtained 
for  regional  and  seasonal  distributions  of  vegetation-produced  surface 
roughness  by  using  an  empirical  relationship  between  plant  height  and 
the  aerodynamic  roughness  parameter.  A  least  square  method  is 
described  and  a  program  is  outlined  by  which  wind  profile  parameters 
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can  be  evaluated  from  detailed  anemometry,  with  the  aid  of  electronic 
digital  computers.  Finally,  a  series  of  controlled  experiments  was 
conducted  to  study  the  effect  of  varying  the  area  density  of  defined 
obstacles  on  the  aerodynamic  roughness  and  other  wind  profile  param¬ 
eters. 

(3)  Existing  models  of  mean -velocity  distribution  in  fully  turbu¬ 
lent  duct  flow  are  critically  reviewed  in  terms  of  analytical  expres¬ 
sions  for  the  length-scale  of  turbulence.  A  new  theoretical  model  for 
the  length-scale  of  turbulence  which  satisfies  the  principle  of  simi¬ 
larity  between  duct  and  atmospheric  boundary  layer  flow  is  given.  A 
comparison  of  the  velocity  distribution  in  duct  flow  predicted  by  the 
model  for  smooth  or  rough  walls  shows  satisfactory  agreement  with 
observations. 

The  major  result  of  the  study  is  a  generalized  mathematical  model 
for  steady-state  flow  in  a  neutral  atmospheric  boundary  layer.  Evi¬ 
dence  is  also  presented  which  indicates  the  von  Kdrmdn  relation  is  a 
constant  (0.428)  only  for  fully  developed  flow  in  ducts  with  rough 
walls.  A  suggestion  is  made  that  the  turbulent  viscosity  is  adjusted 
so  that  the  "effective"  Reynolds  number  is  less  than  a  critical  value 
(approximately  1,000),  and  is  independent  of  the  Reynolds  number. 

The  investigator  concludes  that  certain  aspects  of  his  new  mathe¬ 
matical  expression  for  the  mean  length-scale  of  turbulence  merit 
closer  investigation.  The  interpretation  of  a  characteristic  parameter 
(m)  for  duct  and  atmospheric  boundary  layer  flow  as  indicating  "degree 
of  freedom"  (or  number  of  dimensions)  makes  it  similar  to  the  structure 
of  probability  functions.  This  feature  may  help  to  consolidate  a  dual 
picture  of  turbulence,  in  building  a  bridge  between  engineering  and 
theoretical  concepts. 

(4)  The  new  theoretical  form  for  the  length-scale  of  turbulence 
as  a  function  of  distance  from  the  lower  boundary  is  applied  to  a  baro- 
tropic  atmosphere.  As  a  consequence,  all  characteristics  of  the 
boundary  layer  are  determinate  upon  specification  of  three  parameters: 
geostrophic  wind,  Coriolis  force,  and  aerodynamic  roughness.  The 
theory  is  compatible  with  the  conventional  logarithmic  wind  profile  in 
the  surface  boundary  layer  and  computed  wind  spirals  compare  satis¬ 
factorily  with  observed  wind  spirals  for  the  boundary  layer  of  the 
atmosphere.  Applications  of  the  theory  for  estimating  the  dissipation 
of  energy  in  the  lower  troposphere  are  discussed. 

No  major  conclusions  or  recommendations  other  than  those  stated 
above  are  given  by  the  contractor. 

METEOROLOGY  DEPARTMENT 
USAEPG 
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General  Introduction 


Heinz  H.  Lettau 

Department  of  Meteorology 
University  of  Wisconsin 


The  understanding  of  the  kinematics,  dynamics,  and  energetics 
of  the  natural  small-scale  processes  near  the  earth' s  surface,  and 
the  structure  of  the  lower  atmosphere,  is  the  prerequisite  for  any 
successful  application  of  miaometeorological  research.  The  follow¬ 
ing  brief  statement  of  the  general  technical  requirements  of  the  con¬ 
tract  work  summarizes  the  principal  objectives  of  the  research 
program: 

"Experiments  shall  be  performed  by  the  contractor  to  study  the 
horizontal  variations  of  micrometeorological  factors  as  they  relate 
to  surface  characteristics.  The  bulk  of  the  experimental  effort  shall 
be  an  intensified  Investigation  of  a  limited  geographical  area,  under 
a  variety  of  weather  conditions.  To  the  extent  possible,  the  experi¬ 
mental  plot  shall  encompass  a  natural  variety  of  soil  types,  topo¬ 
graphic  features,  and  surface  characteristics  to  include  wooded 
areas,  lakes,  cultivated  land  (in  part  irrigated)  and  urban  areas. 
Experimental  measurements  shall  be  made  within  different  types  of 
air  masses  (including  Continental  Arctic,  Continental  Polar,  and 
Maritime  Tropical)  to  include  data  with  snow  cover  and  various  de¬ 
grees  of  solar  heating  and  vegetative  growth.  Instrumentation  shall 
be  both  ground  based  and  airborne.  Analyses  of  the  data  shall  be 
performed  in  order  to  evaluate  the  interrelationship  of  the  horizontal 
variation  of  micrometeorological  parameters  with  the  state  and  fea¬ 
tures  of  the  surface,  and  to  establish  mathematical  models  that  de¬ 
scribe  and  permit  the  prediction  of  the  variability  of  atmospheric 
boundary  layer  structure  in  a  quantitative  manner. " 

The  work  accomplished  during  the  first  two  years  under  the 
contract  has  been  devoted  to  the  following  subtasks: 
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Heinz  H.  Lettau 


1.  Ground-based  instrumentation 

2.  Airborne  instrumentation 

3.  Miaometeorological  measurements,  and 

4.  Evaluation  and  data  analysis 

Progress  has  been  made  in  all  subtasks.  A  list  of  titles  of  pub¬ 
lished  scientific  papers  follows  in  which  results  obtained  from  re¬ 
search  sponsored  by  the  contract  are  reported: 

J.  A.  Dutton,  "Space  and  Time  Response  of  Airborne  Sensors  for 
the  Measurement  of  Ground  Parameters:"  this  research  was  co¬ 
sponsored  by  ONR  under  Contract  Nonr  (1202)07  and  published  as  an 
ONR  Technical  Report  No.  1,  1959. 

H.  H.  Lettau,  "Wind  Profile.  Surface  Stress,  and  Geostrophic 
Drag  Co-efficients  in  the  Atmospheric  Surface  Layer,"  was  published 
by  Academic  Press,  Inc.,  New  York,  in  "Advances  in  Geophysics," 

Vol.  6,  1959. 

W.  L.  Pelton,  K.  M.  King,  and  C.  B.  Tanner,  "An  Evaluation  of 
the  Thornthwaite  and  Mean  Temperature  Methods  for  Determining 
Potential  Evapotranspiration, "  Agronomy  Journal,  V.oJ,  52,  p.  387-395, 
1960. 

W.  L.  Pelton  and  C.  B.  Tanner,  "Energy  Balance  Data  Hancock, 
Wisconsin,"  Soils  Bulletin  2,  University  of  Wisconsin,  June,  1960. 

W.  L.  Pelton  and  C.  B.  Tanner,  "Potential  Evapotranspiration  Esti¬ 
mates  by  the  Approximate  Energy  Balance  Method  of  Penman,"  Journal 
of  Geophysical  Research,  65,  No.  10,  p.  3391,  I960. 

C.  B.  Tanner  and  S.  M.  Robinson,  "Black-Body  Function  Sigma-T^' 
January,  1959,  University  of  Wisconsin,  Soils  Department  Report. 

H.  H.  Lettau,  "A  Theoretical  Model  of  Temperature  Variations  at 
the  Surface  of  an  Orbiting  Satellite,"  Journal  of  Geophysical  Research, 
66.  No.  11.  p.  3693-3698. 

C.  B.  Tanner,  "Energy  Balance  Approach  to  Evapotranspiration 
from  Crops,"  Proceedings,  Soil  Science  Society  of  America.  Vol.  24, 
p.  1-9,  January,  1960. 

C.  B.  Tanner,  "Radiant  Energy  Exchange  in  a  Cornfield,"  Agronomy 
Journal,  Vol.  52,  p.  373-379,  1960. 
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The  work  under  the  contract  has  been  advanced  on  a  broad  front. 
For  the  present  annual  report  it  was  decided  to  concentrate  on  evalu¬ 
ation  and  data  analysis,  and  to  discuss  subtasks  which  could  be 
attacked  with  a  minimum  of  required  new  instrument  development.  A 
report  on  instrumentation,  new  flight  investigations,  and  other  sub¬ 
tasks  will  be  submitted  in  due  time. 

The  main  characteristic  of  what  may  be  referred  to  as  the 
"Wisconsin-approach"  to  miaometeorological  research  is  the  empha¬ 
sis  on  (a)  controlled,  or  genuine  micrometeorological  experiments, 
and  (b)  the  relationship  between  small-scale  and  large,  or  synoptic- 
scale  meteorological  processes.  The  concept  of  controlled  micro- 
meteorological  experimentation  was  discussed  in  Section  2.2  of  the 
Final  Report,  Contract  DA-36-039-SC-80063,  entitled  "Research 
Problems  in  Micrometeorology,"  by  Heinz  H.  Lettau,  dated  September 
1959.  In  the  present  report,  a  first  series  of  thermal  response  ex¬ 
periments,  involving  artificial  cycling  of  natural  radiation  fluxes,  is 
desaibed  in  the  paper  by  J.  Lienesch,  and  a  first  series  of  wind- 
stress  experiments,  involving  a  controlled  modification  of  surface 
roughness,  in  another  paper  by  J.  Kutzbach.  The  results  of  these 
experiments  are  highly  interesting.  It  can  be  expected  that  they  will 
open  new  aspects  of  micrometeorological  research.  The  thermal  re¬ 
sponse  experiment  deals  with  a  problem -complex,  which  includes 
the  special  problem  of  heat  diffusion  in  a  non-homogeneous  cjonduc- 
tor  and  also  that  of  the  balance  between  heat  conduction  and  black- 
body  radiation,  under  extreme  conditions. 

Concerning  the  relationship  between  small-scale  and  large  to 
synoptic-scale  atmospheric  processes,  it  is,  perhaps,  surprising 
that  the  controlled  wind  stress  experiment  (reported  in  the  paper  by 
Kutzbach)  has  not  only  contributed  to  the  clarification  of  the  action 
of  an  ensemble  of  roughness  elements,  but  also  has  opened  aspects 
of  large-scale  application  of  a  typically  small-scale  experiment.  The 
key  for  this  is  the  use  of  the  concept  of  the  geostrophic  drag  coeffi¬ 
cient,  which  is  also  discussed,  and  directly  applied  to  synoptic- 
scale  problems,  even  numerical  weather  prediction  problems  (as  far 
as  dissipation  of  mechanical  energy  by  friction  is  concerned),  in 
Paper  No.  9.  The  other  papers  of  this  report  fit  into  the  frame  in  that 
they  deal  essentially  with  a  first  attempt  to  establish  the  regional 
and  continent-wide  basis  for  estimating  the  aerodynamic  roughness 
parameter  of  the  earth' s  surface.  This  is  an  important  step  in  the 
evaluation  of  the  interrelationship  of  the  horizontal  variation' of  , 
micrometeorological  parameters  with  the  state  and  features  of  the 
surface,  and  towards  establishment  of  mathematical  models  that  de¬ 
scribe,  and  permit  the  prediction  of  the  variability  of  atmospheric 
boundary  layer  structure  in  a  quantitative  manner. 
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Experimental  Investigation  of  the  Thermal  Response  of  the 
Air-Soil  System  to  Controlled  Radiation  Pulses^ 


James  H.  Lienesch 

Department  of  Meteorology 
University  of  Wisconsin 


Abstract.  The  results  of  an  experimental  investigation  of 
the  transformation  of  solar  energy  at  a  grass  surface  to 
sensible  heat  fluxes  into  air  and  soil  are  presented.  The 
dependence  of  the  response  of  the  air-soil  system  to  arbi¬ 
trary  cycling  of  incoming  solar  energy  by  means  of  a  mov¬ 
able  shading  device  is  examined  as  a  function  of  period. 
Comparison  with  Lettau'  s  theoretical  model  is  presented 
and  probable  causes  for  departure  from  this  model  are 
discussed.  ^ 

2. 1  Introduction 


The  incident  solar  radiation  upon  the  surface  of  the  earth  is 
transformed  into  heat  pulses  into  the  air  and  soil,  which  follow  the 
yearly  and  daily  variations  of  the  sun' s  intensity.  These  pulses 
are  harmonic  in  character  and  produce  sinusoidal  oscillations  in  the 
temperatures  of  the  two  mediums  (air  and  soil).  Lettau  (1951,  1952) 
has  presented  a  theory  which  describes  the  temperature  distribution 
in  terms  of  the  available  radiation  energy  and  the  thermal  parameters 
of  the  surface,  soil,  and  air.  He  has  described  the  response  of  a 
two-layer  system  to  periodic  heat  pulses.  Lettau  (i.959)  has  out¬ 
lined  an  experimental  method  for  the  determination  of  the  thermal 


*  Part  of  this  worlc  was  submitted  to  the  University  of  Wisconsin 
in  partial  fulfillment  of  the  requirements  for  the  degree  of  Master  of 
Science.  In  addition  to  sponsorship  by  the  U.  S.  Army  Electronic 
Proving  Ground,  Fort  Huachuca,  Meteorology  Department,  the  research 
in  this  paper  has  been  supported  by  the  National  Science  Foundation, 
Washington,  D.  C. ,  under  Grant  No.  G-5786. 
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response,  i.  e. ,  the  ratio  of  surface  temperature  amplitude  (deg)  and 
the  amplitude  of  the  forcing  function  (ly/sec),  for  a  variety  of  fre¬ 
quencies  of  the  forcing  function. 

In  the  spring  of  1959,  Dr.  Baumgartner,  University  of  Munich, 
was  visiting  the  University  of  Wisconsin.  He  prepared  a  study  of 
the  feasibility  of  this  experiment  and  computed  the  response  of  the 
two  mediums  to  short  period  energy  pulses,  on  the  basis  of  Lettau'  s 
model.  An  unpublished  manuscript  of  Baumgartner  (1959)  summarizes 
the  result. 

It  is  the  purpose  of  this  paper  to  report  the  results  of  the  exper¬ 
imental  investigation,  as  outlined  by  Lettau,  designed  to  explore  the 
process  Involved  in  the  transformation  of  a  periodic  energy  wave  at 
the  surface,  and  its  propagation  into  the  air  and  soil.  The  energy 
wave  was  obtained  by  periodically  shading  an  experimental  site  from 
direct  solar  radiation.  The  dependence  of  the  temperature  response 
in  the  two  mediums  on  the  period  of  energy  wave  can  be  observed  by 
varying  the  duration  of  the  controlled  radiation  pulses.  It  was  felt 
that  the  range  of  the  shading  periods  should  cover  an  order  of  magni¬ 
tude  to  obtain  significant  changes  in  the  response  to  the  forcing 
function. 

It  is  desirable  that  the  intensity  of  available  energy  be  nearly 
constant  throughout  the  entire  duration  of  the  experiment.  This  re¬ 
stricts  the  experiment  to  approximately  two  hours,  one  hour  before  to 
one  hour  after  noon  of  a  clear  day.  Since  the  longer  periods  of  the 
energy  wave  should  be  of  the  order  of  ten  times  that  of  the  shorter 
periods,  a  minimum  program  (within  two  hours)  consisting  of  a  suc¬ 
cession  of  two,  four,  ten,  twenty  and  forty  minute  periods  was  se¬ 
lected.  Time  allowed  the  immediate  repetition  of  the  shorter  periods 
so  that  a  typical  test  would  proceed  with  three  two-minute  periods, 
followed  by  two  four-minute  periods,  single  ten,  twenty,  and  forty 
minute  periods,  and  then  another  twenty,  ten,  etc.  ...  to  the  two- 
minute  periods  through  the  reverse  order.  For  an  illustration  of  this 
see  Fig.  1.  The  response  of  the  two  layer  system  could  then  be  ob¬ 
served  in  the  form  of  a  spectrum  for  the  frequencies  of  the  energy 
pulse. 

2.2  Description  of  Site 

The  site  of  the  experiment  was  an  open  field  near  the  flat  top  of 
a  low  hill,  part  of  the  Picnic  Point  Reservation  on  the  campus  of  the 
University  of  Wisconsin.  The  area,  sloping  slightly  (6%)  to  the 
east,  was  surrounded  by  farm  fields  with  the  nearest  stand  of  trees 
approximately  100  yards  distant  eastwards  and  down-hill  from  the 
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site.  The  test  field,  approximately  150  feet  square,  was  seeded 
with  grass  in  the  spring  of  the  year  in  order  to  provide  a  uniform 
surface  cover.  By  the  fall  of  the  year,  the  grass  was  quite  dense 
and  provided  a  relatively  smooth  surface. 

Near  the  center  of  the  field  a  radiation  shield  was  constructed 
to  provide  for  the  periodic  shading  of  the  experimental  area.  It 
was  made  of  aluminum  sheeting,  eight  feet  square.  It  could  be 
pivoted  about  a  horizontal  axis  at  its  lower  side  so  that  it  either 
shaded  the  experimental  area  or  exposed  the  area  to  the  direct 
solar  beam.  The  size  of  the  shaded  area  could  be  varied  by 
changing  the  angle  of  the  shield,  thus  minimizing  the  effect  of 
the  sun  angle.  Heating  of  the  shield  and  the  possible  resulting 
increase  in  long -wave  radiation  from  the  shield  was  reduced  by 
positioning  it  parallel  to  the  solar  beam  during  periods  of  exposure 
of  the  experimental  area^  to  direct  sunlight. 

The  area  affected  by  the  shading  was  mown  prior  to  each 
experiment  so  that  the  grass  height  was  about  two  inches.  The 
grass  in  the  surrounding  area  was  kept  at  a  height  of  about  three 
to  four  inches. 

2.3  Instrumentation 

2.3.1  Instruments  above  the  surface 

Incoming  short-wave  radiation  was  measured  by  an  Epply 
Pyrheliometer  at  a  height  of  approximately  one  meter.  This  in¬ 
strument  was  within  the  shaded  area,  and  exposed  so  that  during 
periods  of  shading  it  measured  the  diffuse  radiation  from  approx¬ 
imately  95%  of  the  sky.  This  pyrheliometer  had  a  time  lag  of 
about  one  minute,  therefore  instantaneous  values  of  short-wave 
radiation  immediately  following  shading  or  exposure  could  only 
be  estimated  from  the  values  indicated  after  the  instrument  came 
to  equilibrium. 

A  ventilated  net  radiometer  of  the  type  described  by  Suomi, 
Fransilla,  and  Islitzer  (1954)  was  used  to  measure  net  radiation. 
This  instrument  was  at  a  height  of  about  one  meter  and  positioned 
to  ’’see"  a  surface  area  immediately  adjacent  to  the  area  contain¬ 
ing  the  soil  instruments.  The  surface  beneath  the  net  radiometer 
was  equivalent  to  that  beneath  which  the  soil  sensors  were  lo¬ 
cated.  The  time  lag  of  this  instrument  was  only  a  few  seconds. 

It  also  was  within  the  area  of  shading. 

The  most  difficult  problem  proved  to  be  an  accurate  method  of 
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determining  air  temperature  close  to  the  grass  surface.  It  would  be 
desirable  to  employ  an  instrument  which  is  completely  free  of  any 
radiation  error.  Since  such  an  ideal  sensor  was  not  available,  an 
attempt  was  made  to  determine  and  eliminate  the  radiation  error  of  - 
the  sensor.  One  possible  method  is  to  paint  one  group  of  thermo¬ 
couples  white  and  a  separate  group,  black.  By  varying  the  number 
of  white  thermocouples  relative  to  the  black  thermocouples  and 
opposing  or  "bucking"  the  voltages  produced,  it  was  hoped  that 
radiation  error  would  be  eliminated.  It  was  found  that  such  an  in¬ 
strument,  upon  exposure  to  radiation,  had  a  lag  time  which  was 
longer  than  could  be  tolerated  in  this  Investigation.  Specifically, 
this  sensor  system  tended  to  produce  an  initial  overshoot  in  voltage 
when  solair  radiation  was,  first  incident  and  an  initial  overshoot  when 
it  was  removed. 

It  seemed  therefore  more  advisable  to  record  the  air  temperatures 
with  the  aid  of  two  separate  thermocouples,  one  painted  black  and 
the  other  white.  Knowing  the  absorptivity  of  the  two  paints  (0.31  for 
white  and  0.88  for  black),  the  extrapolation  to  zero  absorptivity  (i.e., 
elimination  of  radiation  error)  was  done  numerically.  Two  white 
thermocouples  and  two  black  thermocouples  were  alternately  arranged 
on  a  support  which  held  the  sensors  about  two  inches  above  the  soil 
surface. 

Average  wind  speed  for  the  duration  of  a  test  was  measured  by  a 
cup  anemometer  at  a  height  of  approximately  two  feet.  A  hot-wire 
anemometer  was  used  on  occasion  (days  with  large  wind  velocities) 
to  measure  maximum  wind  velocities  prior  to  the  test.  These  meas¬ 
urements  were  the  basis  for  the  estimate  of  the  diffusivity  of  the  air. 

2.3.2  Instruments  below  the  surface 

To  measure  the  soil  temperature  a  thermopile  consisting  of  eight 
thermocouples  was  placed  in  the  soil  at  a  depth  of  0. 5  cm.  The  eight 
thermocouples  were  located  within  a  circle  about  3  inches  in  diameter 
to  provide  an  integrated  value  over  a  small  area.  It  was  not  an  easy 
task  to  position  each  individual  junction  so  that  a  0.5  cm  layer  of 
soil  covered  each  one.  However,  it  is  felt  that  the  average  soil  cover 
was  satisfactorily  close  to  0.50  cm. 

The  difference  in  soil  temperature  between  the  0.5  cm  level  and 
the  2. 5  cm  level  was  measured  directly  by  a  thermopile  with  one  set 
of  ten  junctions  at  each  level.  The  two  sets  of  junctions  were  placed 
in  a  horizontal  line  approximately  2  inches  in  length.  Again  it  was 
difficult  to  insure  that  the  proper  depths  were  maintained  individually 
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but  the  average  levels  for  each  of  ^e  two  sets  of  junctions  should 
be  satisfactorily  close  to  the  assumed  depths  (0.50  cm  and  2.5  cm). 

All  the  thermocouples  were  made  from  No.  30  Leeds  and  Northrup 
thermocouple  wire.  Those  placed  in  the  soil  were  covered  by  a  plas¬ 
tic  jacket  to  prevent  electrical  conduction  between  sensors  through 
the  moist  soil.  The  remainder  of  the  thermocouple  wire  in  contact 
with  the  soil  was  coated  with  Dow  Coming  808  Resin  and  a  varnish 
to  minimize  electrical  leakage  to  the  soil. 

At  a  depth  of  1.5  cm  a  heat  flux  transducer  was  placed.  The 
transducer  itself  is  composed  of  a  silver -cons tantan  thermopile  ar¬ 
ranged  in  a  thin  phenolic  resin  plate.  This  instrument,  manufactured 
by  Beckman  and  Whitley  Company,  indicates  the  heat  flow  through 
the  plate  by  an  electromotive  force  produced  by  the  difference  in  tem¬ 
perature  between  the  two  sets  of  thermoj unctions,  which  are  parallel 
and  adjacent  to  the  two  sides  of  the  plate.  The  sensitivity  of  the 
transducer  is  temperature  dependent.  Therefore,  a  thermocouple 
junction  in  the  plate  provides  for  measurement  of  the  transducer 
temperature,  which  corresponds  to  an  additional  soil  temperature 
measurement  at  a  depth  of  1.5  cm. 

The  various  soil  sensors  were  installed  in  the  soil  by  digging 
three  shallow  holes,  each  with  one  vertical  wall.  The  sensors  were 
then  placed  in  the  vertical  walls  at  their  respective  depths,  forming 
a  semicircular  airangement. 

As  a  reference  junction  for  all  thermocouples  a  thermos  bottle 
was  placed  below  the  soil  surface  at  a  horizontal  distance  of  approx¬ 
imately  1.5  feet  from  the  sensors.  Oil  was  used  as  a  reference  bath 
because  of  its  low  electrical  conductivity.  The  temperature  of  the 
oil  was  read  by  a  mercury  thermometer  accurate  to  the  nearest  tenth 
of  a  degree.  Each  reference  junction  in  the  bath  was  covered  with  a 
plastic  jacket  to  prevent  contact  between  junctions. 

2.3.3  Recording  and  calibration 

Tentative  experiments  using  a  different  shading  device  were  made 
by  Dr.  Baumgartner  at  the  same  field  site  in  May,  1959.  At  that  time 
a  one-channel  recorder  was  used,  with  switching  devices.  It  became 
immediately  obvious  that  continuous  and  simultaneous  recording  of 
all  elements  on  one  chart  is  a  necessity. 

For  the  series  of  experiments  described  here,  a  Honeywell  906B 
Visicorder  Oscillograph  was  used.  Through  mlnror  galvonometers 


Table  1.  Results  of  thermal  response  experiment  at  the  University  of  Wisconsin  Micrometeorological  Site  during  three  clear 
days  in  Fall  of  I960.  Listed  are  double  amplitudes  of:  (1)  air  temperature  at  5  cm,  (2)  soil  temperature  at  0.5  cm. 
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this  Instrument  can  record  up  to  12  signals  simultaneously  and  con¬ 
tinuously,  without  interference,  on  photographic  paper  which  develops 
upon  exposure  to  ultraviolet  light.  A  total  of  eight  galvonometers  was 
used,  one  for  each  instrument:  ( 1 )  the  Epply  pyrheliometer,  (2)  the 
net  radiometer,  (3)  the  white  air  temperature  thermocouples,  (4)  the 
black  air  temperature  thermocouples,  (5)  the  soil  temperature  thermo¬ 
pile  at  0.5  cm,  (6)  the  thermopile  for  measuring  the  difference  in  soil 
temperature  between  the  0.5  cm  level  and  the  2.5  cm  level,  (7)  the 
soil  heat  flux  transducer  at  1.5  cm,  and  (8)  the  temperature  of  the  soil 
heat  flux  transducer.  Trace  Identification  was  accomplished  through 
occasional  systematic  shorting  of  the  circuits  in  a  prescribed  se¬ 
quence.  The  paper  transport  of  the  available  instrument  was  incon¬ 
veniently  fast.  A  timer  was  constructed  which  allowed  the  paper 
drive  mechanism  to  operate  for  only  one  second  in  every  ten  seconds. 
Thus  the  record  obtained  indicated  the  voltage  corresponding  to  a 
simultaneous  sampling  rate  of  once  every  ten  seconds. 

The  recorder  was  located  in  a  house  trailer  approximately  80  feet 
from  the  experimental  area  and  connected  to  three  8 -conductor  mag¬ 
netically  insulated  cables  in  a  plastic  tube  6  inches  below  the  soil 
surface.  A  junction-box  10  feet  from  the  sensors  was  used  to  connect 
the  cables  to  the  2 -conductor  cables  from  each  sensor.  Grounding 
was  provided  at  the  trailer  housing  the  recorder. 

The  Epply  pyrheliometer  and  the  net  radiometer  were  calibrated  in 
early  fall  of  1960  with  the  use  of  a  secondary  standard,  an  Abbot 
Silver-Disk  pyrheliometer.  The  calibrations  of  these  instruments  on 
two  separate  days  were  within  4%  for  the  net  radiometer,  and  2%  for 
the  pyrheliometer. 

The  calibration  for  the  heat  flux  transducer  was  provided  by  the 
manufacturer.  The  calibration  for  the  thermocouples  was  obtained 
from  the  Standard  Conversion  Tables  for  Leeds  and  Northrup  Thermo¬ 
couples. 

\ 

2.4  Experimental  Results 


2.4.1  Raw  data 

Due  to  the  nature  of  the  experiment  it  was  necessary  to  confine 
the  investigation  to  cloudless  days  with  the  wind  direction  between 
southwest  and  northwest.  After  completion  of  the  instrumentation,  and 
before  the  winter  season,  only  three  such  days  occurred,  on  October  17, 
October  24,  and  November  3,  1960.  On  October  24  the  wind  was  very 
light;  on  October  17,  gusty  to  20  mph;  and  on  November  3,  gusty  to 
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25  mph.  Owing  to  rainy  days  in  between,  the  soil  moisture  condi¬ 
tions  also  were  different,  as  is  discussed  in  Section  2.6.  2,  in  con¬ 
nection  with,  soil  heat  conductivity  data. 

The  recorded  values  were  read  from  the  photographic  paper  and 
were  re  plotted  for  each  of  the  three  days.  Trends  were  then  elimi¬ 
nated  from  the  results.  This  was  done  in  the  following  manner.  En¬ 
velopes  encompassing  the  maximum  and  minimum  values  of  ampli¬ 
tudes  for  the  sequence  of  radiation  pulses  were  constructed.  The 
actual  data  was  then  expressed  as  the  numerical  difference  between 
the  two  envelope  curves. 

Records  for  the  three  days  did  not  show  any  great  dissimilarities. 
The  amplitude  values  were  not  constant  because  of  changes  in  solar 
radiation,  air  and  soil  parameters  from  one  day  to  the  other.  But,  in 
general,  the  character  of  the  amplitude  response  to  the  periodic  heat 
pulses  as  a  function  of  period  was  found  to  be  satisfactorily  repro¬ 
ducible. 

-  ^ 

Table  1  gives  the  response  of  the  measured  elements  to  the 
various  frequencies  of  the  forcing  functions  for  the  three  days.  The 
values  presented  in  this  table  are  equivalent  to  twice  the  amplitude 
of  the  response,  i.  e. ,  they  are  the  difference  between  maximum  and 
minimum  points  occurring  during  the  sequence  of  radiation  pulses. 

2.4.2  Analysis  of  October  17  data 

The  data  obtained  on  October  17  is  fairly  representative  of  the 
data  obtained  on  the  remaining  two  days.  Figure  1  is  a  composite 
showing  the  response  of  the  system  to  the  controlled  pulses  forced 
upon  the  experimental  area  on  October  17.  The  following  discussion 
applies,  in  general,  also  for  the  data  obtained  on  October  24  and 
November  3. 

The  values  for  solar  short-wave  radiation  and  for  the  total  net 
radiation  give  an  indication  of  the  amount  of  energy  potentially  avail¬ 
able  at  the  surface  for  transformation  into  sensible  heat.  These 
curves  also  illustrate  the  controlled  cycling  of  the  radiation  pulses. 
The  relatively  small  over-all  variation  in  the  solar  radiation  was 
insignificant  for  the  results  of  the  experiment.  The  decrease  in  the 
amount  of  net  radiation  at  the  end  of  the  test  period  is  explained  by 
the  gradual  march  of  the  shaded  area  across  the  experimental  site  as 
a  result  of  the  changing  sun  angle.  The  net  radiometer  "saw"  this 
gradual  advance  of  unshaded  surface  because  of  its  relatively  large 
field  of  view. 
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In  contrast  to  the  soil  variables  the  amplitudes  of  the  air  temper¬ 
ature  show  a  relatively  weak  dependency  on  the  frequency  of  the 
forcing  function.  The  method  by  which  the  air  temperature  was  ob¬ 
tained  was  discussed  in  Section  2.3. 1.  It  amounts  to  the  elimination 
of  the  radiation  error  with  the  aid  of  two  sensors  having  known  but 
relatively  large  radiation  errors.  It  was  thought  that  possibly  the  ab- 
sorptivities  assumed  for  the  two  paints  were  not  accurately  enough 
known.  Therefore,  extrapolation  to  zero  absorptivity  was  tested  as¬ 
suming  ±10%  variation  in  the  absorptivities  of  the  two  paints.  No 
significant  change  in  the  result  was  obtained.  Convective  cooling 
was  apparent  as  is  illustrated  by  the  dispersion  of  the  air  temperature 
data  points  in  Fig.  1.  However,  it  is  not  felt  that  the  cooling  of  the 
sensors  by  the  wind  was  mainly  responsible  for  the  result  obtained. 

The  reliability  of  the  air  temperature  measurement  in  this  experiment 
is  still  an  open  question. 

The  curves  representing  the  soil  temperature  and  the  difference 
in  temperature  between  the  0.5  cm  level  and  the  2.5  cm  level  follow 
the  forcing  function  with  relatively  short  phase  lag,  but  with  a  char¬ 
acteristic  transformation  from  the  originating  rectangular  step-function 
to  sinusoidal  shape.  In  view  of  this  transition,  phase  lags  will  not 
be  discussed.  The  emphasis  will  be  on  amplitude  relations.  At  the 
beginning  of  the  test,  the  general  cooling  trend  shown  by  the  soil 
temperature  is  hardly  affected  by  tht  3olar  radiation  pulses.  When 
the  duration  of  shading  reached  two  minutes  or  more,  the  effect  was 
great  enough  to  overcome  the  initial  cooling  trend.  At  the  end  of  the 
test,  the  variations  in  soil  temperature  to  the  shorter  pulses  are  quite 
recognizable. 

The  curve  representing  the  temperature  difference  between  the 
0.5  cm  and  2.5  cm  s^oil  level  shows  the  response  to  all  pulses  of  solar 
radiation.  For  the  20-minute  heating,  the  AT-amplitude  is  not  as 
large  as  might  be  expected  because  the  heat  flux  reaches  the  lower 
set  of  junctions,  thus  diminishing  the  temperature  difference  between 
the  two  levels.  At  the  end  of  the  20-mihute  shading  (at  absolute  time 
equal  to  82  minutes),  the  figure  shows  that  the  temperature  at  2.5  cm 
is  slightly  larger  than  that  at  0.5  cm,  indicating  a  significant  heat 
loss  in  the  layer  of  soil  close  to  the  surface. 

The  curve  representing  the  heat  flux  in  the  soil  resembles  the 
AT-curve  quite  well.  The  flux  plate  was  inoperative  until  minute  52 
of  the  test.  This  instrument  recorded  an  outward  flux  at  the  end  of 
both  the  20-minute  and  10-minute  shadings.  The  difference  in  temper¬ 
ature  between  the  0.5  cm  and  2.5  cm  levels  does  not  show  this  reversal 
for  the  10-minute  shading,  but  this  is  tolerable  considering  the  differ¬ 
ence  between  measurement  at  the  depth  of  the  flux  plate  (1.5  cm)  and 


Thermal  Response  to  Controlled  Pulses 


15 


a  mean-layer  depth  between  levels  of  0.5  and  2.5  cm.  The  shadings 
of  duration  below  10  minutes  did  not  permit  heat  flux  toward  the  sur¬ 
face  at  the  1.5  cm  level. 

2.5  Discussion  of  Errors 


The  method  of  recording  the  simultaneous  response  of  the  various 
elements  required  that  the  data  be  read  from  the  multichannel  visi- 
corder  chart.  '-A  certain  Amount  of  error  is  present  in  reading  the  paper 
trace.  Also  there  is  the  tendency  to  manufacture  deviations  in  the 
trace  when  such  are  desired.  However,  it  is  felt  that  such  errors  are 
small  and  in  no  way  invalidate  the  result  of  the  experiment;  with  the 
possible  exception  of  the  one-  and  two-minute  shading  durations. 

The  largest  error  persistent  throughout  the  experiment  is  undoubt¬ 
edly  in  the  measurement  of  air  temperatures.  The  presence  of  natural 
fluctuations  of  wind  produced  unsteady  signals  from  the  two  types  of 
sensors  measuring  air  temperature. 

In  general,  the  results  of  this  experiment  are  not  so  dependent  on 
the  absolute  values  obtained  but  rather  on  the  relative  responses  of 
all  the  instruments  to  the  same  forcing  functions. 

2.6  Data  Evaluation  and  Comparison  with  Theory 

2.6.1  Albedo  estimates 

The  simultaneous  measurements  of  solar  radiation  and  net  radia¬ 
tion  provide  a  basis  for  estimates  of  the  surface  albedo.  The  equation 
for  net  radiation  permits  us  to  express  the  albedo  [a)  of  the  surface  as 
follows: 


AN  =  AS(l-«);  or,  a  =  1  -  AN/AS  ,  (1) 

where  AN  is  the  amplitude  of  the  net  radiation  (i.  e. ,  the  difference 
in  net  radiation  between  shaded  and  exposed  conditions)  and  AS  is 
the  amplitude  of  the  solar  radiation.  The  values  of  albedo  for  the 
three  days  are  not  the  same,  although  the  values  computed  at  different 
times  for  the  individual  days  are  very  consistent.  A  relationship 
seems  to  exist  between  the  albedo  and  the  moisture  content  of  the 
soil.  On  October  24,  the  soil  contained  the  smallest  amount  of,  mois¬ 
ture  and  the  corresponding  mean  albedo  and  its  standard  deviation  was 
23  ±  1.6%  .  On  October  17,  with  larger  moisture  content,  the  albedo 
was  15  ±  5; 8%  and  on  November  3,  with  the  largest  moisture  content, 

9  ±  2.3%  .  It  is  reported  in  the  literature  that  the  albedo  of  various 
soils  decreases  with  the  moisture  content. 
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2.6.2  Determination  of  soil  heat  conductivity 

The  heat  conductivity  (ly/min  per  deg/cm)  of  the  soil  was  ob¬ 
tained  graphically  by  plotting  the  instantaneous  values  of  soil  heat 
flux  (ly/mln)  at  the  1.5  cm  level  versus  half  the  difference  in  temper¬ 
ature  between  the  0.5  cm  and  2.5  cm  levels.  The  heat  conductivity 
(\)  is  defined  in  the  equation 

\  =  -F/T'  ,  or  F  =  -XT'  (2) 

where  F  is  the  heat  flux,  and  T'  the  vertical  temperature  gradient 
("C/cm).  Fig.  2  shows  the  data  points  obtained  from  shading  experi¬ 
ments  on  the  three  different  days.  The  slope  of  the  lines  of  best  fit 
result  in  the  following  X  values:  0.00122  for  October  24,  0.00146 
for  October  17,  and  0.00191  for  November  3,  where,  for  convenience, 
the  X  units  are  ly/sec  per  deg/cm. 

The  soil  at  the  site  was  identified  as  a  silt  loam.  Carter  (1951) 
has  found  that  X  for  various  silt  loams  ranges  from  0.0006  for  dry 
conditions  to  0.0040  for  wet  conditions  (30%  moisture).  This  indi¬ 
cates  that  the  soil  heat  conductivity  increases  with  increasing  mois¬ 
ture  content.  The  work  of  Patten  (1909)  substantiates  this  conclusion. 

Table  2  shows  the  relationship  between  the  observed  \  and  the 
total  amount  of  rainfall  for  the  period  of  seven  days  preceding  each 
test  day.  The  actual  soil  moisture  content  was  not  determined  at  the 
time  of  the  test;  however,  this  table  provides  an  indication  of  the 
amount  of  moisture  present  in  the  soil.  An  increase  of  the  X  value 
with  increasing  soil  moisture  is  concluded. 


Table  2.  Soil  heat  conductivity  and  total  rainfall  for  the  seven  days 
preceding  the  dates  of  the  experiment. 


Date 


Rainfall 


Conductivity 


October  24 

0.05 

October  17 

0,25 

November  3 

2.50 

inches 

0.00 

inches 

0.00 

inches 

0.00 

122  ly/sec  per  deg/cm 
146  ly/sec  per  deg/cm 
191  ly/sec  per  deg/cm 


2.6.3  Extrapolation  of  soil  response  to  the  surface 

Since  there  was  no  independent  measurement  of  the  temperature 
of  the  soil  surface,  it  became  necessary  to  estimate  this  value  from 
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the  observed  temperature  at  the  0.5  cm  depth.  In  addition,  the  value 
of  the  soil  heat  flux  at  the  surface  is  also  required  and  may  be  calcu¬ 
lated  knowing  the  heat  flux  at  1.5  cm.  For  simplicity,  the  well-known 
classical  solution  of  the  thermal  diffusion  equation  for  homogeneous 
soil  will  be  used.  The  equation  for  the  two  above  computations  is 

log  (AX,/^X  )  =  z  •J  n/2K  (3a) 

e  z 

where  AXq  is  the  amplitude  (of  temperature  or  flux)  at  the  surface, 
AXz  is  the  amplitude  (of  temperature  or  flux)  at  the  level  z,  n  is 
frequency  =  2ir/T  where  T  is  the  period,  K  is  thermal  diffuslvity 
=  K/c(cm^  /sec),  c  is  volumetric  heat  capacity  (cal/deg  per  cmM* 

An  average  value  of  c  =  0.4  has  been  assumed,  and  \- values 
for  the  individual  days  as  derived  from  Fig.  2. 

To  determine  the  soil  temperature  at  the  surface,  equation  (3a) 
becomes 

log  (AToATcs)  =  0.5><fn/2K  (3b) 

For  the  heat  flux  at  the  surface,  equation  (3a)  gives 

log^(AFg/AF,.5)  =  1.5'yv^  (4) 

where  AT0.5  A  F,.5  are  the  observed  amplitudes. 

Figure  3  illustrates  how  the  amplitude  of  the  heat  flux  in  the  soil 
at  1.5  cm  varies  with  the  frequencies  of  the  radiation  pulse.  The  am¬ 
plitude  of  the  soil  heat  flux  (AFo)  has  been  divided  by  the  value  of 
the  net  radiation  amplitude  (AN)  in  order  to  facilitate  a  comparison  of 
data  for  the  three  days.  The  data  presented  in  Table  1  wasvused  for 
the  evaluation  of  AFq  described  above. 

The  ratio  of  the  air  temperature  amplitudes  to  the  estimated  soil 
surface  temperature  amplitudes  is  shown  in  Figure  4,  which  indicates 
that  this  ratio  is  very  large,  especially  for  the  energy  pulses  of 
higher  frequency.  This  result  is  due  to  the  fact  that  the  air  tempera¬ 
ture  amplitude  was  found  nearly  independent  of  the  frequency  of  the 
forcing  function  as  shown  in  Fig.  1.  Possible  explanations  for  this 
result  are  presented  in  the  following  section. 

2.6.4  Comparison  with  theory 

A  theoretical  model  for  the  response  of  the  two-layer  system 
(soil-air)  to  harmonic  energy  waves  was  discussed  by  Lettau  (1951). 


Figure  4 
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He  has  expressed  the  soil  surface  temperature  amplitude  as  a  func¬ 
tion  of  frequency,  and  the  amplitude  (uq  )  of  the  energy  available  for 
transformation  to  sensible  heat  only.  If  this  model  is  applied  to  nat¬ 
ural  soil,  the  net  radiation  amplitude  must  be  reduced  by  the  portion 
of  energy  used  for  evaporation.  In  Figure  5,  the  three  theoretical 
curves  illustrate  the  soil  surface  temperature  amplitude  as  a  function 
of  the  period  of  the  forcing  function.  These  curves  are  labeled  with 
the  values  assumed  for  the  soil  heat  conductivity  (X)  and  adiabatic 
mixing  velocity  (w*)  which  measures,  in  Lettau'  s  model,  the  diffu- 
sivity  of  the  air. 

Another  parameter  in  Lettau 's  model  is  the  aerodynamic  surface 
roughness  (Zq).  A  constant  value  of  Zq  =  0.5  cm  was  used.  The 
values  for  X  conrespond  to  existing  soil  conditions  for  the  three 
days;  the  values  of  vr*  are  estimated  from  measurements  of  the  wind 
speed  at  two  feet  above  the  ground  and  the  maximum  wind  speeds  de¬ 
termined  prior  to  the  tests  by  the  hot-wire  anemometer.  Within  toler¬ 
able  error  limits,  the  curve  labeled  w*^  =  10  cm/sec  applies  to 
October  24,  w*  =  20  cm/sec  to  October  17,  and  w*  =  40  cm/sec 
to  November  3. 

Fig.  5  also  shows  the  experimental  results  obtained  on  the  three 
days,  in  a  plot  of  the  amplitude  ratio  of  soil  surface  temperature  to 
net  radiation.  Note  should  be  made  of  the  different  scales  used  for 
the  ordinates  of  experimental  and  theoretical  curves.  Within  tolerable 
error  limits  the  general  trend  of  the  experimental  results  is  similar  to 
that  predicted  by  the  theory.  Consequently  the  discrepancy  in  the 
magnitudes  between  the  predicted  and  experimental  values  can  be  at¬ 
tributed  to  the  evaporation  process  which  is  not  considered  in  the  the¬ 
oretical  model.  A  quantitative  discussion  of  evaporation  effects  is 
presented  in  the  following  section. 

The  theoretical  amplitude  of  soil  heat  flux  at  the  surface  relative 
to  that  of  available  energy  for  conduction  of  sensible  heat  (AFq/uo) 
is  shown  in  Figure  6.  The  upper  curve  is  computed  for  soil  heat  con¬ 
ductivity  equal  to  0.00122  (ly/sec  per  deg/cm)  and  an  adiabatic 
mixing  velocity  of  10  cm/'sec,  i.  e. ,  the  values  corresponding  to  con¬ 
ditions  of  October  24;  the  other,  for  X  =  0.00191  and  w*  =  40  corre¬ 
sponding  to  November  3.  The  curve  for  October  17,  for  X  =  0.00146 
and  w*  =  20 ,  would  lie  between  these  two. 

Fig.  6  illustrates  the  difference  between  the  theoretical  values 
and  those  obtained  experimentally  on  the  three  days.  The  amplitude 
ratio  of  soil  heat  flux  at  the  surface  to  net  radiation,  as  determined 
experimentally,  is  in  general  approximately  20%.  However,  for  the 
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assumed  soil  and  air  parameters  for  the  three  days  and  the  frequen* 
cies  of  the  forcing  function  used  in  the  experiment,  the  theory  pre¬ 
dicts  that  more  than  70%  of  the  energy  available  for  sensible  heat 
transfer  will  be  directed  into  the  soil.  The  disaepancy  between  the 
values  determined  experimentally  and  those  predicted  by  theory  are 
attributed  here  also  to  the  evaporation  process  occurring  at  the  soil 
surface. 

2.7  Conclusions 

2.7.1  Apparent  deviation  from  two-layer  theory 

Lettau'  s  theoretical  model  assumes  that  the  lower  medium  (soil) 
is  homogeneous,  separated  from  the  upper  medium  (air)  by  a  horizontal 
surface  of  discontini|ity.  Possibly  the  site  upon  which  the  experiment 
was  performed  did  not  represent  satisfactorily  a  simple  two-layer  sys¬ 
tem.  With  a  grass  height  of  approximately  two  Inches,  an  additional 
layer  might  be  postulated  between  the  soil  and  air.  The  presence  of 
a  third  layer,  consisting  of  a  mat  of  grass,  quite  possibly  could  have 
contributed  to  the  unexpectedly  large  air  temperature  amplitudes,  and 
their  small  variation  with  the  period  of  the  incident  radiation. 

Owing  to  the  late  season  (low  radiation  intensities),  preliminary 
experiments  performed  on  a  grass-free  surface  were  inconclusive  in 
establishing  a  better  understanding  of  the  processes  involved.  On  the 
other  hand,  it  can  be  doubted  that  the  method  of  elimination  of  radia¬ 
tion  errors  from  the  air  temperature  sensors  was  satisfactory.  A  more 
refined  technique  for  the  measurement  of  air  temperature  will  supply 
more  dependable  data  to  resolve  the  existing  uncertainties. 

2.7.2  Evidence  of  evaporation  processes 

In  Section  2.6.4  the  discrepancy  between  the  theoretical  and  ex¬ 
perimental  values  shown  in  Fig.  5  was  attributed  to  the  evaporation 
process.  The  ratio  existing  between  the  two  curves  for  individual 
days  indicates  that  only  approximately  20%  of  the  net  radiation  has 
been  converted  to  sensible  heat.  This  suggests  that  80%  of  the  net 
radiation  energy  is  being  used  for  evaporation.  In  view  of  ample  soil 
moisture  at  the  test  site,  this  ratio  agrees  fairly  well  with  that  ob¬ 
tained  on  occasion  (days  with  relatively  high,  but  absolutely  rather 
moderate,  soil  moisture)  by  Suomi  and  Halstead  at  O'Neill,  Nebraska. 
Reference  is  made  to  H.  Lettau  and  B.  Davidson  (1957),  Vol.  II, 

Tables  7. 3. a. ,  see  Table  3.  Tanner  and  Pelton  (1960)  found  that 
evaporation  from  a  good  alfalfa-brome  cover  may  be  greater  than  or 
equal  to  90%  of  the  net  radiation.  With  advective  heating  they  found 
that  evaporation  may  actually  exceed  net  radiation. 


AMPLITUDE  RATIO  OF  SOIL  HEAT  FLUX  AT  2*0  (AFq.  LY/MIN) 


SHADING  PERIOD  (MINUTES) 


Thermal  Responses  to  Controlled  Pulses 
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Table  3.  Net  radiation  (No)  and  evaporation  (Eq)  in  millily/min, 
averaged  for  hours  10  A.  M.  to  2  P.  M.  at  O'  Neill,  Nebraska,  ac* 
cording  to  Suomi  and  Halstead. 

_ Suomi _  Halstead 

Date  N,  E,  100E,Al,  N«  E,  100  E,/N, 

August  9,  1953  746  397  53  826  464  56 

August  14,  1953  -780  449  58  857  543  64 

The  estimate  of  evaporation  following  the  original  approach  out¬ 
lined  by  Lettau  (1952)  is  a  promising  aspect  of  the  thermal  response 
experiment.  The  data  in  Fig.  6  confirm  the  result  that  approximately 
80%  of  the  net  radiation  energy  was  used  for  evaporation. 

The  understanding  of  the  evaporation  process,  too,  requires  an 
investigation  concerning  the  validity  of  the  two-layer  system  dis¬ 
cussed  above.  Postulating  a  third  layer  of  grass  greatly  increases 
both  the  surface  area  for  exchange  of  sensible  heat  and  evaporation. 
Further  experiments  are  necessary  to  define  the  processes  occurring 
over  a  grass  surface. 

2.8  Recommendations 

To  perform  an  experiment  aimed  at  investigating  the  two-layer 
conduction  processes,  it  will  be  necessary  to  insure  that  truly  two 
mediums  only  are  present.  Therefore,  such  an  investigation  will  be 
best  performed  on  a  soil  surface  which  is  free  of  vegetation,  thus 
eliminating  a  possible  third  or  buffer  layer.  To  provide  data  for  com¬ 
parison  with  the  theory  of  two-layer  heat  conduction,  a  further  re¬ 
striction  is  advisable.  As  reported  in  the  body  of  this  report  the 
process  of  evaporation  apparently  involves  a  large  fraction  of  net 
radiation,  not  considered  in  the  theoretical  study.  Evaporation  is 
eliminated  by  performing  a  similar  experiment  on  a  moisture-free 
medium.  Dry  concrete  or  asphalt  would  conveniently  provide  such  a 
surface. 

A  more  detailed  investigation  of  the  evaporation  processes  as  re¬ 
lated  to  the  period  of  the  energy  wave  would  be  interesting.  Compar¬ 
isons  with  direct  methods  of  determining  the  evaporation,  perhaps, 
performing  shading  experiments  on  lysimeters  (a  sensitive  weighing 
system)  or  one  permitting  active  control  of  the  moisture  flux  are  prom¬ 
ising  possibilities. 
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Further  investigations  similar  to  that  reported  here  should  incor¬ 
porate  an  independent  measurement  of  soil  surface  temperature,  either 
by  contact  thermometers,  or  bolometers.  In  addition,  the  measure¬ 
ment  of  air  temperature  without  radiation  error  is  a  necessity.  An 
aspirated  thermocouple  surrounded  by  a  high  reflecting  cylinder 
might  be  useful  but  has  the  undesirable  features  of  casting  a  shadow 
upon  the  underlying  surface  and  disturbing  the  air  by  ventilation. 

Experiments  along  the  lines  discussed  here  could  add  greatly  to 
an  understanding  of  the  transformation  of  energy  occurring  at  the  sur¬ 
face  of  the  earth  which  is  one  of  the  most  important  micrometeorolog- 
ical  processes. 
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Derivation  of  Roughness  Parameters  from  Wind  Profile  Data 
Above  Tall  Vegetation 

Ernest  Kung 

Department  of  Meteorology 
University  of  Wisconsin 


Abstract.  Wind  profile  measurement  above  tall  vegetation 
(including  corn,  brush  fields,  orchards,  and  fir  forest)  re¬ 
ported  in  the  literature  are  reanalyzed.  The  logarithm  of 
the  roughness  parameter  Zq  is  found  to  inaease  in  propor¬ 
tion  to  the  logarithm  of  vegetation  height;  however,  the 
probable  error  of  this  regression  line  is  relatively  large. 

The  uncertainties  of  the  analysis  are  discussed,  and  ex¬ 
plained  by  certain  unreliabilities  of  anemometry. 

3.1  Introduction 

The  roughness  length  (zq,  cm),  the  zeropoint-displacement 
(d,  cm),  and  the  surface  stress  (tq,  dynes/cm^)  are  the  three  wind 
profile  parameters  in  the  study  of  the  atmosphere  near  the  ground.  In 
specifying  the  degree  of  roughness  of  a  given  surface,  we  must  con¬ 
sider  at  least  two  factors:  the  height  of  protuberance,  and  spacing 
between  individual  elements.  Thus,  the  wind  profile  parameter  will 
be  determined  by  the  structure  and  nature  of  vegetation  cover  at  the 
surface. 

Deacon  (1953)  has  quoted  roughness  parameter  values  for  various 
natural  surfaces  which  are  summarized  in  Table  1.  His  study  referred 
to  surfaces  without  natural  vegetation  and  surfaces  of  low  vegetation. 
In  this  study  it  was  attempted,  through  the  analysis  of  wind  profile 
data  reported  in  the  literature,  to  extend  Deacon' s  tabulation  to  sur¬ 
faces  of  various  kinds  of  relatively  tall  vegetation. 
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Table  1.  Roughness  parameters  of  various  surfaces  adopted  from 
Deacon  (1953). 

Type  of  Surface  Zo  (cm) 

Smooth  mud  flats  0.  001 

Smooth  snow  on  short  grass  0. 005 

Sea  0. 02 

Desert  near  Karachi  0.  03 

Snow  surface^  Canada  0. 10 

Mown  grass  surfaces 

Grass  length  1.  5  cm  0.  2 

Grass  length  3.  0  cm  ..  0.  7 

Grass  length  45  cm:  V2  =  2  m/sec  2.  4 

V2  =  6  to  8  m/sec  1.  7 

Long  grass,  60  to  70  cm  high: 

V2  =  1.  5  m/sec  9.  0 

V2  =  3.  5  m/sec  6.  1 

V2  =  6.  2  m/sec  3.  7 

Downland  at  Porton:  Winter  1  to  2 

Summer  2  to  4 
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3.2  Methods 


In  the  fully-rough  regime,  the  wind  velocity  profile  under  neutral 
conditions  of  stability  is  related  to  the  three  wind  profile  parameters ' 
by  the  equation  (Lettau,  1957) 

V  =  V*  ln[(z  +  d  +  Zo)/zo]  (cm/sec), 

where  V  =  wind  velocity  (cm/sec)  at  height  z ,  and  V*  =  tq/p 

=  2.5 n/  To/p  =  friction  velocity  (cm/sec),  with  p  =  air  density  (g/cm^) 
and  k  s  0.4  =  universal  Karman  number.  For  convenience,  let 
D  =  d  +  Zo  (cm).  ; 

A  curve-fitting  is  most  satisfactorily  achieved  by  a  least  square 
method  (Lettau.  1957)  in  which  the  parameters  are  determined  such 
that  the  sum  of  the  error  squares  is  minimized, 

tj  =  Vj-V*  ln[(Zj  +  D)/zo]  ; 

=  0  ;  and  Z'!  '  minimum  . 

For  practical  computations,  a  numerical  trial -and -error  method 
may  be  used.  One  assumes  an  arbitrary  D-value  which  is  equal  to 
or  smaller  than  zero,  and  repeats  the  computation  for  other  D-value s, 
until  that  D  is  determined  for  which  the  sum  of  the  error  squares  has 
the  least  value.  When  the  bar  indicates  the  mean  for  the  i  observa¬ 
tion  levels  (z^) ,  the  following  abbreviations  are  introduced: 


X.  =  V.  -  V  ;  =  In  (z^  +  D)  -  In  (z  +  D)  ; 


Then,  the  equations, 

£,  =  X  -  V*Y.  ;  and  V*  =  S  /S 
1  i  1  xy  XX 

give  the  deviations  and  their  sum  of  squares  varies  directly  as  a 
function  of  D. 

After  determining  the  least  square  D-value,  the  roughness  param¬ 
eter.  Zo  ,  is  computed  by  the  equation 


In  Zo  =  ln(z  +  D)  -  V/V*  . 


Table  2.  Evaluation  of  wind  profile  data  for  relatively  extreme  surface  roughness  types. 
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Then  the  relation  D  s  d  *1*  Zq  gives  the  zeropoint-displacement,  d , 
and  the  equation 

V*  s  2.5 '/TJTT 


gives  the  surface  stress  parameter,  tq  .  Air  density  p  is  assumed 
to  be  0.0012  g/cm*. 

3.3  Data 

A  search  of  the  literature  revealed  that  wind  profile  measurements 
over  tall  vegetation  are  very  rare.  In  this  computation,  data  from  five 
different  sources  — The  Johns  Hopkins  University  Laboratory  of  Clima¬ 
tology  (1952),  Fons  (1940),  Kepner,  Boelter  and  Brooks  (see  Poppen- 
diek,  1949),  Baumgartner  (1956,  1957),  and  Gisborne  (see  Poppendiek, 
1949)  —  are  used  to  obtain  roughness  parameter  values  for  corn  field, 
brush  field,  orange  orchard,  fir  forest,  and  Idaho  forest.  The 
zeropoint-displacement  and  surface  stress  are  also  computed  in  refer¬ 
ence  to  roughness  parameter. 

As  the  wind  profile  will  be  logarithmic  only  under  neutral  stability 
conditions,  attention  was  paid  to  selecting  wind  profile  data  for  adia¬ 
batic  hours,  if  it  was  possible. 

3.4  Roughness  Parameter  Approximation 


The  Johns  Hopkins  University  (1952)  wind  velocity  data  for  seven 
levels  from  1.6  to  6.4  m  height  at  Shirley  Field  Station,  New  Jersey, 
was  used  to  compute  roughness  parameter  for  a  thickly  planted  corn 
field  before  tasseling  when  the  plant  height  was  about  0.7  to  0.9  m. 
The  data  in  Table  2  refer  to  the  wind  velocity  profile  at  1800  to  1900 
local  time,  July  2,  when  the  lower  atmosphere  was  adiabatic.  The 
computed  roughness  parameter,  Zq  ,  was  6.4  cm  with  d  =  -88.4  cm 
and  To  =  2.2  dynes/cm^  .  This  Zq  value  agrees  closely  with 
Deacon' s  (1953)  value  for  long  grass. 

The  wind  velocity  measurement  by  Fons  (1940)  above  a  brush 
field  for  seven  levels  from  1.5  to  6.4  m  in  California  yielded  a  rough¬ 
ness  parameter  value  of  15.7  cm  with  d  =  -125  cm  and  tq  =  3.8 
dynes/cm^  .  The  brush  was  manzanita  and  snowbrush  of  1.3  m  average 
height. 

Wind  velocities  for  seven  levels  from  4.3  to  15.2  m  were  measured 
above  an  orange  orchard  by  Kepner.  Boelter  and  Brooks  (Poppendiek, 
1949).  Data  were  reported  in  five  groups.  Group  IV  (V  =  1.  36  m/sec) 
and  group  V  ( V  =  2.29  m/sec)  of  relatively  high  mean  velocity  were 
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selected  for  the  analysis.  On  the  average  the  orange  trees  were 

3.3  m  tall  and  had  a  diameter  of  3.3  m.  The  wind  observations  were 
made  at  night  under  temperature  inversion  conditions.  The  roughness 
parameter  value  obtained  from  group  IV  was  24. 1  cm  with  d  =  -332  cm 
and  To  =  0. 4  dynes/cm^  which  seemed  comparable  with  the  values 
obtained  by  Deacon  (Table  1)  for  other  natural  coverages.  However, 
the  computed  roughness  parameter  value  from  group  V  was  122  cm 
with  d  s  -243  cm  and  to  -  3.2  dynes/cm^ ,  which  seemed  unreas¬ 
onably  large.  This  might  be  a  spurious  effect  due  to  the  temperature 
inversion  condition  under  which  these  measurements  were  made. 

Two  sets  of  wind  velocity  data  above  fir  forest  were  selected 
from  Baumgartner' s  work  (1956,  1957).  One  set  represents  the  aver¬ 
age  wind  velocity  profile  for  1700  to  1800  local  time  from  June  28  to 
July  8  of  1952,  and  another  set  represents  the  average  for  1800  to 
2100  local  time  from  July  9  to  July  14  of  1952.  In  both  cases,  the 
data  refer  to  isothermal  condition.  The  tree  height  was  about  5.5  m. 
The  four  anemometer  levels  were  from  8  to  16  m.  The  computed  Zq 
value  for  the  first  set  was  190  cm  with  d  =  -190  cm  and  t©  =4.9 
dynes/cm* ;  for  the  latter  Zq  =  399  cm  with  d  =  -399  cm  and 
To  =  2.4  dynes/cm^  (Table  2).  Both  Zq  values  are  unexpectedly  large 
even  though  the  computation  was  carried  only  to  D  =  0  in  order  to 
avoid  positive  D-values. 

Two  groups  of  Idaho  forest  wind  velocity  data  were  presented  by 
Gisborne  (see  Poppendiek,  1949)  for  three  anemometer  levels  from 

25.3  m  to  47.5  m.  The  average  tree  height  was  about  24.4  m  for 
spruce,  hemlock,  and  Douglas  fir.  The  computed  roughness  parameter 
for  group  I  was  1847  cm,  with  d  =  -1847  cm  and  tq  =83.1  dynes/cm^ 
for  group  II  the  computed  roughness  parameter  value  was  1347  cm  with 
d  =  -1637  cm  and  tq  =  6.  8  dynes/cm^  (Table  2).  In  both  cases,  the 
variation  of  D  was  restricted  to  positive  values  only.  The  unreason¬ 
ably  large  roughness  parameters  suggested  unreliability  of  the  data. 
The  large  tq  value  for  group  I  does  not  seem  to  be  realistic. 

3.5  Reliability  of  the  Data 

The  results  of  roughness  parameter  computation  will  strongly  de¬ 
pend  on  the  reliability  of  the  anemometer  data. 

The  Johns  Hopkins  University  (1952)  data  for  corn  field  include 
seven  anemometer  levels.  Special  care  was  taken  to  use  only  the 
most  reliable  and  matched  anemometers.  The  error  distribution  in  the 
least-square  computation  was  found  to  be  very  even.  To  check  the 
reliability  of  the  roughness  parameter  computation,  the  data  were 
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divided  into  two  parts:  the  first  part  contained  the  four  Z|  with 
i  s  1,  3,  5,  and  7 :  and  the  second  part,  the  three  with  i  »  2. 4.  6. 
The  roughness  parameter  was  computed  for  each  group  separately, 


and  the  results  are  listed  below. 

Levels  (m) 

Zo  (cm) 

d  (cm) 

To  (dynes/cm^) 

1.6,  2.3.  4.8,  6.4 

3.  0 

-106 

1.6 

2.  4,  4.  0.  5.  6 

14.7 

-  57.7 

3.2 

All  .7  levels 

6.4 

-  88.4 

2.2 

The  difference  in  computed  roughness  parameter  values  from  the 
complete  data  was  considerable.  The  probable  error  of  the  result 
from  the  complete  data  may  be  expressed  in  terms  of  a  factor  which 
can  vary  from  j  to  2 . 


The  same  treatment  was  applied  to  Fons'  (1940)  data  for  brush 
field.  The  original  set  contained  seven  levels,  and  the  two  divided 
parts  had  four  and  three  levels  each.  The  results  are  listed  below. 


Levels  (m) 

Zo  (cm) 

d  (cm) 

To  (dynes/cm^ ) 

1.  5,  2.  4,  4.  3,  6.  4 

7.  3 

-135 

2.  5 

1.8,  3.7,  4.6 

87.  6 

-  87.  6  (D=0)  13.  7 

All  7  levels 

15.  7 

-125 

3.  8 

The  differences  in  computed  roughness  parameter  value  are  such 
that  the  probable  error  may  be  expressed  by  a  factor  which  may  vary 
from  i  to  4. 

Both  groups  of  Baumgartner' s  (1956,  1957)  data  for  fir  forest  con¬ 
tained  only  four  anemometer  levels.  On  the  basis  of  the  error  distri¬ 
bution,  the  third  level  seemed  inaccurate.  Therefore  this  anemome¬ 
ter  level  was  omitted  and  Zq  was  computed  for  the  remaining  three 
levels.  The  results  are  listed  below. 

The  group  II  data  showed  some  improvement  by  eliminating  the 
third  level,  though  the  Zq  value  still  remained  larger  than  expected. 
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However,  group  I  data  was  not  affected  by  elimination  of  the 
third  level. 


Levels  (m) 

Zo  (cm) 

d  (cm)  To 

(dynes/cm^ ) 

GrouD  I 

8.0,  9.8, 

12.2, 

16.0 

190 

-190  (D=0) 

4.9 

8.0,  9.8, 

16.0 

192 

-192  (D=0) 

5.0 

Group  II 

8.0.  9.8, 

12.2, 

16.0 

399 

-399  (D=  0) 

2.4 

8.0,  9.8, 

16.0 

203 

-493 

1.4 

Since  the  Idaho  forest  data  by  Gisborne  had  only  three  levels  it 
was  Impossible  to  omit  any  information  from  the  data. 

This  study  leads  to  the  conclusion  that  caution  should  be  used 
in  future  experiments  of  this  nature.  The  considerable  unreliability 
of  previous  data  may  include  the  construction,  exposure,  condition 
of  the  anemometer.  It  is  suggested  that  at  least  seven  levels  are 
required  for  a  reliable  roughness  parameter  approximation.  Wind 
data  at  less  than  four  levels  must  be  considered  inconclusive. 

3.6  Concluding  Remarks 

Plotting  the  roughness  parameter  values  versus  height  of  natural 
vegetation  on  a  logarithmic  graph,  a  straight  regression  line  can  be 
assumed.  It  may  be  stated  that  the  logarithm  of  roughness  parameter 
is  proportional  to  the  logarithm  of  height  of  natural  vegetation.  The 
regression  line  is  expressed  by 

log  Zo  =  -1.24  +  1.19  log  h 

where  Zq  and  plant  height  h  are  in  cm;  see  Fig.  1,  Section  5.4. 

The  surface  stress,  tq  ,  can  be  computed  independently  to  assign 
the  roughness  parameter  value  or  check  the  reliability  of  the  data  for 
the  specific  plant  height. 


Roughness  Parameters  Above  Tali  Vegetation 


35 


References 


Baumgartner,  A. :  Untersuchungen  iiber  den  W&rme  —  und  Wasser- 
haushalt  eines  Jungen  Waldes.  Berlchte  des  Deutchen  Wetter* 
dienstes  No.  28  (Band  5),  1956. 

Baumgartner,  A. :  Beobachtungswerte  und  weitere  Studien  zum 

Warme  —  und  Wasserhaushalt  eines  Jungen  Waldes.  University 
MUnchen  —  Meteorologisches  Institut,  Wissenschaftliche 
Mitteilunaen  No.  4,  1957. 

Deacon,  E.  L. :  Vertical  profiles  of  mean  wind  in  the  surface  layer  of 
the  atmosphere.  Geophysical  Memoirs  2 .  No.  91,  1953. 

Fons,  W.  L. :  Influence  of  forest  cover  on  wind  velocity.  Journal  of 
Forestry  38 .  No.  6,  1940. 

The  John  Hopkins  University  Laboratory  of  Climatology:  Publications 
in  ClimatoloQV  V ,  No.  7,  1952. 

Lettau,  H.  H. :  Computation  of  Richardson  Numbers,  Classification 
of  Wind  Profiles,  and  Determination  of  Roughness  Parameters. 
Section  7.4  in  Vol.  I  of  Exploring  the  Atmosphere*  s  First  Mile 
(H.  Lettau  and  B.  Davidson,  editors),  Pergamon  Press,  London 
and  New  York,  1957. 

• 

Poppendiek,  H.  F. :  Investigation  of  velocity  and  temperature  profiles 
in  air  layers  within  and  above  trees  and  brush.  Contract  N6- 
ONT-275,  Task  Order  VI,  NR-082-036.  Department  of  Engineer¬ 
ing,  University  of  California,  Los  Angeles,  1959. 


4 


Aerodynamic  Drag  In  Tall  Vegetation 

Robert  H.  Burgy* 

Department  of  Civil  Engineering 
and 

Department  of  Meteorology 
University  of  Wisconsin 


Abstract.  This  report  summarizes  the  results  of  a  review 
of  recent  literature  relating  to  drag  forces  and  coefficients 
for  tall  vegetation.  Results  of  wind  tunnel  tests  on  conifer 
saplings,  and  field  tests  on  broad  leafed  trees  are  quoted 
which  show  that  the  drag  coefficient  deaeases  with  inaeas- 
ing  wind  speed. 

4.1  Introduction 

The  study  by  Kung  —  see  Section  3  —  reviewed  the  available  infor¬ 
mation  on  wind  profiles  above  stands  of  shrubs  and  trees.  This  review 
is  directed  mainly  toward  locating  additional  data  on  which  to  base 
calculations  of  roughness  effects  for  trees. 

No  new  wind  profile  data  were  located.  However,  a  group  of 
studies  relating  to  the  effects  of  winds  in  tree  damage  in  forests  from 
atomic  blasts  was  disclosed.  Two  studies  have  been  made,  one  on 
conifers  —  Sauer,  Fons,  and  Arnold  (1951)  — and  the  other  on  hard¬ 
woods  —  Lai  (1955) —  wherein  whole  individual  trees  were  used  to 
measure  drag  in  the  crowns.  Part  of  the  worlc  was  done  in  a  wind 
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tunnel  and  the  bulk  of  the  tests  were  made  on  truck-mounted  tree 
samples  moved  at  constant  velocity  while  drag  and  other  measure¬ 
ments  were  made.  These  tests  led  to  studies  of  crown  characteristics 
of  forest  types  to  relate  aown  areas  to  crown  weights,  with  two  addi¬ 
tional  reports  being  published  on  tree  crown  relationships  *  Storey, 
Fons,  and  Sauer  (1955),  and  Storey  and  Pong  (1957). 

In  hydrodynamic  work  it  is  convenient  to  express  the  frictional 
resistance  of  bodies  immersed  in  flowing  fluids, as  a  "drag  force,"  Df , 
in  appropriate  units.  This  drag  force  is  related  to  the  dynamic  pres¬ 
sure  acting  on  the  body  by  a  suitable  factor,  the  drag  "coefficient," 
which  thus  gives  the  drag  as  a  function  of  dynamic  pressure  on  the 
cross-sectional  area  of  the  body  exposed  to  the  flow,  as  follows: 

Dj  =  CjApV*  /2  ;  or  Cj  =  2Dj/ApV* 

where  =  drag  coefficient,  V  =  velocity,  A  -  silhouette  area  and 
p  =  density  of  fluid. 

The  data  are  for  single  trees  of  various  species  moved  through  the 
air  on  a  vehicle.  The  basic  data  for  the  conifer  tests  were  not  reported 
and  it  is  not  possible  to  make  further  analysis  of  these  data,  since  the 
study  was  concerned  with  shock  waves  from  blasts  and  the  relation¬ 
ships  are  in  nondimensional  form. 

Tree  crown  characteristics  were  developed  to  define  the  several 
types  with  respect  to  crown  weights  and  area.  In  the  limited  scope 
of  this  review  the  writer  has  not  been  able  to  make  more  than  very  . 
generalized  comparisons  to  determine  the  resultant  drag  coefficients 
for  the  two  tree  types.  Some  estimates  based  on  average  values  with 
certain  assumptions  concerning  the  trees  are  made  to  obtain  an  average 
drag  coefficient. 

Another  study  reported  from  Japan  by  Hirata  (1953)  includes  wind 
tunnel  investigations  of  tree  crown  drag  for  "models"  of  forest  trees. 
Drag  coefficients  are  reported  which  are  of  the  order  of  1.0.  The 
study  was  also  directed  toward  wind  damage  effects-  and  is  conse¬ 
quently  of  limited  applicability  to  the  problem  of  forest  drag  under 
normal  conditions. 

the  thought -provoking  question  may  be  asked  concerning  the  re¬ 
lationship  of  values  measured  in  single  tests  as  compared  to  the  com¬ 
posite  conditions  in  a  stand  of  trees.  Certainly  the  drag  forces,  and 
coefficients  are  applicable  for  the  given  velocity.  However,  the  de¬ 
velopment  of  the  boundary  layer  and  the  consequent  redistribution  of 
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the  wind  profile  over  the  forest  will  cause  marked  differences  in 
comparison  with  drag  for  a  single  tree.  Additionally  the  integrated 
effects  of  uniform  or  variegated  vegetation  and  existence  of  both 
primary  and  secondary  canopies  will  be  important  in  micrometeorol- 
oglcal  work;  reference  is  made  to  Kutzbach' s  discussion,  see 
Section  7.5.4. 

The  application  of  mass  transfer  concepts  to  the  forest  for  micro- 
meteorological  investigations  requires  a  considerably  better  under-^ 
standing  of  the  internal  relationships  for  forests  and  tall  vegetation 
than  is  now  available. 

An  additional  topic  reviewed  in  this  study  was  the  effect  of  wind 
speed  on  drag  coefficients.  for  rigid  bodies  of  many  shapes  have 
been  reported.  Flexible  bodies  like  vegetation  have  rarely  been 
studied.  Some  work  has  been  done  on  fabrics,  flags,  tow  targets, 
and  parachutes,  which  may  have  some  bearing  on  the  question;  see 
Hoemer  (1951,  1958). 

Sauer  (1951)  and  Lai  (1955)  recognized  the  effect  of  wind  on  Cd 
and  suggest  that  plant  components  may  respond  to  increased  dynamic 
pressure  by  deflecting, .  with  a  resultant  variation  in  drag  coefficient. 
If  the  plant  tends  to  become  more  streamlined  in  shape  under  higher 
pressures  the  drag  parameter  would  become  commensurateiy  smaller, 
making  a  non-linear  relation  between  squared  velocity  and  stress. 

If,  however,  in  some  plants  the  higher  speeds  will  cause  flutter  or 
flapping  and  oscillations  of  the  leaves  and  branches,  then  the  drag 
may  become  greater  and  the  coefficient  increased. 

Hirata'  s  (1953)  data  show  an  effect  of  the  wind  speed,  increas¬ 
ing  in  some  species  and  decreasing  in  others.  He  concludes  that 
there  is  no  effect  for  the  conditions  of  his  tests  and  attributes  the 
variation  to  errors.  The  trends  are  quite  apparent  and  cause  this 
writer  to  suspect  that  a  relationship  does  exist. 

Observational  evidence  seems  to  support  either  contention  and 
detailed  investigation  must  be  made  to  resolve  the  question.  The 
thought  also  occurs  that  the  oscillations  in  grain  fields,  or  waves 
moving  across  the  stalks  (seen  as  bending  of  the  stalks)  may  be 
partly  due  to  fluctuations  of  Cc[  as  the  stalks  bend,  introducing  a 
harmonic  cycle  to  the  wind  stream,  which  would  rise  and  fall  both 
in  vertical  position  and  in  magnitude  as  the  boundary  layer  responded 
to  a  changing  surface  stress.  Conversely,  it  may  well  be  concluded 
on  the  basis  of  even  fragmentary  evidence  that  any  variation  of  drag 
with  wind  speed  is,  for  practical  purposes,  of  no  significance  in 
micrometeorological  work. 
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4.2  Wind  Tunnel  Tests  on  Conifer  Saplings 

The  data  presented  in  the  reports  relating  to  crown  drag  in  coni¬ 
fers  and  broadleafed  (hardwood)  species  have  been  used  to  calculate 
some  typical  values  of  Cd  for  the  conditions  of  the  tests.  Some  as¬ 
sumptions  are  made  of  the  effective  area  of  crown  exposed  since  the 
reports  use  the  relationship  of  foliage  and  branch  weights  to  correlate 
drag  force  with  a  moment-stress -dynamic  pressure  factor.  There  is 
some  evidence  to  indicate  that  the  use  of  foliage  dry-weight  is  appro¬ 
priate  and  some  correlations  are  available.  For  these  calculations 
the  writer  preferred  to  use  standard  terms.  The  basic  data  are  not 
complete  for  the  conifer  tests  of  whole  trees.  Even  the  wind  tunnel 
tests  are  not  fully  explained,  and  possible  misinterpretation  of  the 
graphs  may  be  inherent  in  the  calculated  values  presented  here. 


For  this  purpose  it  was  believed  appropriate  to  estimate  the  gross 
silhouette  area  of  the  tree  species  based  on  some  average  dimensions 
and  shapes.  This  is  a  reasonable  assumption  since  the  coefficients 
should  be  calculated  on  the  standard  terms  used  in  hydrodynamics 
and  aerodynamics  to  be  meaningful. 


As  noted  above,  the  tests  in  the  wind  tunnel  of  saplings  of  sev¬ 
eral  conifer  species  are  graphed  using  nondimensional  parameters  and 
therefore  with  incomplete  explanation  this  reviewer  is  unable  to  in¬ 
sure  the  meaning  of  the  graphs.  If  a  correct  interpretation  is  made, 
then  the  drag  coefficients  of  white  fir  saplings  for  several  wind  speeds 
are  as  listed  below.  The  average  height  was  0.6  m,  the  average 
width  was  0.3  m.  The  average  dry  foliage  weight  was  180  grams.  A 
gross  silhouette  area  of  0.18  m^  and  an  air  density  of  0.00106  g/cm^ 
was  assumed. 


Velocity 

23 

15 

9 

(m/sec) 

Drag  force 

2.  53 

1.  60 

0.  70 

(10^'dynes) 

Drag  coefficient 

0.  41 

0.  60 

0.  72 

(  -  ) 

It  should  be  noted  that  these  are  '‘apparent"  drag  coefficients 
since  the  area  is  gross  silhouette.  As  the  velocity  increases  the  tree 
tends  to  deflect  downwind  and  this  action  seems  to  result  in  a  more 
compact  silhouette.  It  is  impossible  to  estimate  the  amount  of  this 
effect  on  wind  speed  vs.  drag  since  the  reduced  area  of  the  tree  must 
yield  less  drag  or  a  higher  coefficient  and  using  the  constant  area 
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here  shows  a  reducing  coefficient.  These  two  effects  must  be  par¬ 
tially  compensating. 

The  tests  reported  for  the  conifer  trees  on  vehicles  do  not  include 
sufficient  detail  or  data  to  permit  a  calculation  of  .  Several  at¬ 
tempts  were  made  to  select  average  values  from  the  graphs  and  to  re¬ 
late  these  to  the  wind  tunnel  tests  which  were  also  plotted  with  the 
field  data.  These  attempts  have  proved  fruitless  to  this  writing  since 
the  values  of  thus  obtained  range  from  over  50  to  about  3.  These 
are  believed  to  be  meaningless  due  either  to  inaccuracy  or  probable 
erroneous  interpretation  of  the  data  as  presented.  Inquiries  have  been 
made  to  the  authors  for  more  explanation  of  the  methods  used  and 
meaning  of  the  terms. 

4.3  Field  Tests  on  Broadleafed  Trees 


The  tests  of  the  broadleaf  species  provide  more  dependable  re¬ 
sults  since  the  basic  data  are  reported  for  the  runs  on  the  various 
t^'oes.  The  values  obtained  for  three  speeds  of  tree  moved  on  the 
truck  are  listed  below.  Assumptions  were  made  of  the  tree  crown 
area  and  shape  as  before.  The  total  tree  height  was  12  to  20  m ,  the 
average  crown  height  7.5  m.  The  gross  silhouette  area  was  assumed 
to  equal  13.5  m^. 


Velocity 

23 

10 

6 

(m/sec) 

Drag  force 

133 

51 

24 

(10^  dynes) 

Drag  coefficient 

0.  31 

0.76 

0.  82 

(  -  ) 

It  should  be  noted  that  the  use  of  gross  tree  area  is  not  represen¬ 
tative. of  the  actual  situation  since  the  broadleafed  species  are  rela¬ 
tively  open  in  structure.  The  Cd  values  thus  calculated  show  the 
reduction  with  increase  in  velocity  as  did  the  conifer  saplings. 

There  is  no  rational  way  in  which  this  deflection  of  the  tree  form 
and  streamlining  can  be  identified  from  the  data  presented  in  these 
studies.  The  writer  is  confident  that  a  means  can  be  found  to  measure 
this  effect  and  to  isolate  the  relation  of  wind  speed  and  drag. 

In  both  conifer  and  broadleaf  tests  some  specimens  were  defoli¬ 
ated  and  tested  with  the  drag  being  reduced  from  l/lO  to  l/5  of  the 
whole  crown  drag  values. 


4.4  Addenda 


Data  furnished  by  Lemon  (1960),  in  response  to  the  author' s  re¬ 
quest,  show  relationships  for  three  species  of  aops  —  corn,  alfalfa, 
and  wheat  —  of  the  wind  speed  and  the  roughness  parameter  Zq  .  , 

For  winds  ranging  from  200  to  1000  cm/sec  measured  above  the 
crop,  the  Zq  -values  vary  as  follows: 

Alfalfa  (75  cm)  Zq  decreases  continuously  — 

Wheat  (130  cm)  Zq  rises,  then  falls  off  rapidly  — 

Com  (240  cm)  Zq  increases  continuously  — 

—  all  as  the  wind  speed  increases.  The  displacement  height  d  varies 
inversely  as  the  Zq  value. 

The  curves  are  suggestive  that  for  winds  greater  than  those  shown 
for  the  corn,  this  crop  would  probably  show  a  reversal.  This  is  pos¬ 
sible  since  the  dynamic  pressure  of  greater  wind  will  tend  to  bend  the 
drop  over  to  a  more  streamlined  shape.  Some  opinions  have  been  ex¬ 
pressed  that  crops  such  as  corn  with  long  flexible  leaves  that  flap  • 

probably  would  h^ve  greater  drag  coefficients  under  higher  wind  speeds, 
up  to  the  point  that  physical  damage  to  the  plant  might  occur.  Lemon' s 
data  do  not  cover  this  range  but  do  show  the  increased  Zq  which  evi¬ 
dences  the  greater  drag  at  higher  velocities. 

It  is  interesting  to  note  that  this  information  supports  some  of  the 
contentions  outlined  previously  in  this  report.  On  the  basis  of  these 
relationships  it  would  be  appropriate  to  conclude  that  for  long  stemmed 
and  flexible  bladed  crops  and  vegetation  Zq  does  increase  with  wind. 

That  alfalfa  is  reversed  follows  logically  since  it  is  a  small  leafed 
plant.  A  similar  relationship  undoubtedly  exists  in  the  several  arrange¬ 
ments  of  leaves,  needles,  branches,  etc. ,  in  trees.  The  botanical 
structure  of  tree  leaves  and  needles,  etc. ,  should  be  investigated  in 
the  light  of  these  concepts  regarding  velocity-drag  relationships. 

Sauer  (1951)  and  Lai  (1955)  tend  to  lead  toward  this  in  analyzing  the 
individual  species  which  they  studied. 
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Abstract*  Data  on  the  areal  distribution  of  defined  types 
of  plant  cover  are  evaluated  for  the  region  of  the  State  of 
Wisconsin,  and  for  latitudinal  zones  of  all  continents. 

An  empirical  numerical  relationship  between  plant  height 
and  aerodynamic  roughness  parameter  (zq  )  is  used  to 
estimate  weighted  area-averages  of  log  Zq  .  In  an  at¬ 
tempt  to  establish  a  basis  for  utilizing  micrometeorolog- 
ical  information  in  large-scale  meteorological  problems, 
the  regional  and  meridional  distributions,  and  their  sea¬ 
sonal  variations,  of  vegetative  surface  roughness  is 
discussed. 

5. 1  Introduction 


Frictional  forces  in  the  lower  atmosphere,  and  the  dissipation  of 
mechanical  energy  of  large-to-small-scale  air  currents  depend  prima¬ 
rily  on  the  roughness  structure  of  the  natural  ground.  Micrometeorol- 
ogical  information  indicates  that  the  frictional  effects  at  the  lower 
boundary  can  be  expressed  in  terms  of  the  aerodynamic  roughness 
parameter  (z©,  cm),  which  depends  significantly  on  the  type  of  vege¬ 
tation  cover.  In  a  preceding  report  (9)  a  numerical  relationship  be¬ 
tween  the  roughness  parameter  and  the  plant  height  was  extended  to 
include  several  types  of  tall  vegetation.  In  order  to  obtain  the  rough¬ 
ness  parameter  for  a  region,  in  contrast  to  that  of  a  defined  micro- 
meteorological  site,  it  is  necessary  to  combine  the  above  relationship 
with  data  on  the  horizontal  distribution  of  plant  covers. 

In  the  present  paper,  an  attempt  is  made  to  study  in  detail  the 
probable  seasonal  variation  of  aerodynamic  roughness  for  a  region  of 
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approximately  10*  km*  as  represented  by  the  State  of  Wisconsin. 

Also,  available  ecological  information  for  all  continents  is  used  to 
prepare  the  basis  of  a  world-wide  study  of  the  meridional  distribution 
of  roughness  parameters  for  all  land  surfaces. 

5.2  Method  and  Data 

The  evaluation  of  the  aerodynamic  roughness  of  a  selected  region 
involves  the  following  four  steps: 

(1)  Classify  the  vegetation  into  major  types  and  determine  the  area- 
percentage  of  each  type  (i.  e. ,  the  "land  use"). 

(2)  Estimate  representative  heights  for  each  vegetation  type  with  a 
breakdown  by  month,  or  by  season. 

(3)  On  the  basis  of  (2)  assign  roughness  parameters  to  each  vegeta¬ 
tion  type. 

(4)  Compute  the  regional  mean  of  the  logarithm  of  roughness  param¬ 
eters,  employing  the  percentage  area  of  each  vegetation  type  as 
a  weight  factor,  with  a  breakdown  by  month,  or  by  season. 

For  the  regional  study  of  the  State  of  Wisconsin,  the  data  on  land 
use  were  obtained  from  the  following  sources:  statistics  by  the  Wis¬ 
consin  Crop  and  Livestock  Reporting  Service  (20)  and  the  Conservation 
Department  (10,  19),  and  estimates  by  the  staff  of  the  Forestry  Exten¬ 
sion  Service  of  the  University  of  Wisconsin  (8).  For  the  continent- 
wide  study  the  statistics  of  the  United  Nations  (6)  and  other  sources 
(2,  5,  7,  11,  13,  14.  15,  16,  17)  were  used. 

Phenological  data  and  estimates  of  plant  heights  were  obtained 
from  various  sources  (1,  3,  12,  18),  and  also  as  the  result  of  discus¬ 
sions  with  staff  members  of  the  Departments  of  Agronomy  and  Botany 
at  the  University  of  Wisconsin.  For  the  estimate  of  roughness  param¬ 
eters,  on  the  basis  of  plant  height,  the  numerical  relationship  derived 
by  E.  Kung  (9)  was  used. 

5.3  Monthly  Data  on  Surface  Roughness  for  the  Region  of  Wisconsin 

The  area  of  the  State  of  Wisconsin  was  divided  into  nine  districts 
on  a  county  basis,  in  the  following  manner. 
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District 


Counties  Included 


1.  Northwest  Baxron,  Bayfield,  Burnett,  Chippewa, 

Douglas,  Polk,  Rusk,  Sawyer,  Washburn 


2.  North  Ashland,  Clark,  Iron,  Lincoln,  Marathon, 

Oneida,  Price,  Taylor,  Vilas 

3.  Northeast  Florence,  Forest,  Langlade,  Marinette, 

Oconto,  Shawano 


4.  West 


5.  Central 


6.  East 


7.  Southwest 


Buffalo,  Dunn,  Eau  Claire,  Jackson, 

La  Crosse,  Monroe,  Pepin,  Pierce, 

St.  Croix,  Trempealeau 

Adams,  Green  Lake,  Juneau,  Marquette, 
Portage,  Waupaca,  Waushara.  Wood 

Brown,  Calumet,  Door,  Fond  du  Lac, 
Kewaunee,  Manitowoc,  Outagamie, 
Sheboygan,  Winnebago 

Crawford,  Grant,  Iowa,  Lafayette. 
Richland,  Sauk,  Vernon 


8.  South 


Columbia,  Dane,  Dodge,  Green,  Jeffer 
son.  Rock 


9.  Southeast  Kenosha,  Milwaukee,  Ozaukee,  Racine, 

Walworth,  Washington,  Waukesha. 

According  to  land  use  in  these  districts,  nine  major  categories  of 
vegetation  cover  were  distinguished.  Each  category  involves  either 
one,  or  more  than  two  main  kinds  of  vegetation,  as  follows: 

(1)  Corn  -  (2)  Oats  -  (3)  Hay  -  (4)  Other  field  crops  -  (5)  Natural 
vegetations  -  (6)  Oak  -  (7)  Maple  -  (8)  Aspen  -  (9)  Other  trees.  The 
results  of  the  vegetation  survey  for  each  district  are  shown  in  Table  1, 
in  percentages  of  total  district  area. 

Combining  the  vegetation  cover  survey,  phenological  data  and 
estimates  of  vegetation  height,  the  average  monthly  roughness  param¬ 
eter  was  evaluated  for  each  district  of  Wisconsin.  Results  are  shown 


Table  1.  Percentage  area  covered  with  indicated  vegetation  types  for  nine  districts  of  the  State  of  Wisconsin. 
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in  Table  2.  For  this  evaluation  a  roughness  parameter  for  snow  cover 
of  Zq  -  0. 10  cm  was  adopted  from  Deacon' s  study  (4). 

The  distribution  of  average  roughness  parameters  over  different 
districts  and  their  seasonal  change  can  be  summarized  as  follows. 

The  northern  parts  of  the  state  have  high  annual  Zq  -averages  — 
zo  s  64.96  cm  in  Northeast.  44.62  cm  in  North,  and  38.85  cm  in 
Northwest.  The  zq  -values  decrease  towards  the  Southwest  and 
South.  The  annual  average  Zq  -values  are  Zo  =  8. 12  cm  in  West, 

12.71  cm  in  Central,  5.24  in  East,  3.19  in  Southwest,  4.29  in  South 
and  6.82  cm  in  Southeast.  The  regional  difference  of  characteristic 
roughness  parameter  is  mainly  due  to  heavy  acreage  of  forest  trees 
in  the  northern  parts  of  the  state,  while  in  the  southern  parts  of  the 
state  farm  crops  and  grass -type  natural  vegetation  dominate.  The 
West  and  Central  districts  are  a  zone  of  transition.  Since  the  Zq  - 
value  is  proportional  to  the  1.19th  power  of  the  vegetation  height, 
these  characteristic  differences  in  vegetation  cover  are  emphasized. 
The  seasonal  change  of  roughness  parameters  is  also  significant  for 
the  whole  state  as  well  as  each  district.  This  is  due  to  the  phenolog- 
ical  cycle.  The  lush  summer  vegetation  produces  the  highest  Zo  - 
values.  The  spring  and  fall  values  are  intermediate  between  the 
summer  maximum  and  winter  minimum.  We  also  notice  that  the  sea¬ 
sonal  change  of  roughness  parameter  is  greater  in  the  southern  and 
central  part  of  the  state  than  in  the  northern  parts.  This  is  illustrated 
by  the  following  ratio  of  winter  to  July  Zo  values: 


Districts 

Winter  Zo  / July  Zo  (% ) 

Northwest 

20.7 

North 

24.7 

Northeast 

27.6 

West 

5.6 

Central 

8.0 

East 

3.  7 

Southwest 

3.  4 

South 

2.  9 

Southeast 

4.  0 

State  Average 

9.6 

This  is  caused  by  the  withering  of  field  crops  and  grass -type  vegeta¬ 
tions,  with  snow  cover  over  these  low  vegetation  areas  in  southern 


Table  2.  Seasonal  trend  of  estimated  area-averages  of  the  aerodynamic  roughness  parameter  (Zq.  cm)  for  nine 
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and  central  parts  of  the  state.  The  effect  of  snow  cover  is  much  re¬ 
duced  in  the  forested  northern  parts  of  the  state. 

5.4  Meridional  Distribution  of  Vegetlve  Roughness 

The  land  area  of  the  world  was  studied  separately  for:  (1)  the 
Americas,  (2)  Europe-Africa,  and  (3)  Asla-Oceania  (including  the 
USSR),  with  a  subdivision  into  10-degree  latitudinal  zones,  and 
5-degree  latitudinal  zones  between  30  and  50*N  in  North  America. 

The  estimates  of  the  aerodynamic  roughness  parameter  of  all  zones, 
for  the  seasons,  were  based  on  surveys  of  the  land  use  and  the  dis¬ 
tribution  of  vegetation  cover. 

For  the  United  States  and  Canada  the  main  sources  of  information 
were  the  USDA  Official  Statistics  (15,  16)  and  the  Canada  Yearbook 
(5).  The  data  for  China,  India,  and  Australia  were  taken  from  Shen 
(13),  Time  of  India  (14),  and  the  Yearbook  of  the  Commonwealth  of 
Australia  (2),  respectively,  together  with  a  survey  of  existing  maps 
(7,  17).  The  1957  FAO  Yearbook  (6)  was  used  for  all  other  areas. 
Atlasses  (7,  17)  and  other  miscellaneous  Information  (1,  8,  16,  etc.) 
were  also  considered  for  bridging  gaps,  and  to  make  slight  adjust¬ 
ments  in  a  few  cases  where  the  available  statistics  data  were  insuf¬ 
ficient  or  incomplete. 

The  following  eight  classes  of  continental  vegetation  cover  were 
distinguished:  (1)  Crops  -  (2)  Woodland  and  tree  crops  -  (3)  Meadow 
and  pasture  -  (4)  Forest  -  (5)  Natural  vegetations  -  (6)  Tundra  - 
(7)  Desert  -,  and  (8)  Built-on  and  waste  land.  The  results  of  the 
survey  are  shown  in  Tables  3a  and  3b  in  percentages  of  the  total 
zonal  area,  separately  for  the  three  major  land  masses. 

The  roughness  parameters  for  each  season  and  all  continental 
zones  were  computed  in  the  same  manner  as  the  regional  values  for 
the  State  of  Wisconsin.  Values  of  Zo  =  0.03  cm  for  deserts,  and 
Zo  =0.1  for  snow  cover  were  adopted  from  the  study  of  Deacon  (4). 
The  results  are  listed  in  Tables  4a  and  4b.  and  illustrated  in  Fig¬ 
ures  1,  2,  and  3. 

The  latitudinal  distributions  reflect  the  characteristic  vegetation 
covers  over  the  continents.  In  general,  the  Zq -values  are  relatively 
large  for  both  the  equatorial  zone  and  temperate  to  high  latitudes: 
they  show  minimum  values  about  20  to  30  degrees  from  the  equator. 
This  is  caused  by  the  meridional  sequence  ranging  from  tropical  for¬ 
ests,  to  low  natural  vegetations  and  crops,  then  to  the  woodlands  and 
coniferous  forests  of  the  temperate  zone.  Farther  north  the  Zq -values 


Table  3  a.  Vegetation  cover  In  percent  of  total  zonal  area  —  (1)  Qrops, 

(2)  Woodland  and  tree  crops,  (3)  Meadow  and  pasture,  (4)  Forest, 

(5a)  Natural  vegetation,  (5b)  Natural  vegetation  of  generally  low  types  if 
distinct  from  5a,  (6)  Tundra,  (7)  Desert,  (8)  Built-on  and  waste  land. 


The  Americas 

(1) 

(2) 

(3) 

(4) 

(5a) 

(5b) 

(6) 

(7) 

(8) 

Total 
Area 
10*  km* 

70-60*N 

0.0 

0.0 

15.8 

30.5 

0.0 

0.0 

53.7 

0.0 

0.0 

400.2 

60-50 

7.0 

5.8 

0.8 

63.8 

21.  1 

0.0 

1.5 

0.0 

0.0 

498.0 

50-45 

18.4 

1.3 

21.9 

42,0 

14.8 

0.0 

0.0 

0.0 

1.6 

153.6 

45-40 

27.7 

0.0 

34.4 

30.2 

4.3 

0.0 

0.0 

0.0 

3.4 

213.6 

40-35 

20.5 

0.0 

37.7 

33.5 

6.9 

0.0 

0.0 

1.3 

0.  1 

236.2 

35-30 

14.8 

0.0 

40.7 

37.3 

3.1 

0.0 

0.0 

4.  1 

0.0 

202.9 

30-20 

10.0 

0.0 

31.9 

23.2 

4.6 

0.0 

0.0 

0.0 

30.3 

208.  0 

20-10 

11.3 

0.0 

11.  1 

47.  2 

7.6 

0.0 

0.0 

0.0 

22.8 

56.0 

10-Equ. 

3.8 

0.0 

12.3 

50.  1 

1.7 

0.0 

0.0 

0.0 

32.  1 

257.  4 

Equ.-10*S 

1.  6 

0.0 

6.  4 

79.  3 

8.7 

0.0 

0.0 

3.7 

0.  3 

540.9 

10-20 

2.0 

0.0 

20.2 

33.  3 

28.5 

0.0 

0.0 

4.8 

11.2 

434.  9 

20-30 

2.8 

0.0 

11.7 

48.0 

14.4 

6.  5 

0.0 

4.4 

12.  2 

291.7 

30-40 

22.9 

0.  0 

34.9 

16.  4 

13.  1 

1.3 

0.0 

0.0 

11.  4 

164.9 

40-50 

0.0 

0.0 

56.3 

30.0 

10.0 

0.0 

0.0 

0.0 

2.  8 

70.7 

50-55 

0.0 

0.0 

24.9 

69.8 

5.3 

0.0 

0.0 

0.0 

0.  0 

22.  8 

52 


Table  3b.  Vegetation  cover  in  percent  of  total  zonal  area  —  (1)  Crops, 

(2)  Woodland  and  tree  crops,  (3)  Meadow  and  pasture,  (4)  Forest, 

(5a)  Natural  vegetation,  (5b)  Natural  vegetation  of  generally  low  types  if 
distinct  from  5a,  (6)  Tundra,  (7)  Desert,  (8)  Built-on  and  waste  land. 


Eurooe-Africa 

(1) 

(2) 

-(3) 

(4) 

(5a) 

(5b) 

(6) 

(7) 

(8) 

Total 

Area 

lO^km^ 

70-60*N 

5.9 

0.0 

2.7 

42.6 

0.0 

0.0 

48.8 

0.0 

0.0 

121.4 

60-50 

40.4 

0.0 

27.0 

19.  1 

0.0 

0.0 

0.0 

0.0 

13.  5 

108.0 

50-40 

40.  1 

0.0 

20.1 

25.7 

3.7 

0.0 

0.0 

0.0 

10.4 

183.2 

40-30 

22.9 

0.0 

16.4 

12.6 

6.7 

9.6 

0.0 

31.8 

0.0 

211.6 

30-20 

0.7 

0.0 

5.0 

p 

o 

0.1 

0.0 

0.0 

94.2 

0.0 

440.  9 

20-10 

2.2 

0.0 

10.8 

26.5 

2.2 

6.8 

0.0 

51.5 

0.  0 

802.  3 

10-Equ. 

14.2 

0.0 

18.8 

37.  5 

10.8 

8.4 

0.0 

10.3 

0.0 

619.  9 

Equ.-10*S 

15.4 

8.4 

1.2 

40.8 

34.2 

0.0 

0.0 

0.0 

0.  0 

333.  9 

10-20 

10.  1 

1.  6 

29.8 

40.  9 

17.6 

0.0 

0.0 

0.0 

0.0 

387.  1 

20-35 

3.  2 

0.  0 

65.4 

2.  5 

13.  5 

15.  4 

0.0 

0.0 

0.  0 

276.  0 

Asia-Oceania  (including  USSR) 


70-60‘’N 

5.  3 

0.0 

0.0 

52.  8 

0.0 

0.0 

41.9 

0.0 

0.  0 

716.  4 

60-50 

26.  9 

0.  0 

18.9 

34.2 

20.0 

0.0 

0.0 

0.  0 

0.  0 

848.  7 

50-40 

13.  2 

0.  0 

29.4 

13.  9 

7.8 

0.0 

0.0 

35.7 

0.  0 

1008.  2 

40-30 

15.  1 

0.0 

16.7 

8.8 

6.  3 

18.0 

0.  0 

35.^1 

0.0 

852.  2 

30-20 

23.  5 

0.  0 

12.  1 

11.  7 

4.9 

33.  9 

0.0 

13.  9 

0.  0 

765.  2 

20-10 

23.4 

0.  0 

7.  1 

35.  0 

8.6 

18.  3 

0.0 

7.6 

0.  0 

285.  0 

10-Equ. 

16,  6 

0.0 

0.4 

72.  6 

2.  1 

7.  5 

0.0 

0.0 

0.  8 

40.  1 

Equ.-lO^S 

8.  5 

0.  0 

0.0 

49.4 

42.  0 

0.  0 

0.  0 

0.0 

0.  1 

209.  7 

10-20 

50.8 

2.  2 

36.6 

0.  0 

0.0 

0.0 

0.0 

10.4 

0.  0 

145.  9 

20-30 

0.0 

29.9 

26.  7 

13.  3 

0.  0 

0.  0 

0.0 

30.  1 

0.  0 

419.  1 

30-40 

8.7 

2.7 

52.  7 

32.  6 

0.0 

0.0 

0.  0 

3.  3 

0.  0 

200.  3 

40-50 

2.4 

35.  4 

37.  5 

9.9 

1.  2 

13.  6 

0.0 

0.0 

0.  0 

33.  6 
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Table  4a.  Estimated  xonal  averages  of  the  aerodynamic  roughness 
parameter  (zo  $  cm)  for  the  seasons. 


The  Americas 

Spring 

Summer 

Fall 

Winter 

70-60*N 

1.  92 

8.  68 

2.  07 

1.  36 

60-50 

41.70 

141.  3 

86.74 

35.  28 

50-45 

17.86 

54.  06 

24.  97 

4.  00 

45-40 

9.  08 

i 

31.46 

19.  90 

1.28 

40-35 

16.  30 

32.  92 

19.  84 

13.  13 

35-30 

28.  29 

35.  56 

28.  86 

20.  92 

30-20 

15.  25 

16.  75 

14.  64 

12.  28 

20-10 

45.  67 

48.  87 

41.97 

48.  19 

10-Equ. 

77.  45 

79.  20 

75.  32 

78.  85 

Equ.  -lO'S 

418.9 

431.3 

416.  2 

421.6 

10-20 

44.  85 

45.  08 

44.  37 

45.  15 

20-30 

44.  03 

44.  51 

43.  57 

44.  59 

30-40 

13.  80 

10.  91 

15.  30 

17.  71 

40-50 

24.  84 

15.  53 

22.  71 

28.  80 

50-55 


126.  0 


105.4 


116.  9 


133.  1 


Table  4b.  Estimated  zonal  averages  of  the  aerodynamic  roughness 
pardmeter  (zo  i  cm)  for  the  seasons. 

Eurooe-Afrlca 

Spring 

Summer 

Fall 

Winter 

70-60* N 

4. 66 

16.  65 

11.  97 

3.  88 

60-50 

2. 19 

18. 49 

4.  48 

0.51 

50-40 

7. 48  . 

26.  24 

10.85 

2.  90 

40-30 

2.  24 

2.78 

2.31 

1.  83 

30-20 

0.  04 

0.  04 

0.  04 

0.  04 

20-10 

1.  32 

1.  34 

1.  30 

f 

1.  34 

10-Equ. 

32.  23 

35.  09 

28.  98 

34.  50 

Equ.  -10*8 

126.  1 

134.  9 

109.8 

135.  1 

10-20 

59.  50 

62.  22 

54.  29 

64.  40 

20-35  9. 32 

Asla-Oceania  (including  USSR) 

5.  37 

7.  52 

9.  46 

70-60*N 

10.  55 

29.  66 

11.  21 

9.  26 

60-50 

5.  02 

29.  42 

8.  64 

1.  87 

50-40 

0.  67 

1.  99 

1.  06 

0.  21 

40-30 

1.  00 

1.  51 

1.  21 

0.78 

30-20 

4.  08 

5.  08 

3.  69 

3.43 

20-10 

37.  60 

38.  94 

36.  82 

38.  15 

10-Equ. 

352.  6 

361.  5 

347.  4 

356.  2 

Equ.  -10°  S 

187.  2 

193.  2 

188.  3 

194.  0 

10-20 

143.  5 

100.  0 

103.  5 

145.  9 

20-30 

23.  66 

10.  37 

10.68 

10.  94 

30-40 

22.  27 

18.  93 

19.  41 

24.  62 

40-50 

47.  59 

35.  92 

43.  33 

52.  61 

55 


ZONAL  yg/uis  Of  wommutt  of  vchtatiom  coycji 


ARCA-WEMHTCO  LOe.-MCAN  -  COMPUTCO  OY  C.  KUNO,  FMM 
UNITED  NATIONS  YEANBOOK  OF  FOOD  N  AONIC.  STATISTICS  FOB  ISS7 


Fig.  2.  Meridional  profiles  of  the  average  aerodynamic 
roughness  parameter  derived  from  vegetation  cover  on  the 
continents  in  summer. 
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decrease  again,  especially  in  North  America,  due  to  lower  vegetation 
and  extended  tundra  regions.  i 

In  the  southern  hemisphere,  the  latitudinal  distributions  of  zo 
are  nearly  the  same  for  all  three  continental  groups.  In  the  northern 
hemisphere,  there  are  significant  differences.  In  the  low  and  middle 
latitudes  of  Europe  and  North  Africa,  extremely  small  Zo  -values  are 
due  to  extended  low-vegetation  and  desert  areas,  while  lush  vegeta¬ 
tion  in  the  same  latitudes  of  North  America  and  Asia  causes  relatively 
high  values.  At  highest  latitudes  of  Europe  and  Asia  the  roughness 
parameter  is  significantly  larger  than  in  North  America  owing  to  the 
differences  in  northern  extent  of  coniferous  and  other  forests. 

The  seasonal  trends  of  zonal  roughness  parameters  reflect  the 
*  phenological  cycle.  There  is  little  change  in  tropical  regions.  Crop 
growth  and  phenological  sequences  of  vegetation  height  in  middle  and 
high  latitudes  cause  roughness  maxima  in  summer  and  minima  in  win¬ 
ter.  It  can  be  noted  that  the  autumn-wilting  of  low-type  vegetations, 
and  winter  snow  covering  them,  affects  significantly  the  seasonal 
changes  in  roughness  parameter. 

5.5  Discussion,  and  Concluding  Remarks 

The  evaluation  of  regional  and  zonal  roughness  parameters  is 
based  on  several  assumptions,  and  the  results  must  be  considered 
tentative  and  restricted.  In  the  survey  of  land  use  and  vegetation 
cover,  a  great  number  of  different  sources  were  employed.  Since  the 
statistical  data  were  frequently  found  incomplete,  certain  adjustments 
had  to  be  made.  Extrapolations,  supported  by  economical  maps,  were 
necessary,  especially  for  Australia,  South  America,  and  the  USSR. 

Vegetation  height  was  considered  the  major  factor  in  causing 
aerodynamical  roughness.  The  forests,  therefore,  play  an  important 
role.  It  must  be  expected  that  the  density  of  forest  trees  is  an  addi¬ 
tional  factor.  Due  to  the  lack  of  pertinent  information  this  factor 
could  not  be  considered  here.  The  problem  can  be  solved  only  by 
detailed  experimental  investigations  in  relatively  small  areas.  Simi¬ 
lar  restrictions  are  imposed  by  the  adoption  of  Deacon' s  zq  values 
over  desert  and  snow.  However,  it  is  believed  that  only  relatively 
minor  revisions  of  the  zonal  distribution  of  aerodynamic  roughness 
parameters  will  be  necessary. 

The  method  of  computing  the  representative  area-means  of  rough¬ 
ness  parameters  as  the  weighted  average  of  (log  Zo)-values  is  sug¬ 
gested  by  Lettau'  s  (21)  findings  that  the  logarithm  of  the  surface- 
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Rossby  number  (which  is  proportional  to  l/zo)  determines  the  drag 
coefficient  of  large-scale  atmospheric  currents.  The  meteorological 
and  climatological  significance  of  the  estimates  of  regional  and  zonal 
roughness  parameters  can  only  be  tested  by  a  comparison  of  actual, 
with  computed  friction  effects  in  the  lower  atmosphere.  A  relatively 
direct  and  readily  available  climatological  measure  of  actual  friction 
effects  is  the  mean  ratio  of  observed  wind  speed  (at  anemometer  level 
of  a  synoptic  or  climatological  station)  to  regional  geostrophic  speed; 
reference  is  made  to  Lettau  (21).  A  continent-wide  study  of  this  ratio, 
together  with  an  evaluation  of  vegetive  roughness  parameters  and 
subsequent  drag  coefficients,  is  in  progress.  The  investigation  will 
include  a  study  of  the  combined  effects  of  plant  cover  and  topography 
on  the  aerodynamic  roughness  of  continent  surfaces. 
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A  Method  for  Machine  Computation  of 
Wind  Profile  Parameters 
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Abstract.  A  least-square  method  is  described  by  which 
the  three  parameters  of  the  wind  profile  in  the  adiabatic 
surface  layer  of  the  atmosphere  can  be  determined  with 
the  aid  of  electronic  digital  computers. 

6.1  Introduction 


Under  adiabatic  conditions,  the  wind  profile  in  the  atmospheric 
surface  layer  is  determined  by  three  parameters,  which  are  the  rough 
ness  length,  the  zero  point  displacement,  and  the  surface  stress. 
When  reference  is  made  to  Lettau  [  1  ],  the  three  parameters  are  re¬ 
lated  by  the  following  equation: 

i 

V  =  )c"Mto  /p)^  In  [  (z  +  d  +  zo  )/zo  ]  (1) 

where 

k  =  K^rm^n  number  =  0.38 
To  =  surface  stress  (dyne/cm^) 
p  =  air  density  (g/cm^) 
z  =  nominal  height  of  anemometer  (cm) 
d  =  zero  point  displacement  (cm) 

Zo  =  roughness  length  (cm) 

V  =  mean  wind  speed  (cm/sec)  at  height  z. 
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Let  V*  =  k’^To/p)^  *  friction  velocity,  (cm/sec"*) , 

D  -  d  -t*  zo  »  (cm)  ; 

then,  equation  (1)  can  be  rewritten  as: 

V  =  v*ln((2  +  D)Ao]  .  (2) 

Given  at  least  three  observations  of  V  and  z ,  we  wish  to  solve 
(2)  for  V*,  D,  and  zo  •  We  therefore  subscript  the  observed  vari¬ 
ables  V  and  z ,  and  seek  a  least-squares  solution  to 

Vj  =  V*  In  ( (Zj  +  D)/zo]  .  1  <  1  <  N  .  (3) 

where  N  denotes  the  number  of  observation  levels.  The  least-squares 
method  requires  that 

f  e{  ■  f  {V  -V*ln[(z  +D)/z,]}* 

1*1  *  1=1  ‘  ‘ 

be  minimized. 


The  derivation  to  follow  will  be  greatly  simplified  if  we  use 
vector  notation.  A  vector  will  be  simply  an  ordered  array  of  N  ele¬ 
ments;  considerations  of  magnitude  and  direction  will  not  enter 
directly.  Vector  quantities  will  be  distinguished  by  a  bar  (e.  g. , 

V,  F,  G),  and  their  scalar  product  will  be  indicated  by  a  dot: 

-  N 

F-  G  =  Z  (FjG  )  , 

i=l 

where  and  G^  are  arbitrarily  defined  arrays  of  N  elements  eacJk.. 
Let  __ 

W  =  {In  Zo  ,  In  Zo  ,  . . . ,  In  Zo }  , 

X  =  (In  (zi  +  D) ,  In  (Z2  +  D) ,  ....  In  (z  +  D)}  , 

—  N 

V  =  {V,,  Vz . Vj^}  . 

B=  {(2, +  Dr',  (Zz  +  D)-',  ....  (Zjj  +  D)-*}. 

T  =  {1.  1 . 1}  . 

Then  (3)  becomes 

V  =  V*(X  -  W)  , 


(4) 
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and  the  vector,  of  the  values  Is 

r  =  V- V*(X- W)  .  (5) 

Thus, 

N  _ _  _ _  _ 

Z  =  £  •  r  =  V*  V  -  2V*V*  (X-W)  +  V**(X-W)-(X-W).  (6) 

The  condition  that  the  quantity  (6)  be  minimized  is  that  the  partial 
derivatives  of  (6)  with  respect  to  V*,  D,  and  Zo  be  zero.  Remem¬ 
bering  that 

illl  ::  /I  J  illL 

Bzo  'zo  '  81n  Zo  * 

we  have 


^  =  (^^)[2V*VI-2V**(X-W)-r]  =0.  (7) 

or 

v-T  =  V*  (x-w)-T  .  (8) 

If  we  let  Xg  =  N-'{X-7.  j(-T . X-T}.  which  is  the  vector  of 

means  of  the  elements  of  X,  and  use  a  similar  convention  for  V 

and  B^,  we  may  write  the  vectors  of  adjusted  values  of  V,  X,  and 
B,  as 


V  =  V  -  V 

a 

X  =  X -X 

a 

b  =  B  -  B  . 
a 

Now,  by  (8)  we  have 

W  =  W=X-V*  V 
a  a  a 

Hence,  from  (5), 

7  =  ( V- V  )  -  V*(X-X  )  =  V  -  V*x  , 
a  a 

and 


(9) 


(10) 

(11a) 


Since 


£  •  £  =  vv  -  2V’*‘ v •  X  +  •  X 

9(  £  •  £  )/  8V*  =  0,  V*=(v*x)/(x'x), 


(11b) 


(12) 


66 


Stephen  M.  Robinson 


and  r-r  =  V- V  -  (v‘x)V(x-x)  .  (13) 

Applying  the  final  condition,  8(7*7  )/8D  =  0,  we  see  that  since 
V*  X  should  never  be  zero, 

g(D)  =  (v*x)  (x*b)  -  (x*x)  (vb)  =  0  .  (14) 

The  term  g(D)  in  equation  (14)  is  a  function  of  one  unlcnown  only 
which  is  D.  Equation  (14)  can  be  attaclced  by  conventional  iteration 
techniques.  Perhaps  the  easiest  of  these  to  use  —  because  it  requires 
no  derivatives  —  is  the  secant  method  of  Jeeves  [  Z] .  Given  two  suc¬ 
cessive  approximations  to  D,  say  D  ,  and  D  ,  this  method  defines 

n-i  n 

D-  =(D  ,g(D)-Dg(D  ,)] /[g(D  )  -  g(D  ,)]  .  (15) 

n+1  ^  n-1  n  n  n-1  ^  n  n-l  ^ 

The  secant  method  is  of  convergence  order  1.62  .  This  means  that 
the  increase  in  the  number  of  significant  figures  at  each  step  is  ap¬ 
proximately  1.62  times  the  previous  increase.  The  secant  method 
has  been  found  most  effective  in  practice.  The  unknown  D  may  be 
determined  to  as  close  a  tolerance  as  one  wishes;  after  doing  so, 
one  finds  V*  and  Zo  by  substitution,  using  equations  (10)  and  (12). 

6.3  Sample  Computation 


Some  wind  data  for  an  adiabatic  surface  layer  given  by  Lettau  [1, 
p.  333]  have  been  used  to  compute  a  sample  set  of  wind  profile 
parameters.  The  data  are: 


i 

(cm.  sec’* ) 

Zj  (cm) 

1 

837.  5 

640 

2 

759.2 

320 

3 

673.7 

160 

4 

585.  3 

80 

5 

490.2 

40 

The  computation  proceeds  as  follows:  Let  us  choose  two  initial 
guesses  for  D-say,  -20  and  -15  cm  —  and  call  them  D.i  and  Do  . 
Then,  machine  computation  resulted  in  the  following  sequence: 
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Di  =  -12.684 
Di  =  -10.636 
D,  =  -  9.  779 
D4  =  -  9.  554 
Ds  =  -  9.534 
Dfc  =  -  9.  534  . 

Thus,  within  six  steps  the  Iteration  has  converged  to  a  final  D-value 
within  a  tolerance  of  0.001  cm.  Upon  determining  the  values  of 
and  Zo  from  (10)  and  (12) ,  respectively, 

D  =  -9.  534  cm 
V***  =  115.  1  cm/sec 
Zo  =  0.433  cm. 

Using  hand  computation  Lettau  [  1,  p.  335]  had  obtained  for  the 
same  sample  profile  0=  -9.  5  cm,  and  Zo  =  0.442  cm. 

It  should  be  noted  here  that  the  value  of  Zo  ,  being  an  exponential 
function,  Is  very  much  subject  to  error.  A  relatively  small  variation 
In  D  and/or  V***  can  produce  a  large  variation  in  Zo .  One  must  there¬ 
fore  be  careful  when  doing  the  bade  substitution. 

6.4  Windprofile  Program 

The  program  (see  Appendix)  will  read  in  two  estimates  to  D,  along 
with  a  decimal  tolerance  T  >  0 ,  and  N  pairs  of  data,  each  pair  con¬ 
sisting  of  a  level  zi  and  a  wind  speed  Vi.  After  reading  these,  the 
program  computes  and  prints  the  values  of  V^r  D,  and  zo  necessary 
to  produce  a  least-squares  fit,  along  with  the  final  error  sum  of 
squares  (R^ ).  The  program  is  written  in  the  IBM  FORTRAN  programming 
language,  and  can  be  compiled  and  run  on  any  computer  able  to  use 
that  language. 

Special  features  are  as  follows:  (a)  If  the  program  encounters  a 
non-positive  argument  of  a  logarithm  function,  or  if  it  attempts  to 
divide  by  zero,  the  computer  will  print: 

ERROR  AT  ’"tH  ITERATION.  TRY  OTHER  INITIAL  GUESSES. 
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The  asterisk  will  be  replaced  by  a  number  indicating  the  number  of 
successful  Iterations  made  before  the  error  occurred.  The  next  set 
of  data  will  be  read  Immediately,  (b)  If  sense  switch  1  (or  its  equiv¬ 
alent)  is  up,  the  value  of  D  at  each  iteration  will  be  printed,  along 
with  the  number  of  that  iteration.  This  feature  is  useful  for  keeping 
track  of  the  progress  of  the  solution. 

The  card  format  (for  each  set  of  data)  is  as.  follows: 


Card  1: 


Cols.  1-6  ±xx.  XX 
7-12  ±xx.xx 
13-20  ±  .xxxxxx 


(the  1st  estimate  to  D,  or  D.i) 
(the  2nd  estimate  to  D,  or  Do) 
(the  value  of  the  tolerance  T ) 


Cards  2  through  (N  +  l): 

Cols.  1-8  +  XXX.  XXX 

9-16  +  XXX.  XXX 


Card  (N  +  2): 

Cols.  1-8  -999. 999 

9-16  -999.999 


(wind  speed  V^) 
(anemometer  height  Zj) 


(dummy  variables  for 
control  purposes) 


To  run  the  program,  simply  load  the  compiled  deck  followed  by 
the  data  cards.  The  printing  will  take  place  automatically. 

Answer  format  (all  answers  are  in  floating-point  form): 

V*  D  Zo  R^ 
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Appendix 


PROGRAM  WPRF 

C  PROGRAM  TO  FIND  WINDPROFILES  BY  ITERATION 

C  1  MARCH  1961  STEPHEN  M.  ROBINSON 

C 

DIMENSION  V(50),  Z(50) 

3  FORMAT  (2F6.  2,  F8.  6) 

4  FORMAT  (2F8.  3) 

5  FORMAT  (I  3,  4E15.  6) 

6  FORMAT  (lOH  ERROR  AT  14,  39HTH  ITERATION.  TRY 

6  C  OTHER  INITIAL  GUESSES) 

7  FORMAT  (El 5.  6,  13) 

ER=(10.  )**(-5) 

20  N:0 

SV=0 
S2V=0 

READ  3,  DO.  DN.  T 
J=0 

150  READ  4,  Ul,  U2 

IF  (Ul)  10,  11,  11 
11  N=N+1 

Z(N)=U2 
V(N)=U1 
SV=SV+U1 
S2V=S2V+U1**2 
GO  TO  150 
10  AN=N 

S=0 
D=DO 

13  SX=0 

S2X=0 
SVX=0 
SB=0 
SVB=0 
SXB=0 

DO  14  1=1,  N 
Y=(Z(I)+D) 

IF  (Y)  23,  23,  102 

23  PRINT  6,  J 

GO  TO'20 
102  X=LOGF(Y) 

SX=SX+X 
S2X=S2X+X«*2 
SVX=SVX+X*V(  I ) 
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B=l.  A 
SB>SB+B 

SVB«SVB+B*V(1) 

14  SXBsSXB+X«B 

G»(SVX-(SV*SX)/AN)*(SXB-(SX*SB)/AN)-(S2X- 
C  (SX**2)/AN)*(SVB-(SV*SB)/AN) 

IF  (S)  15, 16, 15 

16  S=l. 

G0=G 

19  D=DN 

GO  TO  13 

15  GN^ 

M+1 

El  :  (D0*GN-G0*DN) 

E2  :  (GN-GO) 

22  DP  =  E1A2 

IF  (ABSF(DP-DN)-T)  17,  17, 18 
18  DO  =  DN 

G0=GN 
DN  =  DP 

IF  (SENSE  SWITCH  1)  100,  101 

100  PRINT  7  DP,J 

101  GO  TO  19 

17  VST  =  (SVX-(SV*SX)/AN)/(S2X-(SX**2)/AN) 
ZNT:«XPF((SX-{  1.  AST)*SV)/AN) 
R2=(S2V-(SV**2)/AN)-VSTe((2.*SVX-VST*S2X)- 
SX*(2.  ♦SV-VST*SX)/AN) 

PRINT  5,  J,  VST,  DP.  ZNT,  R2 

GO  TO  20 

END 

END 
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Investigations  of  the  Modification  of  Wind  Profiles 
by  Artificially  Controlled  Surface  Roughness^ 


John  E.  Kutzbach 

Department  of  Meteorology 
University  of  Wisconsin 


Abstract.  A  series  of  controlled  wind  profile  experiments 
were  conducted  on  the  ice  of  Lake  Mendota  under  condi¬ 
tions  of  near  neutral  thermal  stratification.  The  logarith¬ 
mic  model  for  the  variation  in  wind  speed  and  height  is 
assumed.  Changes  in  the  wind  profile  parameters  for 
controlled  variations  in  the  area  density  (number  per  unit 
area)  of  ordinary  bushel  baskets  are  examined.  The  den¬ 
sity  of  the  obstacles  ranged  from  one  basket  per  48.5  m^ 
to  one  basket  per  0.4  m^.  Within  this  range,  the  aero¬ 
dynamic  roughness  parameter  was  proportional  to  basket 
density  raised  to  the  1.1th  power,  the  zero  point  displace¬ 
ment  was  proportional  to  basket  density  to  the  0.29th 
power,  and  the  surface  stress  was  proportional  to  basket 
density  to  the  0.49th  power.  These  results  are  compared 
with  similar  experiments  conducted  in  ducts  and  wind 
tunnels.  The  dependence  of  the  geostrophic  drag  coeffi¬ 
cient  on  the  Surface -Rossby  Number  is  discussed. 

7.1  Symbols  and  Definitions 


7.1. 1  Geometric  characteristics  of  the  obstacles  and  their  distribution 
h  =  height  of  bushel  basket,  or  in  general,  obstacle  height  (cm) 
S  =  average  spacing  distance  between  bushel  basket  centers  (m) 


*Part  of  this  work  was  submitted  to  the  University  of  Wisconsin 
in  partial  fulfillment  of  the  requirements  for  the  degree  of  Master  of 
Science. 
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j.  1.  Lake  Mendota.  Madison,  Wisconsin showinQ  the  location  of  the 
partment  of  Meteorology.  University  of  Wisconsin,  micrometeorological 
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S.  A.  s  specific  area  of  obstacle  (m^ )  s  ( S  =  reciprocal  of  area 
density  of  obstacles 

A  =  non-dimensional  ratio  of  specific  area  over  the  lateral  sil¬ 
houette  area  of  the  roughness  element 

R  =  upwind  roughness  fetch  (m)  =  length  of  upwind  trajectory 
over  roughness  field 

r  =  downwind  roughness  fetch  (m)  s  length  of  downwind  (or 
crosswind)  trajectory  over  roughness  field 

^  -  angular  width  of  roughness  field  (degrees) 

X  =  distance  downwind  from  leading  edge  of  roughness  field  (m) . 

7.1.2  Dimensional  variables  and  wind  profile 'parameters 

V(z)  =  mean  velocity  profile  (cm/sec) 

V  =  mean  velocity  at  control  anemometer  (cm/sec) 
con 

To  =  surface  stress  (dynes/cm^ ) 
p  =  air  density  (gm/cm^ ) 

K  =  K^m^n  constant 

i 

V’*'  =  k“*(to/p)*  =  friction  velocity  (cm/sec) 
z  =  nominal  height  (cm) 
zo  =  aerodynamic  roughness  parameter  (cm) 
d  =  zero  point  displacement  (cm) 


7.1.3  Flow  parameters  and  stability  parameters 


C  =  (to/p)^/V  =  geostrophic  drag  coefficient 

g,o 

Roo  =  V  /zo  f  =  Surface  Rossby  Number 
—  g»o 

Ri  =  (g/T  )  Oe/az)/OV/8z)^  =  g  AzAe/T„{AVF 
m 

=  Richardson  number  referred  to  the  geometric  mean  height 


z(z  +  Az)  of  two  observation  levels  z  and  z  +  Az 


Fig.  2.  Schematic  illustration  of  the  shape  and  size  of  the  rough¬ 
ness  field,  the  dimensions  of  the  obstacles,  and  the  sequence  of 
obstacle  densities  (drawn  to  scale)  of  wind  profile  experiments  on 
the  ice  of  Lake  Mendota,  1961. 


—  denotes  first  basket  array, 

S.  A.  =  14.7  m 

©  +© 

—  denotes  second  basket  array. 

S.  A.  =  7.7  m 

©  +©  +  • 

—  denotes  third  basket  array, 

S.  A.  =  4.0  m 

% 
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(Ri)'  =  Yj  Ri(zJ/E  stability  parameter  (l/m) 

i  ^  i  * 

7.1.4  Additional  variables  from  duct  flow 

k  =  equivalent  height  of  sand  roughness  elements  (cm) 
s 

N  s  h/zo  (or  k  /zo  )  s  Nikuradse'  s  roughness  ratio 
s 

i 

V*  =  '  shearing  velocity  at  rough  wall  (cm/sec) 

#  i 

V  =  (t  /pr  “  shearing  velocity  at  smooth  wall  (cm/sec) 

7.1.5  Additional  variables  from  atmospheric  flow 

V  s  geostrophic  wind  speed  at  surface  (cm/sec) 

giO 

u  =  horizontal  component  of  mean  velocity  (cm/sec) 

w  =  vertical  component  of  mean  velocity  (cm/sec) 

or  =  direction  of  surface  wind  (degrees) 

a  =  direction  of  geostrophic  wind  at  the  surface  (degrees) 
g.o 

7.2  Introduction 


The  functional  relationship  between  the  physical  and  vegetational 
structure  of  the  earth/air  interface  and  the  basic  parameters  of  the 
low-level  wind  profile  is  an  important  meteorological  problem.  The 
pertinent  parameters  of  the  problem  (roughness  length,  zero  point  dis¬ 
placement,  and  surface  stress)  have  been  measured  at  a  variety  of 
micrometeorological  sites  under  diverse  atmospheric  conditions.  This 
information  has  limited  usefulness  for  the  physical  understanding  of 
the  processes  involved  because  the  configuration  of  natural  roughness 
elements  ’is  usually  poorly  defined. 

Lettau  (1959)  has  outlined  the  concepts  of  "passive"  and  "active" 
control  in  micrometeorological  field  experiments.  He  defines  the 
terms  in  the  following  manner:  ".  .  .  passive  control.  .  .  means  that 
the  field  experiment  is  made  over  a  natural  surface  which  is  carefully 
chosen  according  to  its  physical  structure,  and  scheduled  at  times 
when  the  weather  conditions  are  defined  and  satisfy  certain  predeter¬ 
mined  principles.  .  .  .Active  control  requires  a  deliberately  planned 
modification  of  at  least  one  of  the  external  parameters  of  physical 


Fig.  3.  The  10-mlnute  mean  wind  profiles  associated  with  the  basket 
distributions  of  Fig.  2.  Solid  lines  and  circles  indicate  increasin<^ 
obstacle  density  -  broken  lines  and  squares  indicate  decreasing  ob¬ 
stacle  density. 


O  ""  denotes  profile  over  unmodified  ice 

®  -  denotes  profile  over  first  basket  array,  S.  A. 

©  —  denotes  profile  over  second  basket  array,  S.  A. 

#  —  denotes  profile  over  third  basket  array,  S.  A. 


14.  7  m' 
7.  7  m' 
4.  0  m' 
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process  In  the  lower  atmosphere."  Possibilities  of  exerting  some  de¬ 
gree  of  control  on  natural  processes  In  the  lower  atmosphere  deserve 
more  attention  In  meteorological  research. 

The  primary  purpose  of  this  paper  Is  to  examine  changes  In  the 
wind  profile  parameters  for  controlled  variations  In  the  area  density 
(number  per  unit  area)  of  well-defined  roughness  elemehts.  The  na¬ 
ture  of  the  experiments  is  discussed  In  Section  ,7. 3.  Section  7.4 
deals  with  the  effect  of  obstacle  density  (number  per  unit  area)  on 
the  wind  profile  parameters  and  the  dependence  of  the  geos  trophic 
drag  coefficient  on  the  Surface  Rossby  Number. 

Some  supplementary  results  are  presented  In  Section  7.5.  These 
Include  an  experimental  determination  of  the  Kdrmdn  constant,  the 
vertical  motion  pattern  over  the  roughness  field,  the  variation  In 
wind  speed  In  the  air  layer  between  the  obstacles,  and  an  estimate 
of  the  individual  drag  coefficient  of  the  obstacle. 

Section  7.6  contains  an  analysis  of  possible  errors  and  a  discus¬ 
sion  of  data  evaluation.  Particular  attention  Is  given  to  the  growth  of 
the  "internal  boundary  layer"  downwind  over  the  roughness  field. 
Conclusions  are  presented  in  Section  7.7. 

7.3  The  Nature  of  the  Experiments 

7.3.1  Passive  control 

The  smooth,  frozen  surface  of  Lake  Mendota  at  Madison,  Wis¬ 
consin,  is  of  sufficient  extent  (approximately  40  km^  surface  area)  to 
allow  the  establishment  of  a  representative,  low-level  wind  profile. 
Lettau  (1959)  has  suggested  as  a  "rule  of  thumb"  that  the  highest  in¬ 
strument  level  should  not  exceed  l/50th  of  the  distance  from  signifi¬ 
cant  discontinuities  in  the  surface  structure.  The  highest  anemometei 
level  was  3.4  m  above  the  ice  and  the  upwind  fetch  over  the  unob¬ 
structed  ice  was  at  least  2000  m  during  all  experiments  at  the  micro- 
meteorological  site  located  400  m  off  Second  Point  (see  Fig.  1). 
Therefore,  Lettau' s  suggested  criterion  for  representative  profiles  is 
more  than  satisfied. 

Each  experiment  consisted  of  a  pre-planned  sequence  of  6  to  10 
ten-minute  mean  wind-profiles  extending  over  a  2-  to  3 -hour  period. 
Synoptic  conditions  were  chosen  such  that  moderate  to  strong,  steady 
winds  could  be  expected  for  several  hours. 

Ideally,  the  investigation  should  be  restricted  to  cases  of  neutral 
thermal  stratification  (adiabatic  lapse  rate).  In  practice,  it  must  be 


PROFILES- BUSHEL  BASKET  EXPERIMENT 
ft>-MINUTE  KIMS.  1961  JANUAMfV  27.  a<29-»40 


1961.  The  controlled  obstacle  parameter  was  the  specific  area  (S.  A.  ); 
see  Fig.  8  for  illustration  of  obstacle  array.  Only  one  run  was  made  at 
each  obstacle  density.  No  control  anemometer  was  used;  therefore,  wind 
speeds  are  presented  as  the  ratio  of  the  wind  speed  at  a  given  level  to  the 
wind  speed  at  150  cm. 
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realized  that  adiabatic  lapse  rates  occur  in  the  atmospheric  surface 
layer  as  transient  states  only.  Slightly  diabatic  conditions  were 
permitted.  Richardson  number  requirements  will  be  discussed  in 
Section  7.6.2. 

The  turbulence  produced  by  an  obstacle  may  be  associated  with 
mechanical  effects,  thermodynamic  processes  due  to  differential 
heating,  or  both.  The  effect  of  differential  heating  has  been  mini- 
mized  in  these  experiments  by  the  presence  of  generally  overcast 
skies. 

7.3.2  Active  control 

Approximately  500  commercial  bushel  baskets  placed  with  their 
open  ends  down  were  used  for  the  controlled  modification  of  rough¬ 
ness  at  the  air-ice  interface.  It  was  necessary  to  remove  the  handles 
(leaving  1  cm  stubs)  to  prevent  the  baskets  from  sliding  on  the  ice. 
Three  or  four  people  distributed  the  baskets  by  hand. 

In  a  typical  experiment,  an  initial  number  (approximately  100)  of 
baskets  were  distributed  in  a  fan-shaped  area  (of  approximately  60* 
angular  width)  upwind  of  the  main  anemometer  .mast.  The  upwind  ex¬ 
tent  of  the  roughness  fieid'(R)  was  usually  50  m  or  less.  To  eliminate 
any  edge  effect  two  to  four  rows  of  baskets  were  placed  completely 
around  the  main  anemometer  mast.  The  outside  radius  (r)  of  this 
small  circle  of  baskets  was  approximately  5  m.  For  a  given  area,  the 
density  of  the  baskets  was  systematically  doubled  two  (or  three) 
times.  The  general  shape  and  size  of  the  array,  the  dimensions  of 
the  baskets,  and  their  distribution  are  illustrated  in  Fig.  2. 

Ten-minute  mean  wind  and  temperature  profiles  were  obtained 
before  and  after  each  surface  modification  as  well  as  over  the  ice  at 
the  beginning  and  end  of  each  experiment.  This  was  done  to  insure 
reproducibility  of  the  results,  and  to  eliminate  trends  or  over -all 
variations  due  to  large-scale  meteorological  factors.  Wind  speeds 
are  presented  as  the  ratio  of  the  wind  speed  at  a  given  level  to  the 
wind  speed  at  a  control  anemometer  located  at  the  3.4  m  level  approx¬ 
imately  10  m  upwind  of  the  roughness  field. 

Fig.  3  illustrates  the  wind  profiles  associated  with  the  basket 
distributions  of  Fig.  2.  There  is  a  consistent  systematic  change  in 
wind  shear  as  obstacle  density  is  first  increased  (solid  lines)  and 
then  decreased  (broken  lines).  All  wind  profile  parameter  computa¬ 
tions  are  based  on  the  average  of  the  two  wind  profiles  obtained  for 
each  obstacle  density. 
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The  results  discussed  in  the  following  sections  refer  to  observa¬ 
tions  obtained  on  five  different  days  in  the  period  January  through 
March  1961.  Since  the  experimental  procedure  and  the  external  con¬ 
ditions  varied  somewhat  from  day  to  day,  the  date  of  the  experiment 
is  retained  for  purposes  of  identification  on  all  graphs  and  tables. 

The  wind  profiles  for  the  five  experiments  are  presented  in  Figs.  3, 

4,  5,  6,  and  7.  Each  experiment  is  briefly  described  in  the  caption 
and  the  values  of  the  controlled  parameter,  the  .wind  speed  at  the 
control  anemometer,  and  the  computed  values  of  zero  point  displace¬ 
ment  and  aerodynamic  roughness  parameter  are  summarized  in  the  key. 

The  range  of  obstacle  densities  used  in  the  study  is  illustrated 
in  Fig.  8.  To  avoid  the  continual  use  of  fractions,  it  is  convenient 
to  refer  to  "basket  specific  area"  S.  A.  (m^ ).  As  illustrated,  S.  A. 
ranged  from  48.5  to  0.4  m^  which  is  equivalent  to  the  statement  that 
the  obstacle  density  ranged  from  one  basket  per  48.5  m^  to  one  bas¬ 
ket  per  0.4  m^ . 

Attention  is  called  to  the  variation  in  upwind  fetch  (R)  over  the 
roughness  field  (see  Fig.  8).  This  was  unavoidable  because  of  the 
limited  number  of  baskets  available.  The  possibility  exists  that  the 
shortest  roughness  fetch  (R  =  18  m)  was  not  sufficient  for  the  estab¬ 
lishment  of  a  representative  wind  profile.  This  point  will  be  discussed 
in  the  section  on  error  analysis. 

Detailed  quantitative  data  on  basket  distribution,  computed  wind 
profile  parameters,  and  Richardson  numbers  are  summarized  in  Table  3 
in  the  Appendix. 

7.3.3  Instrumentation 

It  was  anticipated  that  conclusive  results  could  only  be  obtained 
with  a  sufficiently  large  number  of  anemometers,  spaced  more  densely 
in  the  vertical  than  is  customary  for  conventional  micrometeorological 
observations.  Therefore,  the  eleven  anemometers  were  located  at 
levels  of  10,  20,  30,  40,  60,  80,  100,  160,  180,  320,  and  340  cm 
on  the  main  anemometer  mast.  The  relative  spacing  of  the  anemometer 
levels, (determined  within  a  tolerance  of  ±1  mm)  was  fixed;  however, 
the  entire  mast  could  be  moved  vertically  —  its  usual  zero  position  be¬ 
ing  approximately  3  cm  above  the  ice.  As  a  result,  the  nominal  height 
z  (i.  e. ,  the  measured  distance  from  the  ice  surface  to  the  center  of 
the  anemometer  cups)  is  obtained  by  adding  3  cm  to  each  of  the  mast 
levels.  A  remote  anemometer  at  a  height  of  340  cm  was  used  for  con¬ 
trol  purposes. 


The  controlled  obstacle  parameter  was  the  roughness  fetch  (R).  The  specific  area 
(S.  A.)  was  kept  constant  at  2.0  . 
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The  anemometer  (the  3-cup  rotor  variety  with  cups  of  light-weight 
molded  plastic,  reinforced  with  aluminum  rings)  were  designed  and 
manufactured  by  C.  W.  Thornthwaite  Associates.  Shielded,  aspirated, 
thermocouple  probes  at  20.  40.  80,  160,  and  320  cm  above  the  ice 
provided  the  data  for  the  temperature  profiles.  Reference  is  made  to 
Steams  (1961)  for  a  detailed  description  of  the  Department  of  Meteor¬ 
ology.  University  of  Wisconsin,  micrometeorological  installation  on 
Lake  Mendota. 

7.4  Effect  of  Variations  in  Obstacle  Density  on  the  Wind  Profile 
Parameters 


;  It  is  well  known  that  under  adiabatic  conditions  in  the  atmospheric 
surface  layer  the  mean  wind  profile  is  satisfactorily  determined  by  the 
aerodynamic  roughness  parameter  (zo ),  the  zero  point  displacement 
(d),  and  the  surface  stress  (tq  ).  The  three  parameters  are  related  by 
the  following  equation: 

1 

V(z)  =  k“*(to/p)*  In  [(z  +  d  +  Zo)/zo] 

(1) 


=  V*  In  [  (z  +  d  +  Zo  )/zo] 


where  k  is  the  Kirmdn  constant  and  V***  is  the  friction  velocity. 

Lettau  (1952,  pp.  60  and  75;  1957,  p.  334)  discusses  the  neces¬ 
sity  of  introducing  the  zero  point  displacement  because  of  the  unlike¬ 
lihood  that  the  true  reference  level  of  the  logarithmic  model  of  the 
wind  profile  coincides  exactly  with  the  nominal  reference  level  (from 
which  z  is  measured)  immediately  below  the  anemometers.  Notice 
that  the  boundary  condition  is  that  V  =  0  where  z  =  -d . 

A  least-square  error  technique  for  the  determination  of  the  three 
wind  profile  parameters  on  electronic  digital  computers  is  described 
in  Section  6  by  Robinson  and  has  been  used  in  the  analysis  of  all 
wind  profile  data.  The  method  is  essentially  the  same  as  that  de¬ 
scribed  by  Lettau  (1957).  The  computer  prints  the  values  of  V*, 
d  +  Zo  ,  and  Zo  necessary  to  produce  a  least-squares  fit  for  equation 
(1),  along  with  the  final  sum  of  squared  errors.  This  information  is 
presented  in  tabular  form  in  the  Appendix  (see  Table  3). 

The  computed  values  of  Zo ,  d,  and  tq  will  be  presented  graphic¬ 
ally  as  functions  of  a  non-dimensional  variable  (A)  defined  as  the 
ratio  specific  area  (S.  A.)  over  lateral  silhouette  area.  The  latter  con¬ 
cept  is  used  in  aerodynamic  drag  studies  and  is  defined  as  the  frontal 


WIND  PROFOS- BUSHEL  BASKET  EXPERIMENT 
m  OHiiiun  runi.  itn  marcm  m.  get 


Fig.  7.  Wind  profile  experiment  on  the  ice  of  Lake  Mendota, 
March  18,  1961.  The  obstacle  height  (h)  was  doubled  by  in¬ 
verting  the  lower  basket  and  placing  another  basket  on  top  of 
it.  The  roughness  fetch  (R)  and  the  specific  area  (S.  A.)  were 
kept  constant  at  R  =  26  meters  and  S.  A.  =  4.0  m^  . 
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(crosswindi  or  lateral)  area  exposed  by  the  body  to  the  flow 
(Schlichtlng,  1955,  p.  15).  Because  of  the  near  equality  between 
the  silhouette  area  when  the  basket  is  viewed  laterally  (0. 1185  m^), 
and  from  above  (0.1385  m^),  the  quantity  lOO/A  can  also  be  inter¬ 
preted  as  the  percentage  of  the  horizontal  roughness  field  covered  by 
obstacles.  The  area  covered  by  baskets  ranged  from  0.25  to  30%  of 
the  total  area  of  the  roughness  field  when  calculated  in  this  manner  — 
the  exact  range  was  from  0*3  to  35%.  The  use  of  the  non-dimensional 
variable  A  allows  the  inclusion  of  the  results  from  the  variation  in 
obstacle  height  (March  18,  1961)  with  the  rest  of  the  data. 

A  linear  least- square  error  method  was  used  to  obtain  a  power 
law  relating  the  computed  values  of  Zo ,  d ,  and  tq  to  a  non-dimen¬ 
sional  variable  A  over  the  limited  range  of  obstacle  densities  used 
in  this  study  (3.4  <  A  <  410).  It  will  be  discussed  later  that  these 
power  law  relationships  cannot  be  expected  to  hold  over  the  range  of 
all  possible  obstacle  densities,  but  they  permit  us  to  compare  the  re¬ 
sults  of  this  study  with  similar  experiments  conducted  in  ducts  and 
wind  tunnels. 

7.4.1  Aerodynamic  roughness  parameter 

The  variation,  in  the  aerodynamic  roughness  parameter  (zo )  with 
obstacle  density  from  the  results  of  the  present  study  is  indicated  in 
Fig.  9.  Notice  that  Zo  varied  by  more  than  two  orders  of  magnitude 
from  its  initial  value  over  the  ice.  The  line  of  best  fit  shows  that  Zo 
is  proportional  to  A“**^  A  decrease  in  Zo  is  indicated  for  the  final 
increase  in  obstacle  density;  however,  this  may  be  spurious  due  to 
scatter  in  the  data. 

Schlichting  (1937)  has  reported  on  roughness  experiments  in  the 
laboratory  using  obstacles  of  various  shapes  and  sizes  fastened  to 
the  bottom  of  a  water-filled  duct  under  fully  developed  turbulent  con¬ 
ditions.  He  has  described  the  hydrodynamic  roughness  of  these  ob¬ 
stacles  in  terms  of  kg  ,  the  diameter  of  tightly  packed  grains  of  sand 
that  would  produce  an  equivalent  hydrodynamic  roughness.  To  com¬ 
pare  Schlichting' s  data  with  our  own  we  consider  that  Zq  equals  h/N 
for  fully  turbulent  flow  in  rough  ducts  when  h  =  kg  .  Nikuradse  (1933) 
found  N  =  30  assuming  a  logarithmic  velocity  distribution  for  the 
entire  duct  cross  section,  and  k  =  0.40.  In  Section  8  Lettau  describes 
a  new  mathematical  model  for  the  universal  velocity  distribution  law 
under  fully  turbulent  conditions  in  rough  ducts  for  which  Nikuradse' s 
measurements  yield  k  =  0.428  and  N  =  23.9.  The  expression 
Zo  =  kg/23.9  is  used  to  transform  Schlichting '  s  values  of  ks  into 
values  of  Zo  .  This  is  done  in  Fig.  10,  where  the  ordinate  is  a 


Fig.  8.  The  range  of  obstacle  densities  (drawn  to  scale)  and 
the  general  shape  and  size  of  the  roughness  field  for  three 
controlled  wind  profile  experiments  on  the  ice  of  Lake  Mendota, 
1961. 


Case  I 

January  27 

S.  A  =  48.5, 

23.0, 

11.7 

Case  II 

February  15 

S.  A.  =  14.7, 

7.7, 

4.0 

Case  III 

February  24 

S.  A.  =  3.0, 

1.5, 

0.8 ,  0.4 
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non-dimensional  variable  defined  as  the  ratio  of  the  aerodynamic 
roughness  length  (Zq)  to  the  geometric  height  of  the  obstacle  (h) . 

The  similarity  between  the  line  of  best  fit  from  the  basket  experi¬ 
ments  and  the  data  from  Schlichting' s  experiments  is  apparent.  No¬ 
tice  that  Zo/h  for  Schlichting' s  spheres  reaches  a  maximum  at 
s...Arvalues.Qt4PPrQ^imately  3  and  decreases  for  further  increases  in 
cbstacle  density. 

More  recently,  Sayre  (1961)  has  desaibed  the  effect  of  rough¬ 
ness  spacing  on  open  channel  flow  for  roughness  elements  consisting 
i  of  sheet-metal  baffles.  He  notes  the  general  tendency  for  his  spac¬ 
ing  parameter  (hydrodynamic  roughness/obstacle  height)  to  increase 
with  inaeasing  roughness  density.  Using  our  notations,  his  obstacle 
density  varied  from  A  approximately  equal  to  67  to  A  approximately 
equal  to  4 . 

The  reason  for  the  observed  variations  in  aerodynamic  and  hydro- 
dynamic  roughness  can  be  explained  as  follows.  The  addition  of  ob- 
i  Stacies  at  an  initially  smooth  surface  causes  increased  resistance  to 
the  flow  near  the  surface.  If  the  obstacle  density  approaches  satura¬ 
tion  (100%  of  the  total  area)  and  if  there  is  little  separation  between 
"  obstacles  (i.  e. ,  parallel  sides)  another  smooth  surface  is  established 
at  the  height  of  the  obstacles.  Therefore,  the  optimum  roughness  may 
occur  at  an  obstacle  density  less  than  saturation. 

It  is  perhaps  instructive  to  compare  certain  expressions  that  have 
been  used  to  relate  Zq  to  h  with  the  plots  of  Zo/h  versus  A  illus¬ 
trated  in  Fig.  10.  From  observations  obtained  over  snow,  grassland, 
and  beet  fields,  Paeschke  (1937)  has  found  Zq  =  h/7.35  .  In  Section  3 
Kung  shows  that  Zg  =  h*'*’/ 17.4  for  tall  vegetation  ranging  from  grass 
types  (Zq  -  h/lO)  to  forest  trees  (Zg  =  h/4).  The  ratios  given  by 
Nikuradse  and  Lettau  have  been  discussed.  The  values  of  Zg/h  ob¬ 
tained  by  these  investigators  are  illustrated  along  the  left-hand  margin 
of  Fig.  10.  The  expressions  of  both  Paeschke  and  Kung  seem  to  be 
valid  only  for  the  maximum  obstacle  densities  used  in  this  study  (i.e.. 
A'=  3.4  to  7  or  roughly  30  to  14%  of  the  surface  area  covered  with 
baskets),  and  for  the  region  containing  the  peak  in  the  roughness 
curve  of  Schlichting' s  spheres.  The  ratios  of  both  Nikuradse  and 
Lettau,  which  pertain  to  sand  type  roughness  in  duct  flow,  seem  to 
be  valid  only  when  A  *  16  for  the  baskets  and  A  »  1.3  or  30  for 
the  spheres.  Clearly,  aerodynamic  roughness  is  a  function  of  ob¬ 
stacle  density  as  well  as  obstacle  height. 


Fig.  9.  The  aerodynamic  roughness  parameter  (Zo )  versus  A  for 
five  controlled  wind  profile  experiments  on  the  ice  of  Lake  Mendota, 
1961.  The  solid  line  indicates  the  linear  least-sguare  error  fit. 

The  broken  line  indicates  an  asymptotic  approach  to  the  average 
value  of  Zo  over  the  unmodified  ice. 
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7.4.2  Zero  point  displacement 

i 

The  values  of  d  obtained  over  the  unmodified  ice  are  very  small 
(approximately  2  cm)  and  are  presented  in  Sectioa  7.6.  Recalling  the 
boundary  condition  for  equation  (1)  and  noting  that  all  values  of  d 
obtained  from  wind  profiles  over  obstacles  are  negative  and  that 
|d|  <  h,  we  conclude  that  V  s  0  for  the  logarithmic  model  of  the 
wind  profile  at  some  height  between  the  surface  of  the  ice  and  the 
tops  of  the  obstacles.  Actual  wind  profiles  and  computed  solutions 
to  equation  (1)  from  the  actual  wind  profile  .data  are  plotted  on  a 
linear  height  scale  in  Fig.  11  to  Illustrate  the  zero  point  displace¬ 
ment.  The  cases  shown  refer  to  the  wind  profile  over  the  unmodified 
ice  and  the  wind  profile  over  the  maximum  obstacle  density  (S.  A.  = 
0.4  m* )  on  February  24. 

The  variation  in  d  with  obstacle  density  is  shown  in  Fig.  12. 

The  ordinate  is  a  non-dimensional  variable  defined  as  the  ratio  -d 
over  h .  The  line  of  best  fit  shows  that  -d/h  is  proportional  to 
A-0.29  increases  with  obstacle  density  but  at  a  signif¬ 

icantly  smaller  rate  as  Zo  increases  with  obstacle  density. 

Neither  of  these  power  law  relationships  should  be  expected  to 
hold  for  further  increases  in  obstacle  density.  The  quantity  -d/h 
must  approach  unity  asymptotically  as  obstacle  density  approaches 
100%  ,  while  ZqM  may  have  a  maximum  for  some  obstacle  density 
and  be  decreasing  as  the  obstacle  density  approaches  100%  (as  in 
the  case  of  Schlichting' s  spheres). 

7.4.3  Surface  stress 

The  variation  in  surface  stress  tq  with  obstacle  density  is  illus¬ 
trated  in  Fig.  13.  The  ordinate  is  the  square  of  a  non-dimensional 
variable  defined  as  the  ratio  os  the  friction  velocity  at  the  modified 
rough  surface  to  the  friction  velocity  at  the  ice  surface  —  (V*/Vfce)^. 
This  expression  reduces  the  ratio  tq  for  the  obstacle  field  over  to  for 
the  ice  if  k  is  a  true  constant.  Notice  that  to  over  the  maximum 
density  is  nearly  an  order  of  magnitude  greater  than  to  over  the  ice. 
The  line  of  best  fit  shows  that  is  proportional  to 

A  heavy  dashed  line  extends  from  the  line  of  best  fit  at  low  obstacle 
densities  to  indicate  that  (VVvJ'^g)^  must  approach  unity  asymptotic¬ 
ally  as  obstacle  density  becomes  small. 

Comparison  with  Schlichting' s  data  is  possible.  He  lists  values 
of  shearing  velocity  v*  for  the  rough  wall  of  his  duct  and  v*^  for  the 


Fig.  10.  The  ratio  of  the  aerodynamic  roughness  parameter  (zo) 
to  the  obstacle  height  (h)  versus  A.  The  solid  line  indicates 
the  line  of  best  fit  from  the  basket  study.  The  broken  lines  re¬ 
fer  to  recomputed  data  from  the  laboratory  e^xperiments  of 
Schlichting  (using  Zo  =  kg/23.9).  The  ratios  of  Zo/h  obtained 
by  various  investigators  are  indicated  along  the  left-hand  margin 
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smooth  wall*  for  a  number  of  different  obstacles  and  a  variety  of  duct 
flow  speeds.  Using  this  data*  values  of'(v‘^f/v*g)^  have  been  com¬ 
puted  for  spheres  and  spherical  segments  as  functions  of  A  and  are 
Illustrated  in  Fig.  13.  The  values  pertain  to  flow  speeds  at  the  cen¬ 
ter  of  the  channel  of  sio  ±Z0  cm/sec.  The  rate  of  increase  Ih  sur¬ 
face  stress  with  obstacle  density  indicated  by  Schlichting'  s  data  is 
less  than  that  indicated  by  the  baskets.  A  maximum  surface  stress 
is  indicated  for  the  spheres  at  A  *  3 . 

Moore  (1951)*  in  a  series  of  wind  tunnel  tests*  determined  aero¬ 
dynamic  drag  as  a  function  of  the  density  of  small  wooden  cubes  ar¬ 
ranged  in  a  diagonal  square  pattern.  He  states  that  a  curve  of  drag 
versus  spacing  exhibited  a  broad  maximum  extending  from  a  relative 
spacing  of  3  to  6  times  the  side  of  the  cube*  which  corresponds  to 
16  <  A  <  49  in  our  notation.  This  effect  was  not  observed  in  the 
basket  experiments.  ‘  A  slight  decrease  in  does  occur  for 

the  final  Increase  in  obstacle  density:  however*  this  may  be  spurious 
due  to  scatter  in  the  data. 

7.4.4  Geostrophic  drag  coefficient  —  Surface  Rossby  Number 

Attention  is  called  to  the  drastic  changes  in  the  appearance  of 
the  wind  profiles  and  in  the  values  of  the  computed  profile  parameters 
that  are  produced  by  the  controlled  modification  of  the  surface  rough¬ 
ness.  This  is  particularly  noticeable  in  the  profiles  of  February  24, 
1961  (see  Figs.  5  and  12).  Wind  shear  between  50  and  150  cm  is 
tripled*  wind  speed  immediately  above  the  obstacles  is  more  than 
halved,  Zq  is  inaeased  by  more  than  two  orders  of  magnitude  and  to 
by  nearly  one.  Ellison  (1957)  has  made  the  following  statement: 

"The  only  effect  of  a  change  in  the  value  of  zq  is  to  superimpose  a 
uniform  translation  on  the  whole  flow  without  modifying  its  internal 
mechanism."  Clearly*  our  experiments  show  that  a  change  in  z©  does 
have  more  complex  consequences  than  a  simple  translation.  There  is 
an  additional  adjustment  between  surface  stress  and  surface  wind 
speed.  This  leads  us  to  search  for  a  truly  ambient  wind  speed,  or 
reference  speed. 

Lettau  (1958)  has  suggested  that  the  horizontal  pressure  gradient, 
or,  consequently,  the  geostrophic  wind  in  the  atmospheric  surface 
layer,  assumes  the  role  of  an  over-all  forcing  function  similar  to  that 
of  the  pressure  head,  or  the  center  speed  in  duct  flow,  or  the  ambient 
speed  in  wind  tunnel  experiments.  He  has  shown  the  dependency  of 
the  geostrophic  drag  coefficient, 


Fig.  11.  Actual  wind  profiles  (solid  lines)  and  computed  solutions 
to  equation  1  (broken  lines)  over  the  unmodified  ice  (I )  and  over  the 
maximum  obstacle  density  ( II  —  S.  A  =  0.4  mM  for  a  controlled  wind 
profile  experiment  on  the  ice  of  Lake  Mendota,  February  24,  1961. 
The  zero  point  displacement  (d)  for  each  of  the  computed  profiles  is 
indicated  along  the  left-hand  margin. 
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C  =  {To/p)Vv„  ^  (2) 

g*o 

on  what  he  defines  as  the  Surface  Rossby  Number, 

SSo  =  V„  /z,  f  .  (3) 

Lettau  obtained  an  empirical  relationship  relating  C  and  Roo 
from  mlcrometeorologlcal  data  available  at  that  time.  More  recently, 
Lettau  (see  Section  9)  has  computed  the  geostrdphic  drag  coefficient 
of  theoretical  wind  spirals  (a  generalized  concept  of  Ekman  spirals) 
as  functions  of  goo.  The  present  theoretical  and  1958  empirical  rela¬ 
tionships  along  with  the  data  from  the  basket  experiments  are  pre¬ 
sented  in  Fig.  14. 

The  shape  of  the  experimental  curve  is  steeper  than  that  of  the 
1958  empirical  relationship,  which  in  turn  is  steeper  than  the  theo¬ 
retical  curve.  This  result,  however,  should  be  expected.  In  the  the¬ 
ory,  full  adjustment  between  wind  and  stress  distribution  throughout 
the  entire  planetary  boundary  layer  is  assumed.  In  view  of  the  average 
thickness  of  this  layer,  upwind  fetches  over  uniform  terrain  of  the 
order  of  100  km  are  required.  In  the  selection  of  conventional  micro- 
meteorological  sites,  actual  fetches  over  uniform  terrain  will  normally 
be  between  100  m  and  10  km.  The  roughness  fetches  in  our  experi¬ 
ments  were  always  less  than  100  m.  Clearly,  such  short  distances 
will  not  allow  the  wind  profile  to  establish  a  balance.  The  establish¬ 
ment  of  a  balance  requires  a  change  in  the  direction  of  the  surface 
wind  —  namely,  an  increasing  component  towards  lower  pressure  to 
compensate  for  increasing  frictional  drag.  No  changes  in  wind  direc¬ 
tion  were  evident  in  our  experiments. 

The  discrepancy  between  the  January  27  data  and  the  rest  of  the 
points  in  Fig.  14  is  probably  due  to  an  over-estimate  of  Vg  q  .  The 
method  used  in  the  determination  of  Vg^o  values  obtained  are 

summarized  in  the  Appendix  (Table  5). 

7.  5  Supplementary  Results 

A  determination  of  tq  independent  of  equation  (1)  is  desirable. 
Lettau  (1961)  has  outlined  such  a  method  based  on  continuity  of  hori¬ 
zontal  momentum  for  two-dimensional,  steady  state  flow.  The  gener¬ 
alized  equation  of  conservation  reduces  to  the  form, 

^  (p^iu)  (P‘3W)  =  0  (4) 

where  q  =  conservative  property  of  the  fluid  per  unit  mass. 
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Integrating  with  respect  to  height  from  the  surface  to  some 
height  z ,  and  interchanging  the  sequence  of  integration  (with  re¬ 
spect  to  z )  and  differentiation  (with  respect  to  x ),  we  obtain 

8 

=  (pqw)^  P**"***  • 

Setting  q  =  1  in  (5),  and  noting  that  for  a  horizontal  solid 
boundary  (pw)o  =  0 ,  we  have 

8 

=  ■  to  Jo 

Equation  (6)  can  be  used  to  determine  the  vertical  velocity  at 
various  heights  (z)  in  a  truly  two-dimensional  flow. 

Substituting  (6)  into  (5)  gives 

z  z 

‘^'>’*>0=  -^zfxL  ^  to 

Select  a  level  H  such  that  for  z  >  H  we  have  8q/8x  =  0  for 
any  x.  Then, 

a 

to  Jo 

The  result  of  averaging  over  the  finite  distance  R  =  X2  -  Xi  is 
Xj _ X2  j  H 

(pqw)o  =  R  4  p[uj(q2  -  Qjj)  -  ui(qi  -  qjjlJdz .  (9) 


Finallyi  for  the  case  of  horizontal  momentum,  substitute  q  =  u  , 

q  =  U  (such  that  dU/ax  =  0). 
rl 

Xi _ Xi  X, _ Xj  p  H 

To  =  -(puw)^  ■  R  4  [U2(U‘U2)'“l(U-Ul)]dZ 

P 

=  -  [U(Ui -u,)-(ui -u,  )]dz.(10) 


This  solution  reduces  to  the  vonKdrmdn  relation  (see  Sutton,  1953, 
p.  55)  for  the  special  case  when  Ui  =  U  at  Xi  =  0  and  U2  =  u  at 
X2  =  X,  which  is  true  for  the  development  of  a  flat-plate  boundary 
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layer  downstream  from  the  leading  edge,  in  an  initially  uniform  wind 
tunnel  flow. 

7.5. 1  Experimental  determination  of  the  Kdrmdn  constant 

Combining  the  result  of  equation  (10)  with  the  slope  of  the  wind 
profile  V*  from  equation  (1)  we  obtain  a  direct  estimate  of  the  K^rmdn 
constant , 

Xl _ X2  I 

K  =  (  TO  /p)*/V*  ,  (11) 

Equation  (11)  is  used  to  estimate  k  during  three  experiments  in 
which  less  than  db3%  change  in  wind  speed  occurred  at  the  upper¬ 
most  anemometer  level  (see  Figs.  5,  6,  and  7).  The  method  is  not 
exact  because  the  velocity  distribution  above  the  highest  anemometer 
level  is  unknown.  The  K-values  obtained  are  listed  in  Table  1. 

The  average  of  the  nine  K-estimates  is  0.41  with  a  standard  devi¬ 
ation  of  0.07  .  The  average  value  is  in  satisfactory  agreement  with 
accepted  values  and  the  spread  is  not  unreasonable  for  a  calculation 
of  this  type. 

The  numerical  evaluation  of  equation  (10)  was  accomplished  by 
plotting  the  profiles  on  a  linear  height  scale  and  reading  values  of 
U2  (wind  profile  over  baskets)  and  Uj  (wind  profile  over  unmodified 
ice)  at  5  cm  height  intervals.  The  velocity  defect  (u2  -  Ui )  was  ob¬ 
tained  directly  from  a  velocity  defect  profile,  thus  eliminating  the 
effect  of  systematic  errors  in  measurement. 

7.5.2  Vertical  velocity  pattern  over  the  roughness  field 

As  previously  mentioned,  equation  (6)  allows  the  deduction  of 
the  vertical  motion  pattern  over  the  roughness  field.  The  experiment 
of  February  19  is  well  suited  for  this  purpose  because  the  roughness 
fetch  R  was  varied.  The  resulting  values  of  mean  vertical  velocity 
apply  to  mean  distances  of  10,  20,  and  24  meters  from  the  leading 
edge  of  the  roughness  field  and  indicate  the  rate  at  which  the  mean 
streamlines  are  displaced  upward. 

Assuming  w  ^  ,  the  value  of  n  was  obtained  for  each  of 

three  levels  from  a  log-log  plot  of  the  data  (see  Fig.  15).  The  ex¬ 
ponent  n  varied  from  -0.7  to  -0.6  to  -0.5  for  w  evaluated  at  1,  2, 
and  3.4  meters,  respectively,  above  the  surface.  Vertical  motions 
of  3  to  4  cm/sec  are  indicated  at  heights  of  2  m  and  distances  of  15 
to  20  m  downwind  from  the  leading  edge  of  the  roughness  field. 


Fig.  14.  Geostrophic  drag  coefficient  versus  Surface  Rossby 
Number.  The  short  broken  lines  refer  to  the  values  obtained 
during’  five  controlled  wind  profile  experiments  on  the  ice  of 
Lake  Mendota,  1961. 
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If  we  assume  that  at  and  above  some  level  the  mean  flow  pattern 
is  undisturbed  by  the  presence  of  the  roughness  field,  momentum  con¬ 
tinuity  for  two-dimensionali  steady  state  flow  requires  an  Increase 
in  horizontal  wind  speed  below  this  level  to  compensate  for  the  de¬ 
crease  in  horizontal  wind  speed  within  and  immediately  above  the 
roughness  field.  This  effect  is  observed  in  the  profiles  of  February  15 
(see  Fig.  3). 

7.5.3  Wind  speed  between  the  obstacles 

The  wind  speed  in  the  air  layer  between  the  roughness  elements 
is  strongly  dependent  upon  the  position  of  the  nearest  obstacles.  The 
baskets  near  the  anemometer  mast  were  arranged  such  that  the  mast 
occupied  the  position  that  would  ordinarily  be  occupied  by  an  obstacle. 
As  a  result,  the  indicated  speeds  should  be  quite  representative. 

Fig.  16  illustrates  the  ratio  of  wind  speed  at  20  cm  within  the  layer 
of  air  between  the  baskets  (V20 )  to  the  wind  speed  at  the  same  level 
over  the  unmodified  surface  (Vzoflce)  as  a  function  of  A.  Notice 
that  V2o/V2o,lce  ^s  0.5  when  only  15%  of  the  roughness  field  is 
covered  with  obstacles. 

7.5.4  The  drag  coefficient  of  the  individual  basket 

The  drag  coefficient  of  an  individual  bushel  basket  was  not  deter¬ 
mined  directly.  Schlichting  (1937)  gives  an  expression  for  determining 
the  drag  coefficient  (Cf)  of  an  individual  roughness  element  from 
measurements  of  the  geometry  of  the  obstacle  field,  the  fluid  velocity 
at  the  obstacle  height,  and  the  friction  velocity  at  the  rough  wall 
(basket  roughness  field)  and  the  smooth  wall  (unmodified  ice). 

A  nearly  constant  value  of  Cf  =  0.45  was  obtained  for  the  bushel 
basket  using  Schlichting' s  equation.  This  was  in  contrast  to  the  be¬ 
havior  of  Schlichting' s  spheres  which  show  a  decrease  in  Cf  from 
approximately  0.89  at  A  =  100  to  approximately  0.20  at  A  =3. 

7.6  Error  Analysis 

7.6.1  Anemometry 

Simultaneous  comparison  of  the  twelve  Thornthwaite  anemometers 
on  a  horizontal  bar  was  made  immediately  before  and  after  each  ex¬ 
periment.  The  results  were  averaged  and  correction  terms  applied  to 
individual  anemometers.  In  general,  corrections  were  less  than  1%. 
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In  an  attempt  to  obtain  a  more  detailed  profile,  it  was  necessary 
to  stagger  the  anemometers  on  the  mast.  Slight  changes  in  wind  di¬ 
rection  caused  the  anemometers  more  nearly  "upwind"  of  the  mast  to 
slow  down.  The  resulting  spread  in  measured  wind  speeds  at  adja¬ 
cent  levels  is  obvious  oh  the  profiles  of  February  24  and  February  19 
(see  Figs.  5  and  6).  This  spread  introduces  some  uncertainty  into 
the  wind  profile,  but  the  use  of  the  least-square  technique  in  the 
analysis  of  the  wind  data  avoids  any  serious  m^interpretations. 

7.6.2  Effect  of  thermal  stratification 

The  Richardson  number,  Ri,  was  computed  at  four  levels  from 
averages  of  temperature  samples  obtained  every  30  seconds  and  from 
the  10-minute  average  wind  shear.  These  values  are  listed  in  Table  3 
in  the  Appendix.  A  practical  criterion  for  thermal  stratifications  suf¬ 
ficiently  close  to  neutral  is  [Rimn  |  <  0.003 . 

Lettau  (1957)  defines  a  bulk  parameter  (Ri)'  as  a  measure  of  the 
over-all  convective  stability  of  a  layer. 

(Ri)'  =  2lRi(2,)/Z  *,  •  <12) 

i  1 

For  practical  purposes  (Ri)'  expressed  in  units  of  10“^ per  meter  can 
be  considered  numerically  equal  to  Ri  times  10"^  at  z  =  100  cm. 

For  the  Great  Plains  Turbulence  Field  Program,  using  the  bulk  param¬ 
eter  (Ri)'  Lettau  (1957,  p.  331)  has  classified  as  weak  lapse  (Ri)' 
between  -9  and  -5,  as  neutral  (Ri)'  between  -4  and  +4,  and  as 
weak  inversion  (Ri)'  between  +5  and  +6.  These  terms  are  used  in 
a  relative  sense  because  the  over-all  wind  speed  during  the  program 
was  large. 

The  values  of  Rin3  and  (Ri)'  from  this  study  have  been  averaged 
for  all  wind  profiles  over  the  baskets  and  for  the  wind  profiles  over 
the  ice.  The  results  are  summarized  in  Table  2.  The  presence  of  the 
baskets  increases  the  wind  shear  sufficiently  to  reduce  Rins  and 
(Ri)'  to  near  neutral  conditions.  The  extent  of  reproducibility  of  the 
wind  profile  parameters  over  the  ice  is  also  noted  in  Table  2. 

Hourly  observations  of  cloud  cover,  temperature,  wind  direction 
and  wind  speed,  and  station  pressure  have  been  obtained  from  the 
United  States  Weather  Bureau  located  at  Truax  Field,  Madison,  Wis-  - 
consin.  These  observations,  along  with  computed  values  of  air  density 
are  summarized  for  each  experiment  in  Table  4  in  the  Appendix. 


Fig.  16.  The  ratio  of  the  wind  speed  at  20  cm  within  the  layer  of 
air  between  the  obstacles  (V20)  to  the  wind  speed  at  20  cm  over  the 
unmodified  ice  (V2o,ice)  versus  A  for  five  controlled  wind  profile 
experiments  on  the  ice  of  Lake  Mendota,  1961.  The  broken  line  in¬ 
dicates  the  limiting  values  of  (V2o/V2o,ice)  • 
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7.6.3  Representativeness  of  the  wind  profiles 

There  is  some  question  concerning  which  anemometer  levels  to 
include  in  the  wind  profile  parameter  computations. 

Anemometers  below  and  Immediately  above  the  basket  height  re-' 
fleet  the  flow  over  and  around  an  individual  obstacle  and  not  the 
statistical  effect  of  all  obstacles.  A  number  of  computations  were 
made  on  all  profiles  omitting  the  lowest  anemometers  one  by  one. 

It  was  found  that  systematic  variations  in  the  profile  parameters 
could  be  obtained  only  by  excluding  the  anemometer  located  immedi¬ 
ately  above  the  baskets  (31,  32  cm),  and  all  levels  within  the  layer 
of  air  between  the  baskets.  This  result  was  considered  for  all  com¬ 
putations. 

The  possibility  that  in  some  cases  the  anemometers  at  the  high¬ 
est  levels  might  be  outside  the  modified  layer  was  di.scussed  in  Sec¬ 
tion  7.4.  The  problem  is  somewhat  akin  to  the  growth  of  an  "internal 
boundary  layer"  referred  to  in  fluid  dynamics.  An  attempt  was  made 
on  February  19  to  Investigate  this  effect  by  systematic  variation  of 
the  upwind  extent  of  the  roughness  field  while  keeping  the  density 
distribution  constant.  The  resulting  profiles  are  indicated  in  Fig.  6 
and  are  discussed  below. 

R  =  20  meters:  A  "kink"  appeared  in  the  profile  above  the  160  cm 
level.  This  suggests  that  higher  levels  were  not  yet  modified’ by  the 
roughness  field.  A  similar  "kink"  was  observed  in  the  profiles  of 
February  24  when  R  =  18  m  (see  Fig.  5). 

R  =  40  meters:  Further  modification  (approximately  a  5%  decrease 
in  wind  speed)  occurred  in  the  middle  portion  of  the  wind  profile.  The 
"kink"  in  the  profile  was  not  as  sharp  as  on  the  previous  run  and  the 
slope  of  the  profile  above  the  160  cm  level  increased.  This  indicates 
that  modification  extended  above  the  160  cm  level. 

R  =  48  meters;  Further  modification  (Ipss  than  a  5%  decrease  in 
wind  speed)  occurred  in  the  middle  and  upper  portion  of  the  profile. 
Unfortunately,  no  second  run  was  made  for  this  fetch. 

If  all  anemometer  levels  above  the  height  of  the  baskets  are  used 
in  the  computation  of  the  profile  parameters  for  the  case  R.=  20  m , 
we  find  -d  =  44.6  cm  and  Zo  =  0.025  cm  .  If  anemometer  levels  above 
the  180  cm  level  (i.  e. ,  above  the  "kink")  are  omitted,  we  find 
-d  =  28.4  cm  and  Zo  =  0.322  cm.  The  latter  values  are  in  close  agree-,., 
ment  with  the  results  of  the  computations  for  R  =  40  m  in  which  all 
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levels  were  Included.  Using  anemometer  levels  located  above  the 
modified  layer  will  lead  to  overestimating  the  -d  value  and  under¬ 
estimating  zo  .  The  decrease  in  -d  and  the  Inaease  in  Zo  for 
R  s  48  m  must  be  considered  uncertain  because  only  one  run  was 
made  for  this  roughness  fetch. 

The  following  convention  was  adopted: 

1.  Anemometer  levels  above  180  cm  were  not  included  in  the 

analysis  of  the  profile  parameters  for  R  <  20  meters. 

2.  All  anemometer  levels  were  included  in  the  analysis  of 

the  profile  parameters  for  R  >  20  meters. 

As  previously  mentioned,  Lettau  (1959)  has  suggested  as  a  "rule  of 
thumb"  that  the  highest  instrument  level  should  not  exceed  «  l/SOth 
of  the  upwind  distance  from  significant  discontinuities  in  the  surface 
structure.  On  the  basis  of  the  appearance  of  the  profiles,  instrument 
levels  at  heights  =  l/lOth  the  upwind  fetch  over  the  baskets  have 
been  included.  As  a  result,  the  parameters  obtained  may  not  be  com¬ 
pletely  representative  of  the  surface. 

The  problem  of  the  growth  of  the  internal  boundary  layer  must  be 
examined  more  closely  in  future  experiments. 

7.7  Conclusions 

/ 

The  value  of  controlled  modification  of  surface  roughness  on  wind 
profile  structure  has  been  demonstrated  by  the  results  of  this  study. 

The  systematic  increase  in  obstacle  density  on  the  ice  surface 
results  in  a  tripling  of  the  wind  shear,  an  increase  of  more  than  two 
orders  of  magnitude  in  zo ,  and  an  increase  of  nearly  one  order  of 
magnitude  in  tq  .  In  addition  to  providing  direct  information  on  the 
variation  of  the  wind  profile  parameters  with  the  geometry  of  the  sur¬ 
face,  the  technique  has  obvious  applications  to  the  study  of  various 
other  micrometeorological  problems. 

Future  plans,  as  outlined  by  Lettau  (1959),  include  the  use  of 
obstacles  with  different  albedo  to  absorb  or  reflect  solar  radiation  — 
thus  producing  varying  correlation  between  the  momentum  sink  struc¬ 
ture  and  the  heat  source  pattern.  Also,  the  vertical  fluxes  of  heat 
and  moisture  and  the  power  spectrum  of  turbulence  near  the  surface 
can  be  observed  as  functions  of  surface  roughness,  wind  shear,  and 
thermal  stratification. 
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Table  3.  Roughness  field,  wind  profile,  and  temperature  profile  parameters  for  20  ten-minute  mean  runs 
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All  listed  values  are  ±5% 

AR  =  distance  between  basket  rows 

AS  (not  listed)  =  distance  between  baskets  in  each  basket  row  =  S.  A./AR  . 


Table  3,  continued — Temperature  profile  parameters 
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Table  3,  continued  — Wind  profile  parameters 
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Table  4.  Hourly  observations  of  sky  cover,  station  pressure,  air  temperature,  wind  direction,  and  wind 
speed,  and  computed  air  densities  for  five  controlled  wind  profile  experiments  on  the  ice  of  Lake  Mendota 
1961. 
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A  Generalized  Mathematical  Model  of 


the  Mean-Velocity  Distribution 
in  Fully  Turbulent  Duct  Flow 


Heinz  H.  Lettau 

Department  of  Meteorology 
and 

Department  of  Civil  Engineering 
University  of  Wisconsin 


Abstract.  Existing  mathematical  models  of  the  mean- 
velocity  distribution  in  fully  turbulent  duct  flow  are 
critically  reviewed,  in  terms  of  analytical  expressions 
for  the  length-scale  of  turbulence  (as  a  function  of 
radial  distance  from  the  center  of  the  duct).  A  new 
theoretical  form  for  the  length-scale  of  turbulence  is 
proposed  which  contains  no  empirical  constant  in  addi¬ 
tion  to  K^rm^n' s  constant  and  which  satisfies  a  principle  ^ 
of  similarity  between  duct  flow  and  atmospheric  boundary- 
layer  flow.  The  resulting  universal  velocity-distribution 
in  cylindrical  ducts,  with  smooth  or  rough  walls,  agrees 
satisfactorily  with  the  observations  reported  by  Nikuradse 
and  Laufer.  It  is  found  that  Kirm^n' s  constant  is  a  true 
constant  (equal  to  0.428)  only  for  fully  developed  flow 
in  ducts  with  rough  walls.  It  is  suggested  that  turbulent 
viscosity  is  adjusted  so  that  an  "effective  Reynolds" 
number  is  less  than  a  critical  value  (in  the  vicinity  of 
1,000),  independent  of  the  actual  Reynolds  number. 

8. 1  List  of  Symbols 

8.1. 1  Fluid  properties,  and  duct  characteristics 
'  p  =  fluid  density 

V  =  kinematic  viscosity  of  the  fluid 
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2R  s  lateral  diameter  of  the  duct 
k  =  average  height  of  the  wall -roughness  elements 
Zq  =  aerodynamic  roughness -length  of  the  wall 

8.1.2  Mathematical  symbols,  subscripts 

'  =  prime,  denoting  differentiation  with  respect  to  y  (or  -x ) 

=  bar,  denoting  an  average  across  the  duct 
0  =  subscript,  indicates  the  value  of  a  quantity  at  the  wall 
c  =  subscript,  indicates  the  value  of  a  quantity  at  the  center 

8.1.3  Dimensional  variables  and  flow  parameters 

r  =  lateral  distance  from  the  center  (independent  coordinate) 
u  =  u(r)  =  mean-velocity  profile 
U  =  =  maximum  velocity  at  the  center 

T  =  T(r)  =  shearing-stress  profile 

I 

u***  =  (t/p)^  =  u’*‘(r)  =  friction-velocity  profile 

}jc 

Uo  =  wall  friction-velocity 

L  =  u’V(9u/8r)  =  L(r)  =  profile  of  characteristic  length-scale  of 
turbulence 

8.1.4  Dimensionless  duct  characteristics,  variables,  and  flow 
parameters 

Re  =  2Ru/v  =  Reynolds  number 
y  =  r/R  =  dimensionless  lateral  distance  from  the  center 
X  =  1  -  y  =  dimensionless  lateral  distance  from  the  wall 
Xq  =  Zq/R  =  aerodynamic  roughness-parameter  of  the  wall  surface 
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N  s  k/zo  -  Nikuradse's  roughness  ratio 
t  =  correction-term  function 

T)  =  (U-u)/uo  =  nly)  =  dimensionless  velocity-defect  profile 
n  -  (U-u)/uo  =  average  dimensionless  velocity-defect 
C  =  uJ/U  =  I/ho  =  wall  friction-velocity  coefficient 
K  =  K^rm4n' s  numerical  constant 

t 

F  =  L/R  =  F(y)  =  dimensionless  length-scale  of  turbulence 
f  =  4to/(puV2)  =  conventional  friction-factor 
a,  c.  m,  n  =  numerical  constants 
8.Z  Introduction 


The  discussion  will  be  restricted  to  fully  turbulent  steady  mean 
flow  in  a  straight  duct  of  uniform  cross  section,  with  a  defined  wall 
roughness.  Such  conditions  are  realized  for  Reynolds  numbers  larger 
than  approximately  10^,  and  distances  of  at  least  40  to  50  diameters 
from  the  entry- of  the  duct.  For  this  flow,  the  lateral  distribution  of 
mean  velocity  will  be  independent  of  the  down-stream  coordinate  (i.e., 
uniformity  of  flow),  and  also  of  the  fluid  viscosity.  For  steady  mean- 
flow  states,  the  Navier-Stokes  equation  reduces  to  an  equality  be¬ 
tween  the  lateral  gradient  of  eddy  shearing  stress  and  the  down-stream 
pressure  gradient  which  equals  a  constant  for  a  given  experiment.  Con¬ 
sequently,  the  shearing-stress  profile  must  be  exactly  linear  in  lateral 
distance  r,  with  =  0  at  the  center;  thus 

T  =  Toy  ;  or,  u*  =  u^NTy  .  (1) 

where  tq/R  is  determined  by  the  constant  pressure-gradient.  Fluid 
compressibility  effects  will  be  neglected. 

The  velocity  profile  must  satisfy  the  boundary  conditions  that,  at 
the  center  of  the  duct,  u^  equals  its  maximum  value  U ,  while  Uq  =  0 
at  the  wall.  In  place  of  u(r)  we  can  consider  the  velocity  defect  U  -  u  , 
or  its  dimensionless  form  Ti(y)  with  the  boundary  conditions  of  r\^-  0 
at  the  center,  and  =  U/uf  =  maximum  value  at  the  wall. 

Let  a  characteristic  length-scale  (L)  of  turbulence  be  defined  with 
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the  aid  of  the  dimensionless  equation 

L/R  =  F(y)  =  uVu’  (2) 

Equations  (1)  and  (2)  determine  uniquely  the  mean-velocity  profile 
u(r) ,  or  ii(y)  in  the  following  relation 

»1'  =  F‘'')y  :  or.  ii  =  F-'y*dy  .  (3) 

The  integral  in  equation  (3)  can  be  evaluated  when  F(y)  is  a  known 
analytical  function  of  y .  It  can  also  be  stated  that  a  defined  func¬ 
tion  F(y)  must  exist  if  a  mean-velocity  profile  is  determined. 

Note  that  only  dimensional  and  purely  kinematic  arguments  are 
necessary  for  the  definition  of  L  in  equation  (2).  This  length  L  is 
basically  the  same  as  the  "mixing  length"  which  was  introduced  by 
Prandtl  (1925)  and  expressed  independently  by  Taylor.  In  the  first 
part  of  his  famous  discussion  of  the  "statistical  theory  of  turbulence," 
Taylor  (1935)  made  the  following  remarks:  "This  length  could  only  be 
defined  in  relation  to  the  definite  but  quite  erroneous  conception  that 
lumps  of  fluid  behave  like  the  molecules  of  a  gas,  preserving  their 
identity  till  some  definite  point  in  their  path,  when  they  mix  with 
their  surroundings  and  attain  the  same  velocity  and  other  properties 
as  the  mean  value  of  the  corresponding  property  in  the  neighborhood." 
Although  Taylor  referred  specifically  to  a  Lagrangian  discussion,  his 
statements  seem  to  have  resulted  in  a  rather  widespread  discredit 
and  a  rejection  of  the  basically  Eulerian  concept  as  used  in  shearing 
flow.  In  fact,  since  equation  (2)  is  purely  kinematic,  the  concept  of 
mixing  need  not  enter.  Objections  can  only  be  raised  against  over¬ 
simplified  interpretation  but  not  against  the  legitimacy  of  the  kinematic 
definition  of  L.  As  an  Eulerian  concept  the  length  L  is  just  as  real¬ 
istic  and  useful  as  that  of  mean-velocity  for  shear-flow  in  ducts. 

In  fact,  the  historical  development  of  the  theory  of  fuHy  furbulent 
duct  flow  can' be  adequately,  documented  by  the  special  explicit  func¬ 
tions  F(y)  used  by  various  investigators.  A  review  of  the  problem,  in 
terms  of  F-functions,  also  establishes  the  basis  for  the  generalization 
of  the  mathematical  approach  to  the  mean-velocity  distribution  pro¬ 
posed  in  this  paper.  The  fact  is  important  that  atmospheric  boundary- 
layer  flow  is  more  closely  related  to  duct  flow  than  to  flat-plate 
boundary  layer  flow,  even  though  the  atmospheric  fluid  is  bounded 
only  on  the  lower  side.  The  physical  reason  for  this  statement  is  that 
for  constant  and  uniform  pressure  gradient  in  the  atmosphere  the  re¬ 
sulting  large-scale  air  flow  is  steady  and  horizontally  uniform,  be¬ 
cause  the  Coriolis  force  balances  exactly  the  forces  due  to  pressure 
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and  friction.  This  leads  to  two-dimensionality  of  air  flow  (i.  e. ,  sig¬ 
nificant  change  of  both  speed  and  wind  direction  with  height)  in  the 
layer  of  frictional  influence,  but  the  two  horizontal  components  of  air 
motion  are  horizontally  uniform,  as  in  duct  flow,  and  the  height  of 
the  friction  layer  is  independent  of  the  down-stream  coordinate.  In 
experimental  flat-plate  boundary-layer  flow,  with  or  without  pressure 
gradient,  the  Ck)riolls  force  is  negligible  (with  the  exception  of  ex¬ 
periments  involving  rapid  rotation  of  the  boundary),  and  the  friction 
force  is  essentially  balanced  by  the  inertia  force;  this  leads  to  a 
down-stream  growth  of  boundary-layer  thiclcness,  and  horizontally 
non-uniform  velocity,  which  makes  this  flow  type  significantly  dis¬ 
similar  to  atmospheric  flow. 

Consequently,  atmospheric  flow  resembles  duct  flow,  although 
only  the  lower  boundary  is  a  material  wall.  The  upper  wall,  or  the 
duct,  is  provided  by  the  geostrophic  control  (i.  e. ,  the  balance  be¬ 
tween  pressure  gradient  and  Coriolis  forces)  in  the  free  atmosphere, 
where  the  geostrophic  wind  blows  perpendicular  to  the  pressure 
gradient.  The  equivalent  of  the  one -dimensional  velocity-defect  in 
duct-flow  is  the  two-dimensional  "geostrophic  departure"  in  atmos¬ 
pheric  boundary-layer  flow.  With  this  in  mind,  some  degree  of  simi¬ 
larity  between  F-functions  for  experimental  duct  flow  and  atmospheric 
boundary-layer  flow  must  be  expected.  In  order  to  apply  duct-flow 
theories  to  the  atmosphere  (see  Section  9),  it  is  convenient  to  discuss 
F  as  a  function  of  distance  from  the  boundary,  i.  e. ,  F(x)  rather  than 
F(y) ,  considering  that  x  =  1  -  y .  A  substitute  for  the  duct-diameter 
is  the  height  of  the  "geostrophic  wind  level"  or  any  related  height 
parameter,  for  example,  the  "geostrophic  displacement  thickness. " 

The  practical  value  or  the  success  of  a  theoretic  model  can  only 
be  judged  on  the  basis  of  a  comparison  between  theoretical  and  well- 
established  empirical  results.  Thus  far,  the  geometrical  form  of  the 
duct  was  not  specified.  In  view  of  the  fact  that  the  largest  amount  of 
empirical  data  available  in  the  literature  concerns  turbulent  flow  in 
round  tubes,  let  us  restrict  the  discussion  to  cylindrical  ducts  of 
diameter  2R.  Readily  observed  quantities  are  the  wall  shear-stress, 
or  u*  (from  measurements  of  the  pressure  head,  pipe  length  and 
radius),  the  average  velocity  u  (from  discharge  measurements),  and 
the  mean-velocity  distribution  u(r),  including  uc  =  U  (by  direct  meas¬ 
urement  of  dynamic  pressure  at  various  positions  r/R).  Most  con¬ 
clusive  are  comparisons  of  the  dimensionless  characteristics  such 
as  Ti(y)  in  place  of  u(r),  including  the  "bulk"  characteristics  r)o 
and  fj  ,  i.  e. ,  the  extreme  and  the  average  velocity-defect.  The 
difference  (no  “  n)  relates  the  average  velocity  u  to  the  wall  shear- 
stress  and,  consequently,  results  in  a  theoretical  expression  for  the 
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conventional  friction-factor  —  reference  can  be  made  to  any  textbook 
of  fluid  mechanics;  see»  for  example*  Rouse  and  Howe  (1953, 
page  140)  — 


4to/(puV2)  *  f  s  8/(t1o  -  r\)^  .  (4) 

A  simpler  measure  of  friction  effects  is  the  wall  friction-velocity  co¬ 
efficient  C  =  lAio  •  By  definition, 

To=C*pU*  .  (5) 

In  equation  (5)  all  conventional  but  unnecessary  numerical  factors  are 
avoided.  This  form  is  especially  useful  in  atmospheric  flow  because 
the  velocity  U  is  determined  by  the  geos^’rophlc  wind  which  can  be 
computed  with  the  aid  of  the  horizontal  pressure  gradient. 

In  fully  turbulent  flow  the  friction-factors  f  and  C  are  independ¬ 
ent  of  the  Reynolds  number  because  the  average  and  the  extreme  mean- 
velocity  defects  are  determined  only  by  the  wall  roughness.  It  has 
been  well  established  for  cylindrical  ducts  that  f  depends  solely  on 
the  ratio  k/R  or  Zq  =  k/N  or  Xq  =  k/NR  ,  for  fully  turbulent  states. 
Nikuradse  (1933)  has  found  N  «  30  for  sandgrain  roughness.  Similar 
results  hold  true  for  the  atmosphere  and  the  C  coefficient.  Lettau 
(1959)  has  shown  that  C  is  a  function  of  a  dimensionless  combination 
(of  geostrobhic  speed.  Coriolis  parameter,  and  Zq  )  which  is  referred 
to  as  the  "Surface  Rossby  Number." 

With  the  exception  of  extremely  rough  walls  in  relatively  narrow 
ducts,  the  order  of  magnitude  of  k/R  is  l/lOO  or  less.  Consequently, 
Xq  *  k/30R  is  very  small  in  comparison  to  unity,  and  terms  containing 
Xq  as  a  factor  can  be  safely  neglected  in  comparison  to  both  In  Xq 
and  numbers  of  the  order  of  unity.  This  is  important  for  the  investiga¬ 
tion  of  the  theoretical  velocity-defect  profiles  Ti(y)  in  the  region  near 
the  wall.  Inasmuch  as  near  the  wall  the  mean  velocity  tends  to  go  to 
zero  in  proportion  to  ln(x/xo ) ,  we  must  consider  as  the  effective  duct 
radius  the  range  from  0<y<l-Xo  orxo^xll.  Examples  of  r)o 
values  obtained  by  putting  y  =  1  -  Xq  and  NTy  =  1  -  Xq/Z,  etc. ,  in 
various  analytical  expressions  for  r){y)  are  presented  in  the  following 
section. 

8.3  Survey  of  Previous  Theories 

In  a  critical  survey  of  existing  theories  of  the  mean-velocity  dis¬ 
tribution  of  fully  turbulent  duct  flow,  it  is  practical  to  organize  them 
according  to  the  specific  F-functions  employed,  regardless  of  whether 
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or  not  the  original  authorities  quoted  their  F-function  explicitly.  For 
convenience  of  comparison,  a  few  of  the  F -functions  discussed  here 
are  illustrated  in  Fig.  1,  and  the  resulting  velocity-defect  profiles  in 
Fig.  2. 

One  of  the  simplest  family  of  analytical  forms  for  F(y)  is  repre¬ 
sented  by  the  "power  law," 

«  -1  ,,  , 

F  =  a  *y  (6a) 

which,  upon  combination  with  equation  (3),  yields  (with  exception  of 
m  =  3/2 ) 

^  2m)/2  ,  _  2a/^3.2m)  .  (6b) 

It  is  historically  interesting  that  one  of  the  oldest  Icnown  solutions, 
namely  that  given  more  than  100  years  ago  by  Darcy  (1858),  is  of  this 
type,  for  the  special  case  of  m  =  0  in  equations  (6).  For  this  assump¬ 
tion  F  is  independent  of  the  distance  from  the  center  of  the  duct,  and 
Darcy's  value  of  the  constant  F  =  l/a  was  0.1212  .  It  is  readily  seen 
that,  if  m  >  0 ,  F  increases  monotonically  from  zero  at  the  center  to 
a  maximum  value  at  the  wall.  For  m  <  0 ,  F  is  smallest  at  the  wall 
but  approaches  infinity  when  y  approaches  zero.  A  special  case  is 
m  =  i  for  which  the  velocity  distribution  is  linear  in  y .  The  excluded 
case,  m  =  3/2  ,  yields  r\  proportional  to  In  y,  which  is  entirely 
unrealistic,  since  empirical  data  suggest,  at  least  for  the  region  near 
the  wall,  proportionality  of  x]  to  ln(l  -  y),  i.  e. ,  In  x. 

In  view  of  the  structure  of  equation  (3)  it  can  be  expected  that 
another  family  of  basically  simple  analytical  solutions  Ti(y)  or  ti(x) 
is  obtained  when  the  F-function  is  a  "modified  power  law," 

_  -lii  i  -1  rn  i  ‘  v 

F  =  a  (1  -y)  y  =  a  X  y*  (7a) 

The  combination  of  (7a)  with  equation  (3)  yields  (with  the  exception  of 
m  =  1  which  will  be  discussed  later) 

Tl  =  'no(l-x  );  or  Tio -Ti  =  u/uy  =  TioX  ,  (7b) 

where  tiq  =  a/(l  -m)  . 

The  simple  power  law  of  mean-velocity  distribution,  equation  (7b),  is 
still  in  use  as  a  convenient  interpolation  formula,  normally  with  an 
m-value  of  approximately  6/7  ,  so  that  the  exponent  1  -  m  in  equa¬ 
tion  (7b)  is  approximately  1/7. 
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The  disadvantage  of  the  power  form,  equation  (7b),  Is  that  both 
constants  In  equation  (7a),  the  exponent  (m)  and  the  factor  (a),  were 
found  to  depend  on  local  conditions.  For  example,  near  the  wall,  m 
varies  with  roughness  parameter  and  also  with  x. 

Prandtl  (1925,  1932)  has  shown  that  most  of  the  discrepancies 
are  removed  if  one  considers  the  case  of  m  =  1  in  equation  (7a)  and 
the  resulting  logarithmic  distribution  of  mean  velocity,  {r\Q  -  ii)  In  x. 
This  is  the  asymptotic  case  of  the  power-form  equation  (7b)  for  1  -  m 
approaching  zero.  Prandtl  also  showed  that  for  m  =  1  in  equation  (7a) 
the  remaining  constant  is  of  a  universal  nature,  i.  e. ,  independent  of 
the  local  geometrical  and  Icinematlc  conditions  of  the  flow  type.  He 
demonstrated  that  this  is  consistent  with  the  assumption  that  the 
"mixing -length"  near  the  wall  will  increase  in  direct  proportion  to  the 
distance  from  the  wall.  All  this  together  can  be  expressed  by  the  re¬ 
quirement  that 


lim  (-8F/dy)  =  constant  =  k  ,  (8) 

y-1 

where  k  is  named  Kdrm^n' s  constant.  For  both,  turbulent  flow  in 
pipes  and  in  the  atmosphere,  the  most  frequently  quoted  value  in  the 
literature  is  k  =  0.40;  reference  is  made  to  Schlichting  (1955);  Hinze 
(1959,  page  516)  uses  k  =  0.41  ;  and  other  authorities  values  which 
maybe  somewhere  between  0.3  and  0.5. 

Ih  view  of  the  uncertainty  of  the  constant,  numerical  relationships 
will  be  derived,  as  far  as  possible,  without  specifying  the  value  of  k  , 
For  example,  the  velocity-defect  profile  will  be  quoted  numerically  as 
KTi  rather  than  .  A  new  appraisal  of  Kdrmdn' s  constant  will  be  pre¬ 
sented  in  Section  8.6,  and  it  shall  be  shown  that  k  is  a  true  constant 
only  for  fully  developed  flow  in  rough  pipes  with  a  "best"  value  of 
K  =  0.428. 

In  the  author's  opinion,  it  is  desirable  that  the  physical  nature 
of  this  constant  be  more  convincingly  demonstrated,  and  its  value 
established  on  the  basis  of  first  principles.  Possibly  this  can  be 
achieved  by  assuming  a  variety  of  values  for  k  and  showing  that  a 
particular  value,  in  the  vicinity  of  0.4  ,  will  lead  to  minimum  condi¬ 
tions  in  the  simultaneous  processes  of  continuous  generation  and 
decay  of  eddy  energy. 

Discarding  the  "power  laws,"  equations  (6)  and  (7),  we  turn  to 
families  of  F-functions' which  satisfy  equation  (8).  Among  the  simplest 
of  such  forms  are  the  "linear  law," 
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F  =  ic(l  -  y)  =  KX  ,  (9d) 

whereupon,  in  equation  (3), 

KT|  =  -2>ry  +  ln(l  -  ln(l  -^/T)  .  {9b) 

with 

KT]o  s  -0.6137  -  In  Xo  ,  (9c) 

and  the  "parabolic  law,"  / 

F  =  k(1  -  y^)/2  =  kx(2  -  x)/2  ,  (10a) 

whereupon,  in  equation  (3), 

KT]  =  -2  tan"Wy  +  In  (1  +^/T)  -  ln(l  -  'sTy  )  ,  (10b) 

.  with 

icq©  =  -0.1845  -  in  x©  (10c) 


Equation  (9a)  was  first  suggested  by  Prandtl.  However,  in  solving 
for  the  mean-velocity  profile,  Prandtl  disregarded  the  shear-stress 
variation  with  r,  and  discussed  the  strictly  logarithmic  profile,  which 
is  a  correct  solution  of  (3)  only  if  instead  of  equation  (9a)  we  consider 
a  "modified  linear  law," 

F  =  k(1  -y)'^y  =  kxn/i-x  ,  (11a) 

whereupon,  in  equation  (3), 

KTi  =  -ln(l  -  y)  with  kti©  =  ln(l/xo)  .  (llb,c) 

Equation  (11b)  has  an  extremely  simple  analytical  structure,  and  gives 
for  the  dimensional  velocity  profile,  as  a  function  of  the  distance  from 
the  wall,  u(x)  =  (uJ'A)  Infx/x©)  ,  This  purely  logarithmic  distribution 
is  taken  frequently  as  a  standard  or  reference,  in  comparison  with  more 
complex  forms,  such  as  equations  (9b)  and  (10b),  or  the  following  the¬ 
oretical  expressions.  The  difference  between  actually  observed  dis¬ 
tributions  KTiQjjg  and  equation  (11b)  is  sometimes  referred  to  as  the 
"correction-term  function; "  see,  for  example,  Hinze  (1959,  page  518). 
In  a  revised  form,  this  concept  will  be  discussed  and  utilized  in 
Section  8.5  , 

It  follows  from  the  above  that  Prandtl  has  implicitly  introduced  a 
modification  of  the  linear  law  —  i.  e.  ,  the  F-function  equation  (11a)  — 
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on  the  basis  of  explicitly  quoting  only  the  resulting  Kn-function, 
equation  (lib).  Kto^n  (1930)  has,  also  implicitly,  introduced 
another  modification  of  the  linear  law, 


F  =  2k(1  -  N/y )'^  =  2k(1  -  ^1  -x)n/  1  -x  ,  (12a) 

whereupon,  in  equation  (3), 

kt|  =  -'sTy  -  ln(l  -  '^Ty)*  with  kt^q  =  -0.3068  -  In  Xq  .  (12b.c) 

Kdrmdn' s  original  and  explicit  starting  point  was  a  certain  "similarity 
hypothesis"  for  the  length-scale  of  turbulence  which  relates  the  F- 
function  to  the  first-  and  second-order  derivatives  of  the  velocity 
profile  in  the  following  equation 

F  =  kti'A"  .  (12d) 

Upon  differentiation  of  equation  (3) 

=  -F‘/F  +  l/2y  .  (3') 

When  F  is  given  by  equation  (12a),  differentiation  with  respect  to  y 
yields 

F'  =  -K  +  F/2y  or  F' /F  =  -  k/F  +  l/2y  ,  (12a') 

and  a  combination  of  equations  (12a' )  with  (3’ )  shows  that  equation 
(12d)  is  satisfied.  Thus,  equation  (12a)  is  necessary  and  sufficient 
for  satisfying  the  requirement  of  Kdrmdn' s  similarity  hypothesis.  How¬ 
ever,  the  correction-term  function  resulting  from  equations  (12b,c)  is 
not  in  good  agreement  with  reality,  as  will  be  demonstrated  in  Section 
8.5  .  Therefore,  the  value  of  Kdrmdn' s  similarity  hypothesis  for  tur¬ 
bulent  shear  flow  is  questionable. 

It  will  also  be  shown  in  Section  8.5  that  none  of  the  above  dis¬ 
cussed  solutions  produces  tolerably  small  and  unsystematic  correction 
terms.  Inasmuch  as  the  F-functions  in  the  region  near  the  wall  are  ad¬ 
justed  to  satisfy  equation  (8),  the  reason  for  the  deficiencies  of  the 
theoretical  n-distributions  must  be  explained  by  unrealistic  behavior 
of  the  F-function  in  the  region  near  the  center  of  the  duct.  Exactly  at 
the  center  —  i.  e.  ,  for  y  =  0  —  K^rmdn' s  and  Prandtl'  s  modified  linear 
laws  —  i.  e.  ,  equations  ( 12a)  and  ( 1  la)  —  give  Fq=0;  the  parabolic 
law,  equation  (10a),  gives  Fc  =  *</2 ;  in  Prandtl' s  original  linear  law, 
equation  (9a),  F^  =  k  .  Empirical  data  seem  to  indicate  that,  in  reality, 
Fc  should  lie  between  0  and  k/2  ,  or,  more  specifically,  near  k/3. 


MEAN- VELOCITY  PROFILE  FOR  TURBULENT  FLOW  IN  CIRCULAR  DUCT 
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In  1932,  Nikuradse  has  referred  to  "an  interpolation  formula  by 
Prandtl"  by  introducing  what  we  can  call  the  "extended  parabolic  law," 

F  =  -  y*)[l  -  c(l  -  y*)]  =  j[l-c-(l-2c)y'  -cy*]  .(13a) 

Equation  (13a)  satisfies  equation  (8)  for  any  value  of  the  constant  c. 
Nikuradse  suggested  the  value  of  c  ^  0.3  so  that 

F  =  k(1  -  c)/2  =  0.35  k  . 
c 

It  is  convenient  to  define,  as  an  abbreviation.  c/(  1  -  c)  =  a^ . 
with  a  =  0.6546  for  c  =  0.30,  or  c/(l-c)  =  3/7.  Using  a^  in 
place  of  c/(l  -  c).  we  can  reformulate  equation  (13a)  as  follows 

2F l-c  +  cy*  ■  l-y*  ^  l  +  a*y*  ’  ' 

whereupon,  in  equation  (3), 


KTl  =  KT)  + 

'  'par 


^  2tan-^(l+V2i5^)  -  Ztan 


‘(1  -V2ay)J 
(13b) 


with 


KT^o  =  0.3046  -  In  Xq 


(13c) 


where  is  defined  by  kt|  in  equation  (10b),  i.  e. ,  for  the  para¬ 

bolic  law.  The  resulting  t] -distribution,  equation  (13b)  is  rather 
complicated  but  it  is  not  surprising  that  it  produces  relatively  small 
correction  terms.  However,  unless  the  new  constant  c  in  equation 
(13a)  can  be  justified  on  the  basis  of  theoretical  reasonings,  or  shown 
to  be  of  universal  nature,  equation  (13b)  must  be  regarded  as  a  semi- 
empirical  interpolation  formula. 


8.4  A  Generalized  Mathematical  Model 


It  has  been  explained  in  Section  8.2  that  the  objective  of  the  gen¬ 
eralization  is  to  establish  an  F-function  which,  basically,  applies  to 
both  turbulent  duct  flow  and  atmospheric  boundary  layer  flow.  The 
fact  that  there  is  no  geometrically  defined  duct  diameter  in  the  atmos¬ 
phere  (because  only  the  lower  boundary  is  material)  has  an  important 
consequence.  Namely,  it  is  imperative  to  utilize  the  F-function  as 
F(x)  rather  than  F(y),  and  it  is  necessary  to  consider  F(x)  not  only 
for  the  limited  region  Xq  <  x  <  1  but  also  for  x  >  1  .  The  physical 
reason  is  as  follows:  even  for  a  perfectly  constant  horizontal  pressure 
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gradient,  the  shear-stress  distribution  in  the  atmosphere  is  not  a 
linear  function  of  the  lateral  distance  from  the  wall.  Owing  to  the 
Coriolis  forces  in  the  atmosphere,  the  horizontal  shearing  stress  is 
a  vector  and  changes  both  intensity  and  direction  with  height.  Both 
horizontal  components  of  the  shearing  stress  asymptotically  approach 
zero  at  considerable  height,  in  relation  to  the  vertical  changes  of  the 
two  horizontal  components  of  wind  shear.  It  is  natural  to  assume  that 
the  length- scale  of  turbulence  can  be  defined  in  atmospheric  boundary- 
layer  flow,  and  that  it  is  independent  of  the  orientation  of  a  horizontal 
reference-frame  of  coordinates.  By  virtue  of  its  concept,  F  would 
then  be  the  same  for  the  two  horizontal  components  and  only  positive 
F-values  can  have  physical  reality,  for  any  x  >  0 . 

In  view  of  these  requirements,  Prandtl' s  and  Kdrmdn'  s  modified 
linear  laws  —  equations  (1  la)  and  (12a)  —  must  be  discarded  because 
both  result  in  imaginary  F-values  when  x  >  1 .  The  parabolic  law  — 
equation  (10a)  —  and  Nikuradse's  extended  parabolic  law  — equation 
(13a)  — are  eliminated  because  both  result  in  negative  F-values  for 
X  >  2  .  The  only  F-function  in  Section  8.2  which  is  positive  for  any 
value  of  X  is  equation  (9a),  F  =  kx  ,  Prandtl' s  original  linear  law. 
However,,  the  resulting  irdistribution  was  not  satisfactory  for  duct 
flow.  For  the  atmosphere,  an  ever-increasing  scale  of  turbulence 
with  height  is  unlikely,  and  not  suggested  by  observational  experience. 

Rossby  (1932)  has  modified  Kdrmdn's  similarity  hypothesis  in  an 
attempt  to  apply  it  to  atmospheric  and  oceanic  turbulence.  He  found 
that  the  length-scale  of  turbulence  should  increase  with  approach 
towards  the  boundary,  which  is  exactly  the  opposite  of  Prandtl '  s 
concept.  In  patching  together  the  two  solutions,  i.  e. ,  by  assigning 
Prandtl' s  form  to  the  lowest  part,  and  Rossby' s  form  to  the  upper  and 
larger  part  of  the  atmospheric  boundary -layer,  Rossby  and  Montgomery 
(1935)  succeeded  in  explaining  several  features  of  actual  wind  profiles. 
However,  their  solution  is  unsatisfactory  owing  to  lack  of  continuity 
in  the  composite  solution,  as  well  as  a  certain  arbitrariness  in  the 
individual  solutions. 

To  a  certain  degree  one  could  draw  a  parallel  to  the  state  of 
knowledge  concerning  the  laws  of  black-body  radiation  at  the  end  of 
the  last  century.  Two  different  asymptotic  expressions  were  derived 
originally  from  classical  physics.  For  small  wave-length  Wien' s 
radiation  law  gave  an  increase  of  intensity  with  wave-length.  For 
large  wave-length,  the  Rayleigh-Jeans  law  of  radiation  gave  a  decrease 
of  intensity  with  wave-length.  It  was  Planck' s  great  achievement  to 
conceive  the  universal  law  which  contained  the  two  previously  known 
solutions  as  asymptotic  cases. 


Unfortunately,  in  contrast  to  classical  physics,  conditions  of 
turbulence  structure  are  still  only  poorly  understood.  We  are  restric¬ 
ted  to  more  or  less  heuristic  accounts  of  the  problem.  It  can  be  con¬ 
cluded  that  the  Immediate  problem  is  to  generalize  the  approach  indi¬ 
cated  by  Rossby  and  Montgomery. 

The  search  for  a  universal  similarity  principle  and  a  subsequent 
F-function  which  is  applicable  to  both  turbulent  duct  flow  and  atmos¬ 
pheric  boundary  layer  flow,  has  led  —  Lettau  (1958)  —  to  another  ex¬ 
tension  of  Prandtr  s  basic  linear  law,  which  may  be  referred  to  as  the 
"extended  linear  law,"  or  the  "linear-hyperbolic  law," 

F  =  k(1  -  y)/[  1  +  m(l  -  =  kx/(1  +  mx^^^^^)  .  (14) 

which,  when  differentiated  once  with  respect  to  x .  gives 

8F/8X  =  r(l  -  +  m  .  (14' ) 

It  can  be  readily  verified  that  equation  (8)  is  satisfied  for  any  value 
of  m  .  F  reaches  its  maximum  value,  for  any  m ,  at  x  =  1 ,  i.  e. , 

(y  =  0),  and  approaches  monotonically  zero  for  increasing  x,  but 
always  F>0  for  0<x<oo. 

Equation  (14)  contains  only  one  additional  number  (m)  in  addition 
to  the  universal  constant  k.  Comparison  with  observational  findings 
suggests  that  m  equals  Z  for  the  one-dimensional  duct-flow  while 
m  equals  4  for  the  two-dimensional  atmospheric  boundary  flow.  The 
fact  is  highly  significant  that  the  parameter  m  seems  to  be  determined 
by  the  flow  structure.  It  means  that  there  remains  no  arbitrary  or 
empirical  constant  in  equation  (14).  For  an  illustration  of  the  two  F- 
functions  see  Fig.  3.  Thus,  for  duct  flow, 

i  i. 

F  =  k(1  -  y)/[l  +  2(1  -  y)^]  rxx/d  +  Zx*)  .  (15a) 

whereupon,  in  equation  (3), 

KTj  =  In  ^  -  2N/y  +  (y  -  il'^y  -  y^  +  i  sin"Wy  ,  (I5b) 

1 

with 

icrio  =  0.2210  -  In  Xq  (15c) 

For  the  region  of  duct  flow,  0  <  y  <  1  -  Xq  ,  the  F' s  in  equations 
(15a)  and  (13a)  are  very  similar  to  each  other.  They  are  identical 
close  to  the  wall;  at  the  center,  equation  (15a)  yields  F^  =  k/3  =  0.333k, 


% 


Fig.  3.  The  non-dimensional  length-scale  (F)  of  turbulence 
as  a  function  of  the  non-dimensional  distance  (x)  from  the 
boundary.  Note  that  the  parameter  m  as  defined  in  the  text 
equals  l/n  in  this  graph.  For  duct  flow,  n  equals  \  (or 
m  =  2  )  and  0  <  x  <  1.  For  the  two-dimensional  atmospheric 
boundary  layer  flow  n  equals  4  (or  m  =  4 ) . 
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which  is  within  5%  of  Nikuradse'  s  F-value  for  y  s  0  ,  =  0.3Sk  . 

This  agreement  is  certainly  within  permissible  limits  of  the  empirical 
determination  of  Nikuradse' s  C9nstant  c. 

A  consequence  of  the  similarity  of  the  two  F-functions  in  the  re¬ 
gion  0  <  y  <  1  is  that  the  resulting  ^-distributions  agree  closely 
with  each  other.  However,  the  mathematical  structure  of  equation 
(15b)  is  obviously  less  complicated  than  that  of  equation  (13b).  For 
an  illustration,  see  Fig.  2,  and  for  numerical  comparisons  see  Sec¬ 
tion  8.5 . 

A  summary  of  the  various  F-functions  discussed  is  represented 
in  Fig.  1 .  It  can  be  assumed  that  those  F-curves  which  show  sharp 
comers  or  other  discontinuities  are  less  realistic  than  those  with 
rounded  forms.  Among  the  latter  group  the  F-function  given  by  equa¬ 
tion  (15a)  is  unique  in  that  it  contains  no  empirical  constant,  and 
satisfies  an  important  similarity  principle  relating  duct  flow  and 
atmospheric  boundary  layer  flow. 

8.5  Average  Velocity  Defect,  and  Comparison  with  Observations 

The  average  of  the  velocity  defect  over  the  aoss- section  in  a 
cylindrical  duct  is 

. 

ti  =  J  iidy*  .  (16) 

0 

Due  to  Xq  «  1  ,  the  integration  can  be  performed  between  the  limits 
of  zero  and  one.  In  order  to  facilitate  the  integration  let  us  consider 
the  following  transformations 

ndy^  =  d(riyM  “  y^dq  =  d[r\{y^  -  1)]  +  (1  -  y^ ) d n  .  (17) 

Now,  consider  that  dn  =  n'  dv  »  whereupon,  with  the  aid  of  equation 
(3), 

1  r  - 

KT)  =  K[Ti(y'-l)]  +  J  k(1  -  y*)y*F‘‘dy  .  (18) 

0 

The  first  member  on  the  right-hand  side  of  equation  ( 18)  is  zero,  be¬ 
cause  of  T)  =  0  for  y  =  0 .  The  second  term  can  be  readily  evaluated 
for  the  F-functions  given  by  equations  (9a,  10a,  11a,  12a,  13a,  and 
15a).  The  resulting  numerical  constants  obtained  for  kt^  are  sum¬ 
marized  in  the  following  listings;  for  comparison,  the  numerical 
constants  (kt]o  +  Inxo)  from  equations  (9c,  10c,  11c,  12c,  13c, 
and  15c)  are  also  included. 
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Fig.  4.  Observed  and  theoretical  variation  of  the  normalized 
velocity-defect  correction  term  (£=  (ri/fi)  -  ln(l  -  y)^^)  as  a 
function  of  the  non-dimensional  distance  (y)  from  the  center 
of  a  cylindrical  duct. 
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For  F  in  equation 


(  9a) 

KTi  =  16/15  =  1.  067  ; 

KT)o  +  In  Xo 

= 

-0.  614 

<10a) 

"  =  4/3  *  1. 333  ; 

M 

-0.  184 

(11a) 

"  =  3/2  =  1. 500  : 

II 

0.  000 

(Ua) 

"  =  77/60  =  1.  283  ! 

M 

= 

-0. 307 

(13a) 

"  =  1.647  ; 

II 

S 

0.  305 

(15a) 

"  =  1. 607  : 

II 

= 

0.  221 

Darcy' s  power  law,  i.  e. ,  m  =  0  in  equation  (6a)  yields  fj  =  8a/21 
=  3.  14.  It  is  significant  that  the  ratio  Kii/(Kfi)  =  r\/r\  is  independent 
of  the  particular  value  of  the  Kdrm^n  constant  k  ,  and  that  all  theoret¬ 
ical  T) -distributions  are  essentially  logarithmic  in  (1  -  y).  Therefore, 
it  is  practical  to  take  the  logarithmic  velocity-defect  distribution, 
i.  e. ,  equation  (11b)  and  the  resulting  Kfj  =  3/2  ,  as  a  standard  or 
reference  distribution. 

Let  us  define  a  "normalized  correction-term  function"  in  the  fol¬ 
lowing  equation 

£  -ln(l  -y)^  .  (19) 

The  average  value  across  the  duct,  F  ,  equals  zero.  Fig.  4  illustrates 
observed  and  theoretical  e  -distributions  as  a  function  of  the  non- 
dimensional  distance  from  the  center  of  the  duct.  The  observational 
distributions  were  obtained  from  data  presented  by  Nikuradse  (1932, 
1933)  and  Laufer  (1954).  They  agree  with  each  other  within  tolerable 
error  limits  although  the  measurements  were  made  for  quite  different 
flow  conditions.  Nikuradse  investigated  flow  of  water  in  smooth  and 
rough  pipes,  and  Laufer  the  flow  of  air  in  a  smooth  tube.  The  theoret¬ 
ical  € -distributions  show  characteristic  differences,  especially  in 
their  behavior  at  y-values  close  to  one,  i.  e.  ,  near  the  wall.  This  is 
explained  by  differences  in  the  values  of  the  numerical  constant 
(additional  to  In  Xq  )  in  the  equations  for  kt^q  which  were  found  in 
Section  4  and  are  summarized  in  the  above  listings. 

An  objective  measure  of  the  agreement  between  observed  and  the¬ 
oretical  velocity  distribution  is  given  by  the  sum,  or  average  of  the 
deviation  squares  (^obs  "  ^theor^^  '  Inspection  of  Fig.  4  shows 
that  the  £ -distribution  based  on  the  extended  linear  F  law  — i.  e.  ,  F 
as  given  by  equation  (15a)  —  evidently  produces  the  least  sum  of  devi¬ 
ation  squares^.,  ~A”similar  statement  is  true  for  F  given  by  the  extended 
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Fig.  5.  Average  velocity-defect  in  cylindrical  ducts  with  defined  wall -roughness 
(R/k)  as  a  function  of  the  Reynolds  number,  recomputed  from  empirical  data  tabu¬ 
lated  by  Nikuradse  (1933). 
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parabolic  law,  equation  (13a)»  which  is  not  shown  in  Fig.  4  because 
it  is  inherently  an  interpolation  formula. 

It  can  be  concluded  that  a  universal  velocity-distribution  applies 
to  turbulent  flow  in  smooth  as  well  as  rough  ducts,  and  that  an  en¬ 
gineering  approximation  to  such  flows  is  obtained  by  equation  (15b), 
i.  e. ,  the  proposed  generalization  of  the  mathematical  expression  for 
the  Icinematic  length-scale  of  turbulence,  as  a  "linear-hyperbolic" 
function  of  the  distance  from  the  wall. 

8.6  K^mdn*  s  Constant,  and  Relationship  between  Geometric  and 

Aerodynamic  Wall-Roughness 

The  existence  of  a  universal  velocity  distribution  implies  that 

is  a  universal  constant.  From  the  discussion  in  Section  8.5  it 
was  concluded  that  the  best  value  for  this  constant  is  Kfj  =  1.607 . 

If  this  is  accepted  it  follows  logically  that  k  will  be  a  universal 
constant  only  when  ^  is  independent  of  flow  conditions.  This  can 
only  be  tested  with  the  aid  of  empirical  data. 

Fig.  5  is  based  on  Nikuradse' s  (1933)  investigation  of  water 
flow  in  rough  pipes.  The  data  show  indeed  that,  within  tolerable 
error  limits,  ^  varies  neither  with  the  Reynolds  number  nor  wall 
roughness.  Nikuradse  gave  n  =  3.75  ,  and  a  recalculation  con¬ 
firmed  this  value,  rj  =  3.754.  Consequently,  for  the  new  form  of 
the  universal  velocity  distribution,  equation  (15b),  the  value  of 
K^rm^n'  s  constant  for  fully  turbulent  flow  is 

K  =  1.607/3.754  =  0.428  .  (20) 

The  most  frequently  .used  value  of  k  =  0.40  would  be  true  when  the 
universal  velocity  distribution  is  assumed  to  be  logarithmic,  i.  e. , 

KT]  =  3/2  and  rj  =  3.75  =  15/4  .  However,  in  view  of  the  facts  il¬ 
lustrated  in  Fig.  4,  there  is  no  justification  for  the  hypothesis  that 
the  universal  velocity  law  is  best  represented  by  the  simple  logarith¬ 
mic  form,  equation  (11b).  The  use  of  k  =  0.4  is  legitimate  only  in 
connection  with  the  logarithmic  distribution,  equation  (11b),  and  will 
lead  to  error  when  used  with  real  data. 

Because  Fig.  4  suggests  a  rather  universal  validity  of  the  same 
velocity-distribution  law,  it  is  interesting  to  discuss,  briefly,  turbu¬ 
lent  flow  in  ducts  with  smooth  walls.  Fig.  6  is  based  on  Nikuradse' s 
(1932)  investigations  of  water  flow  in  smooth  pipes.  From  his  tabula  - 
tions  of  Re,  logio  (l/^/T),  and  u/U  —  considering  that  (tio  -  n)  =  N/8/f  , 
and  u/U  =  (no  ■  h)/^o  ”3  total  of  125  values  of  n  were  recomputed 


AVERAGE  VELOCITY 


Fig.  6.  Average  velocity  defect  in  cylindrical  ducts  with  smooth  walls  as  a  function  of 
the  Reynolds  I  number,  recomputed  from  empirical  data  tabulated  by  Nikuradse  (1932). 
The  scale  at  the  right-hand  margin  indicates  the  value  of  the  Kdrm^  constant  for  the 
universal  velocity-distribution  derived  from  the  extended  linear  F-law. 
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and  plotted  in  Fig.  6.  Another  careful  evaluation  of  n  with  the  aid 
of  Laufer's  (1954)  graphically  presented  u/U-data,  for  air  flow  at 
Re  =  428.000  resulted  in  fi  =  4. 10.  while  =28.4.  Laufer's  rj 
agrees  satisfactorily  with  Nikuradse'  s  data  for  this  region  of  Reynolds 
number. 

The  n-values  for  smooth  tubes  are  significantly  larger  than 
T)  =  3.75  which  was  found  for  flow  in  rough  pipes.  Furthermore,  in 
contrast  to  Fig.  5.  the  data  on  Fig.  6  indicate  that  for  smooth 
pipes  displays  a  relatively  weak  but  significant  variation  with  Rey¬ 
nolds  number.  The  3.75  value  for  rough  pipe  seems  to  be  asymptot¬ 
ically  approached  at  extremely  high  Reynolds  number. 

This  behavior  of  n  must  of  necessity  mean  that  in  smooth-duct 
turbulence  the  Kdrmdn  constant  does  vary  systematically  with  the 
Reynolds  number  and,  therefore,  k  is  not  a  universal  constant.  The 
possible  range  of  k  is  between  0.35  and  0.42.  Its  relationship  to 
Re  can  be  estimated  with  the  aid  of  the  scale  on  the  right-hand  mar¬ 
gin  of  Fig.  6  if  equation  (15b)  is  taken  as  the  universal  velocity 
distribution  law. 

Returning  now  to  fully  developed  turbulent  flow  in  rough  pipes, 
let  us  re-analyze  the  relationship  between  geometric  and  aerodynamic 
measures  of  the  wall  roughness  with  the  aid  of  the  new  velocity- 
distribution  law.  Nikuradse  has  derived  the  empirical  relation 

Tio  =  8.48  +  5.751ogio(RA)  =  8.48  +  2.51n(RA)  .  '  (21) 


Since  Nikuradse  assumed  that  k  =  0.40 ,  equation  (21)  can  be  re¬ 
formulated  as 

Krio  =3.3y2  +  ln(R/k)  =  In  (29.7  RA)  (22) 

The  empirical  equation  (22)  was  compared  with  the  theoretical  rela¬ 
tionship  equation  (11c).  i.  e.  ,  assuming  that  the  universal  velocity 
distribution  is  actually  given  by  equation  (11b);  hence, 

Kilo  =  ln(l/xo)  =ln(NRA)  ,  (11c) 

whereupon  the  value  of  Nikuradse' s  number  N  was  obtained  as 
N  =  29.7  -  30. 


If  equation  (15b)  is  taken  as  the  best  representative  of  the  uni¬ 
versal  velocity  distribution  law,  the  empirical  equation  (22)  must  be 
compared  with  the  subsequent  theoretical  relation 
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KTio  =  0.2210  +  ln(l  -  Xo)  =  ln(1.246/xo)  ,  (15c) 

whereupon  it  follows  that  the  value  of  Nilcuradse'  s  number  is 

N  =  29.7/1.246  =  23.9  ^  24  .  (23a) 

This  establishes  the  non-dimensional  relationship  between  the  geo¬ 
metrical  and  the  aerodynamic  measure  of  wall  roughness, 

Zo  =  k/N  =  k/23.9  .  {23b) 

for  the  new  mathematical  model  of  the  universal  velocity-distribution 
law  in  fully  turbulent  flow  in  rough  ducts. 

With  the  aid  of  equations  (4)  and  (5), 

f  =  8/[-3.23  +  5.38  log,o {R/24k)] *  .  (24) 

C*'=  *)o  »  5.38  log, o(R/24k)  .  (25) 

The  last  two  equations  can  be  used  to  estimate  the  conventional  pipe- 
friction  factor  f  and  the  wall  friction-velocity  coefficient  C  for  the 
new  velocity  law. 

8.7  Turbulent  Viscosity  and  Effective  Reynolds  Number 

The  turbulent  viscosity  is  defined  as  the  quotient  of  actual 
shearing  stress  (or  x/p  )  and  rate  of  shear,  8u/8r  .  It  follows  from 
equations  (1)  and  (2)  that  vr^  must  be  a  function  of  the  radial  dis¬ 
tance  r  from  the  center,  and  that  a  non-dimensional  form  of  this 
function  is  obtained  as 


v^/Ru?  =  F'v/f  .  (26) 

For  any  of  the  F-functions  discussed  in  Sections  8.3  and  8.4,  is 
relatively  small  not  only  at  the  wall  but  also  at  the  center,  and  reaches 
a  maximum  value  in  between.  Let  us  consider  the  average  over  the 
cross-section  of  a  cylindrical  duct, 

1  1  3 

V  /Ruo  =  J  F'/y  dy'  =  2  /  Fy^  dy  .  (27) 

0  0 

For  the  purpose  of  the  following  discussion  it  is  sufficiently  ac¬ 
curate  to  approximate  the  universal  velocity  distribution  by  the  loga¬ 
rithmic  law,  i.  e. ,  F  =  k(1  -  y)V^  and,  consequently,  k  =  0.40. 
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Then,  the  turbulent  viscosity  in  equation  (26)  is  proportional  to 
y(l  -  y)  which  equals  (1  -  x)x  .  The  symmetry  of  this  form  indicates 
that  reaches  its  maximum  value  exactly  at  mid-radius,  where 
X  =  ^  =  y  .  Hence,  in  equation  (27), 

?  =  Ruo  2k/12  =  Ruf/15  .  (28) 

f 

When  F  is  expressed  by  Darcy' s  equation  (6a)  for  m  =  1  ,  the  numer¬ 
ical  factor  in  equation  (28)  would  be  3a/2  *  1/12.4  instead  of  l/l5  . 
It  has  been  pointed  out  by  several  investigators  —  for  example, 

Corrsin  (1957)  — that  it  can  be  interesting  to  utilize  a  suitably  defined 
turbulent  viscosity  for  the  definition  of  an  effective  Reynolds  number 
Reeff  .  In  order  to  preserve  its  characteristic  as  a  bulk  modulus,  let 
us  consider 

ROeff  =  2Ru/v^  =  Re  .  (29) 

With  the  aid  of  equation  (28),  and  remembering  that  u/u?  =  %  -  n. 
equation  (29)  yields 

Re  =  30u/u?  =  30(iio  -  r\)  .  (30) 

eii 

Fully  developed  turbulent  flow  requires  that  the  actual  Reynolds 
number  must  be  in  excess  of  at  least  10,000  ,  but  the  effective  Rey¬ 
nolds  number  as  defined  by  equation  (30)  will  be  less  than  1,000 
which  is  considered  to  be  the  order  of  a  critical  value.  For  example, 
in  Laufer's  experiment  mentioned  in  Section  8.5  and  8.6,  iio  =  28.4  , 
and  fi  =  4. 10  ,  so  that  Reeff  =  729  ,  while  Re  =  428,000  . 

For  ducts  with  significant  wall  roughness  r\  equals  3.75  which 
is  a  truly  universal  constant.  Then,  the  magnitude  of  Re^^f  is 
uniquely  determined  by  the  value  of  no  •  Equation  (30)  shows  that 
Reeff  is  below  1,000  when 

i 

^0  <  (3.33  +  ti)  »  37  =  14.8/0.4  (31a) 

From  the  results  discussed  in  Section  6.  we  substitute  ln(l/xo)  for 
KTio.  with  k  =  0.40,  and  x,  =  lc/30R  ,  whereupon 

ln(l/x„)  <  14.8  »  ln(3.10‘)  (3lb) 

or 

l/30xo  =  (RA)  <  1.0  XlO*  .  (31c) 

The  physical  meaning  of  the  last  inequality  can  be  interpreted  as  fol¬ 
lows.  Fully  turbulent  flow  can  possibly  produce  Re^^^  <  1,000  only 
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when  k/R  is  larger  than  10’  .  This  seems  indeed  to  be  a  reasonably 
close  approximation  to  the  practical  limit  for  which  wall  roughness 
has  been  found  significant:  reference  can  be  made  to  any  textbook  of 
basic  fluid  mechanics  —  for  example*  Rouse  and  Howe  (1953). 

When  the  walls  of  the  duct  must  be  considered  "smooth"  the  ef¬ 
fective  Reynolds  number  in  equation  (30)  will  remain  below  a  critical 
value  of  approximately  1*000  if  rj  is  larger  than  3.75.  This  is 
actually  the  case  as  was  shown  in  Section  8.5 .  Perhaps  it  is  worth¬ 
while  to  investigate  these  relationships  more  thoroughly.  It  appears 
possible  that  the  success  of  a  turbulent  viscosity  is  not  as  super¬ 
ficial  as  it  has  been  often  assumed. 

8.8  Conclusion 


This  paper  is  intended  to  be  primarily  a  heuristic  account  of  the 
specific  problem  of  mean-velocity  distribution  in  fully  turbulent  flow. 
A  generalization  of  a  similarity  principle  is  introduced  and  it  can  be 
concluded  that  certain  aspects  of  the  proposed  new  mathematical  ex¬ 
pression  for  the  mean  length-scale  of  turbulence  merit  closer  investi¬ 
gation.  It  is  perhaps  possible  that  certain  parallels  to  developments 
in  classical  physics,  mentioned  in  Section  8.4,  are  physically  mean¬ 
ingful.  The  structure  of  the  suggested  F-function  (see  Fig.  3)  and  the 
interpretation  of  its  characteristic  parameter  (m)  as  indicating  the 
"degrees  of  freedom"  (or  the  number  of  dimensions)  makes  it  similar 
to  the  structure  of  probability  functions  —  for  example,  the  Maxwell 
distribution  of  molecular  speeds.  This  feature  can  help  to  consoli¬ 
date  the  dual  picture  of  turbulence,  in  building  a  bridge  between  the 
engineering  concepts  (which  are  deterministic,  and  based  on  bulk 
parameters  and  mean-flow  properties)  and  the  theoretician' s  concepts 
(which  are  quantum-mechanistic  and  based  on  turbulence  structure, 
statistical  parameters,  and  spectra). 

Aside  from  these  philosophical  aspects  there  are  several  questions 
which  are -capable  of  more  direct  answers.  For  example,  a  detailed 
study  of  wind  profiles  in  the  atmospheric  boundary  layer,  in  terms  of 
the  generalized  similarity  principle  and  the  subsequent  F-function, 
has  been  completed  and  is  discussed  in  Section  9.  Also  under  way 
are  investigations  of  the  structure  of  shear  flow  in  flat-plate  boundary 
layers. 
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Theoretical  Wind  Spirals  in  the  Boundary  Layer 
of  a  Barotropic  Atmosphere 


Heinz  H.  Lettau 

Department  of  Meteorology 
University  of  Wisconsin 


Abstract.  A  new  theoretical  form  for  the  length-scale  of 
turbulence  as  a  function  of  distance  from  the  lower  bound¬ 
ary  is  proposed  which  contains  no  empirical  constant  other 
than  K^rmdn'  s  constant.  An  integration  of  the  equation  of 
motion  with  height-dependent  friction-velocity  is  performed. 

The  resulting  universal  wind-spiral  solution  for  the  atmos¬ 
pheric  boundary  layer  is  tabulated.  It  can  be  verified  for 
given  sets  of  Vg  =  geostrophic  wind,  f  =  Coriolis  param¬ 
eter,  and  Zo  =  aerodynamic  roughness  parameter  of  the 
ground  surface.  The  solution  is  completed  with  the  analy¬ 
tical  formulation  of  various  coefficients  (geostrophic  drag 
coefficient,  unit-height  coefficient,  etc.  )  as  unique-valued 
functions  of  the  surface-Rossby  number  Roq  =  Vg/zof.  In 
the  atmospheric  surface  layer  the  solution  is  equal  to  the 
conventional  logarithmic  wind  profile.  Computed  wind 
4.  spirals  compare  satisfactorily  with  observed  wind  spirals 
for  the  boundary  layer  of  an  adiabatic  atmosphere.  Agree¬ 
ment  is.  also  found  between  computed  and  empirical  surface- 
stress  relationships.  Applications  of  the  theory  in  esti¬ 
mating  the  dissipation  of  energy  in  the  lower  troposphere 
are  discussed. 

9. 1  Introduction 

The  discussion  is  restricted  to  cases  of  steady  and  horizontally 
uniform  large-scale  air  motion  above  level  terrain  of  uniform  and  de¬ 
fined  surface  roughness.  At  any  point  in  the  lower  troposphere  a 
perfect  balance  between  the  pressure  gradient  force,  the  Coriolis 
force,  and  the  friction  force  is  assumed.  This  excludes  inertia  forces, 
and  vertical  components  of  the  mean  air  motion..  The  resulting  mean 
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Wind  is  a  strictly  horizontal  vector  ( V )  independent  of  horizontal  co* 
ordinates,  but  possibly  changing  direction  and  intensity  from  level 
to  level.  Furthermore,  barotropy  is  assumed,  i.  e.  the  air  density  p 
is  a  unique-valued  function  of  atmospheric  pressure.  This  excludes 
"thermal  winds"  so  that  the  geostrophic  wind  (IL)  is  independent  of 
height. 

An  atmospheric  boundary  layer  will  exist  in  which  the  wind  in¬ 
creases  with  height  from  the  boundary  condition  of  Vo  -  0  at  the 
surface  (z  =  0)  asymptotically  approaching  the  geostrophic  wind  at 
sufficiently  large  distance  from  the  ground  (z  >  H) .  In  comparison 
with  the  height  of  the  atmosphere  (or  the  tropopause),  the  thickness 
(H)  of  the  atmospheric  boundary  layer  is  sufficiently  small  to  justify 
the  neglecting  of  vertical  density  differences:  this  implies  that  terms 
like  V  8p/dz  are  small  in  comparison  with  p  d^dz  and  can  be  disre¬ 
garded  in  the  boundary  layer. 

When  f  s  Coriolis  parameter  and  r  =  vector  of  shearing  stress 
at  height  z ,  the  equation  of  motion  reduces  for  the  above  formulated 
assumptions  to 

.  pfl(V- V  )  =  8t'/8z  .  (1) 

where  complex  vector  notation  Is  used,  1  =  .  It  Is  assumed  that 

the  vector  of  shearing  stress  is  horizontal  and  at  any  level  proportional 
to  the  vector  of  shear  (dV/dz) .  Then  the  factor  of  proportionality  is 
a  scalar  quantity  (function  of  height)  which  defines  the  effective 
dynamic  viscosity  (pj^ .  K  is  referred  to  as  the  eddy  viscosity  or 
diffusivity.  When  8Vg/8z  =  0  ,  the  equation  which  defines  K  can  be 
written  in  terms  of  geostrophic  departure, 

7  =  pKav/az  =  pK8(V- V^)/8z  .  (2) 

It  is  useful  to  transform  the  variables  and  the  equations  into  di¬ 
mensionless  forms.  Let  X  denote  a  parameter  having  the  physical 
dimension  of  a  length;  X  will  be  referred  to  as  the  scale-height  of 
the  boundary  layer.  The  dimensionless  form  of  the  independent  vari¬ 
able  will  be  denoted  by  x  =  z/X  ,  and  differentiations  with  respect 
to  X  by  primes.  Using  the  symbol  T  for  the  horizontal  unit  vector 
in  the  direction  of  the  geostrophic  wind,  we  define  dimensionless 
wind  components  (U,  u,  v),  dimensionless  shearing -stress  components 
(s,  t),  and  a  dimensionless  eddy  viscosity  (y)  by  the  following  four 
identities: 
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T(u  +  lv)  *  V/Xf 
Tu  ■  VgAf 
T(s+lt)  *  'r/pX*!* 

and  y  ■  K/X*f  . 

With  the  aid  of  the  dimensionless  independent  and  dependent  variables, 
the  two  vector  equations  (1)  and  (2)  yield  a  scalar  system  of  four  simul* 
taneous  first-order  differential  equations, 

s  =  y u*  ;  s'  =  -V  ;  t  =  y  v'  ;  t'  =  u  -  U  ,  (3a) 

which  is  equivalent  to  a  scalar  system  of  two  simultaneous  second- 
order  differential  equations, 

s  =yt"  ;  t  =  -y  s"  ,  (3b) 

which  also  is  equivalent  to  one  scalar  fourth-order  differential  equa¬ 
tion, 


s  =  -y(ys")"  .  (3c) 

It  can  be  readily  verified  that  equation  (3c)  applies  for  t  in  the  same 
form  as  for  s .  The  transformation  to  corresponding  differential  equa¬ 
tions  in  u  and  v  is  possible  but  results  in  more  complicated  expres¬ 
sions  than  that  for  s  and  t  . 

The  boundary  conditions  at  large  height  (x  »  1 )  are 

s  =  t  =  s'=t'  =  0  , 


while  at  the  surface, 

s  =  So>0;  t  =  to>0;  (s')o=0;  (t')o  = -U  . 

In  view  of  the  later  discussion  of  conditions  at  rough  boundaries,  with 
a  defined  aerodynamic  roughness  parameter  Zq  of  the  earth' s  surface, 
it  is  useful  to  define  the  independent  variable  in  dimensionless  form 
as  X  =  (z  +  Zq  )/X  .  This  does  not  affect  the  above  definitions  and 
equations.  The  only  consequence  is  that  the  lower  boundary  conditions 
(at  z  =  0)  now  refer  to  x  =  x©  =  Zq  /X,  which  is  a  very  small  number  in 
comparison  to  unity. 
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The  problem  is  to  integrate  the  basic  equations  (1)  and  (2),  or 
one  of  the  equivalent  systems  (3),  in  order  to  obtain  the  vertical  pro¬ 
files  (spirals)  of  geostrophic  departure  and  shearing  stress  in  the 
atmospheric  boundary  layer.  Obviously,  the  solution  depends  ex¬ 
clusively  on  the  analytical  formulation  of  the  height  dependency  of  K , 
or  y  =  y(x) .  The  simplest  possible  case,  that  of  K  =  Kq  »  const  (or 
y'  =0),  results  in  the  classical  Ekman-splral  solution.  For  K  being 
a  linear  or  an  exponential  function,  or  a  power  law  of  height,  closed 
solutions  exist  and  are  given  by  Bessel  functions  (Hankel  functions). 
For  K  being  a  quadratic  function  of  height,  the  closed  solution  is  ex¬ 
pressed  by  elementary  functions  and  is  known  as  the  Rossby-spiral 
solution.  None  of  these,  or  related  solutions,  using  explicit  height 
functions  for  the  eddy  viscosity  K,  has  produced  satisfactory  agree¬ 
ment  with  observational  data. 

The  reason  for  this  lack  of  success  can  be  explained  by  the  fact 
that  eddy  viscosity  K  is  the  product  of  a  length  times  a  velocity. 

Only  the  length  term  can  be  expected  to  be  an  explicit  function  of  the 
distance  from  the  lower  boundary,  while  the  velocity  term  will  depend 
on  the  shearing  velocity,  which,  in  turn,  depends  on  the  solution  of 
the  basic  differential  equations.  Consequently,  the  problem  involves 
a  characteristic  degree  of  non-linearity. 

It  is  the  purpose  of  the  following  discussion  to  derive  a  new  solu¬ 
tion  of  the  atmospheric  boundary  layer  problem.  The  new  theory  is 
based  on  concepts  which  have  been  found  useful  in  the  treatment  of 
fully  developed  flow  in  straight  ducts  —  Lettau  (1961);  see  Section  8. 

9.2  Theoretical  Model  Hypotheses 

The  absolute  value  of  the  horizontal  vector  of  shearing  stress  is 
a  scalar  quantity  t  which  decreases  monotonically  with  height,  from 
the  surface-stress  value,  tq,  to  zero  at  the  top  of  the  boundary  layer. 
The  scalar  shearing  velocity  is  defined  as  ^  r/p  .  It  is  assumed 
that  the  eddy  viscoscity  K  in  a  scalar  form  of  equation  (2)  can  be 
expressed  as 

T/(paV/9z)  =  K  =  i  n/  t/p  ;  thus,  f  =  ^ i/p/av/dz  .  (4a) 

This  establishes  the  defining  equation  of  the  characteristic  length 
scale  of  turbulence,  i  .  Let  the  dimensionless  length  scale  of  tur¬ 
bulence  (f/X)  be  denoted  by  F  =  F(x).  Recalling  that  the  dimension¬ 
less  components  of  the  shearing  stress  vector  were  s  and  t ,  equa¬ 
tion  {4a)  is,  in  dimensionless  form, 
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y  =  F(s*  +  tT  =  FT  =  F(x)  •  T(x)  ,  (4b) 

where  the  convenient  abbreviation  +  t^  is  introduced.  The 

dimensionless  shearing  stress  intensity,  therefore,  can  be  written 
as  T*=  r/pX^f^  . 

The  above  definition  of  the  characteristic  length  scale  of  turbu¬ 
lence  conforms  with  eddy-spectrum  concepts  in  that  i  corresponds 
closely  to  the  largest  eddy  size  of  the  inertial  subrange.  The  scalar 
product  T*  d^/dz  =  TdV/dz  s  pc  (ergs/cm^  per  sec)  represents  the 
rate  at  which  mean-motion  energy  is  locally  converted  to  turbulent 
energy.  If  the  mean  flow  and  turbulence  energy  are  in  steady  states, 
the  local  conversion  rate  pe  must  be  equal  to  the  local  value  of  vis¬ 
cous  energy  dissipation.  It  is  generally  assumed  that  in  fully  devel¬ 
oped  turbulence  this  conversion  process  feeds  mean-motion  energy 
primarily  and  directly  only  into  eddies  of  relatively  large  size  (of 
characteristic  scale  i  )  while  viscous  dissipation  is  important  only 
in' eddies  of  much  smaller  size.  The  physical  quantities  which  de¬ 
scribe  the  latter  process  are  the  above  e  (cm^/sec^)  and  the  molecu¬ 
lar  viscosity  V  (cm^/sec).  Dimensional  reasoning  shows  that  the 
only  possible  combination  of  these  two  quantities  which  has  the  di¬ 
mension  of  a  length  is  {v^/t  With  a  numerical  factor  sufficiently 
close  to  unity  this  relatively  small  length  is  considered  to  represent 
the  lower  limit  of  the  inertial  subrange  of  eddy  sizes.  With  the  effec¬ 
tive  or  eddy  viscosity  (K),  in  place  of  molecular  viscosity  (v),  a  cor¬ 
responding  dimensional  argumentation  results  in  the  length  (KVe)*^^, 
but  this  is  the  same  as  r/p/ (8V/3z)  =  i  and  thus  can  be  con¬ 

sidered  to  represent  the  upper  limit  of  the  inertial  subrange  of  eddy 
sizes. 

A  general  hypothesis  is  introduced  now  as  follows.  F,  the  di¬ 
mensionless  length  scale  of  turbulence  in  equation  (4b)  is  assumed 
to  be  a  universal  function  of  the  dimensionless  distance  from  the 
lower  boundary,  of  the  analytical  form 

F=  Kx/(1  +  n"*x^^”)  ;  so  that  F' /F  =  nx"*  (1  -  x^^^)/(n  +  x^^^) ,  (5a) 

where  the  two  numerical  parameters  are  k  =  0.4  =  universal  Kdrmdn 
constant,  while  n  =  l/2  for  duct  flow  and  n  =  l/4  for  atmospheric 
boundary  layer  flow.  Equation  (5a)  is  an  extension  of  Prandtl'  s  linear 
law  of  the  length -scale  of  turbulence.  For  any  value  of  n  the  F- 
function  has  a  maximum  value  at  x  =  1  ,  while 

lim  F'  =  K  . 
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In  principle,  such  an  extension  of  Prandtl'  s  concept  was  first  sug¬ 
gested  and  discussed  by  Rossby  (1932)i  see  Section  8.4 . 

i 

When  n  s  1/2  ,  i.  e. ,  for  duct  flow, 

3 

F  =  0.4x/(l  +  2xM  ,  (5b) 

and  the  F-maximum  occurs  at  the  center  of  the  duct,  so  that  X  =  R  - 
radius  of  the  duct  which  is  externally  presaibed.  In  view  of  the  ex¬ 
actly  linear  dependency  of  t  on  x,  no  further  hypothesis  is  necessary 
for  the  turbulent  duct-flow  problem,  and  the  mean  velocity  profile  is 
directly  obtained  by  a  closed  integration;  reference  is  made  to  Lettau 
(1961),  see  also  Section  8.4. 

In  the  atmosphere  the  width  of  the  duct  would  correspond  to  the 
height  of  the  geostrophic  wind  level.  However,  this  level  is  not  ex¬ 
ternally  prescribed  but  depends  on  the  flow  conditions  and  must  be 
obtained  as  the  result  of  an  integration  of  the  equation  of  motion. 
Moreover,  in  contrast  to  duct  flow,  the  height-dependency  of  shearing 
stress  is  not  a  priori  known  and  must  also  be  obtained  as  a  result  of 
the  integration.  Consequently,  the  general  hypothesis,  equation  (5a), 
with  n  =  1/4  for  atmospheric  flow, 

5 

F  =  0.4x/(1  +  4xM  ,  (5c) 

must  be  supplemented  by  two  additional  hypotheses. 

Upon  logarithmic  differentiation  of  equation  (4b), 

y'/y  =  F'/F  +  T’/T  .  (6) 

With  F'/F  =  0  at  x  =  1 ,  it  follows  from  equations  (5a-c)  that 
F'  /F  >0  at  X  <  1.  However,  T'  A  will  be  negative  throughout  the 
entire  boundary  layer.  Consequently,  a  level  z*  <  X  must  exist  so 
that  at  X*  =  z*A  we  find  that  y'  /y  =  0  due  to  F'  /F  =  -T'  A  .  In 
other  words,  a  level  z*  must  exist  where  eddy  viscosity  reaches  its 
maximum  value. 

The  first  hypothesis  is  that  x*  is  universally  equal  to  1/2 .  This 
specific  choice  of  the  level  of  maximum  eddy  viscosity  at  half  the 
scale-height  (X)  is  suggested  by  empirical  findings  concerning  both 
aerodynamic  duct  flow  and  atmospheric  boundary  flow. 

The  second  hypothesis  is  that  the  level  of  maximum  eddy  viscosity 
in  atmospheric  boundary  layers  coincides  with  the  level  of  maximum 
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cross-isobar  flow.  In  non-dimensional. terms,  the  cross-isobar  com¬ 
ponent  is  V,  and  a  maximum  value  requires  that  v'  =  0  .  or,  accord¬ 
ing  to  the  third  member  of  equation  system  (3a),  t  =  0 . 

Let  starred  symbols  refer  to  values  at  the  level  of  maximum  eddy 
viscosity,  x*  =  0.5.  From  equation  (5c), 

F*  =  0.2/(l  +  A/V)  =  0.0746  , 

*  1.  X 

(T'/T)  -  2{2*  -  l)/{i*  +  1)  =  0.4322  . 

From  equation  (6),  considering  that  (y‘/y)*  =  0  , 

(T'  /T)*  =  -  (F’  /t)*  =  -0.4322  . 

With  the  aid  of  the  defining  equation,  =  s^  +  t^  ,  and  the.  subse¬ 
quent  form  obtained  by  logarithmic  differentiation,  4T'/T  = 

(2s'  s  +  2t’  t)/(s^  +  t^ )  ,  for  t*  =  0  ,  one  finds  that 

T*  =  n/F  ;  and  2(T' A)*  =  (s’ /s)*  , 


whereupon. 


y*  =  F*  T*  =  0.0746  n/F”  ; 
and 

(s’/s)*  =  -2(F’ A)*  =  -0.8644  . 

In  summary,  three  function-values  are  fixed  relative  to  s* ,  at 
X  =  X*  ,  namely, 

s*  =  s*  ;  t*  =  0 
(s' )*  =  -0.8644  s* 

(s’’)*  =  -t*/y*  =  0  ;  (t")*  =  s*/y*  =  13.40n/F  . 

A  fourth  value  is  at  our  disposal,  which  amounts  to  a  choice  of  the 
ratio  (s' /t’)*  ,  or  (t’)*  . 

9.3  Method  of  Integration 

Given  the  function  values  of  s  and  t,  and  their  derivatives  s' , 
t' ,  s",  and  t",  at  x  =  x* ,  we  want  to  solve  a  two-parameter  family 
of  curves,  s(x)  and  t(x),  by  an  integration  of  equations  {3b).  With 
the  data -summary  at  the  end  of  the  preceding  section,  the  problem  is 
to  find  for  what  values  (if  any)  of  s*  and  (t’)*  we  obtain  a  spiral 
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solution  in  s(x)  and  t(x)  which  converges  monotonically  to  the  point 
(0,  0) ,  according  to  the  boundary  conditions  at  large  x. 

For  the  actual  solution,  the  simple  Cauchy  method  of  numerical  ' 

integration  was  employed,  writing 

s' (x  +  Ax)  =  s’ (x)  +  Ax  [s"(x)  +  s’'(x  +  Ax)]/2  , 

s(x  +  Ax)  =  s(x)  +  Ax  [  s'(x)  +  s' (x  +  Ax)]/2  , 

and  two  corresponding  forms  for  t'  (x  +  Ax)  and  t(x  +  Ax) . 

Given  a  set  of  values  (x,  F,  s,  s',  s",  t,  t',  t",  y),  we  calcu¬ 
late  a  new  set  for  x  +  Ax  by  the  following  procedure: 

(a)  F(x  +  Ax)  is  calculated  directly. 

(b)  Using  the  value  of  s"(x)  as  a  first  approximation  for  s"(x  +  Ax) 
we  calculate  s'  and  s ,  at  x  +  Ax,  by  the  Cauchy  formula. 

(c)  Using  the  value  of  y(x)  we  obtain  t"  =  s/y  at  x,  and  using 
t"(x)  as  a  first  approximation  for  t"(x  +  Ax)  we  calculate  t'  and  t, 
at  X  +  Ax,  by  the  Cauchy  formula. 

« 

(d)  With  the  aid  of  the  above  obtained  values  of  s  and  t  at  x  +  Ax, 
we  calculate  new  values  of  y  and  s"  at  x  +  Ax  by  y  =  F(s^  +  t^)*^*, 

and  s"  =  -t/y  .  ^ 

(e)  In  the  fashion  of  a  true  iteration,  we  repeat  the  procedures  (b) 
through  (d),  using  the  improved  y  and  s" ,  until  the  values  which 

are  generated  at  x  +  Ax  agree  with  the  values  used  to  generate  them. 

Tentatively,  using  hand-computations,  it  was  found  that  con- 
vergency  of  the  spiral  solution  at  large  heights  (x  >  3)  depended 
very  delicately  on  the  proper  choice  of_the  values  of  s*  (or  y*  = 

0.0746'^'s^),  and  (t')*.  The  quality  of  the  convergence  can  be  judged 
by  plotting  the  spiral  curves  in  an  orthogonal  u  versus  v  (or  t'  ver¬ 
sus  s')  diagram.  Another  criterion  is  that  the  value  of  y  must  de¬ 
crease  monotonically  with  x,  at  any  x  >0.5. 

A  program  was  written  for  the  IBM-704  using  FORTRAN  II.  A  con¬ 
stant  increment  of  Ax  =  0.1  was  employed.  The  iteration  was  required 
to  Iceep  the  difference  between  generated  and  assumed  values  for  each  * 

step  within  0.0001%  .  The  numerical  integration  was  performed  from 
X  =  X*  =  0.5  up  to  X  =  4.5  ,  requiring  40  steps.  During  July  and  August, 
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1960,  a  total  of  72  runs  was  made,  with  a  systematic  variation  In 
values  for  y*  and  (t*)* .  The  following  example  Illustrates  the  narrow 
tolerance  with  which  the  quality  of  convergence  was  Judged.  When 
y*  was  kept  equal  to  0.8318,  (t')*  =  -93.284  resulted  In  slight  over¬ 
shooting,  (t’)*  =  -93.291  In  slight  undershooting  of  the  desired  con¬ 
vergence  point  (0, 0),  at  X  larger  than  3.5 .  Values  of  (t’)’^  larger 
or  smaller  than  the  above  mentioned  figures,  and  also  any  slight 
variation  of  the  y***  value,  produced  spiral  distortions,  and  spurious 
minimum  values  of  y  at  large  x.  It  was  concluded  that  the  converg¬ 
ence  criterion  has  absolute  and  objective  significance. 

After  finding  the  correct  values  at  x  =  x*  =  0.5  ,  that  is, 

s*  =  124.41  ,  (s')*  =  -107.54  , 

t*  =  0.00,  and  (f  )*  =  -  93.29  , 

the  solution  was  completed  by  a  numerical  Integration  of  equations 
(3b)  downward  from  x  =  0.5  to  x  =  0.0050  .  Eleven  steps  with  a 
constant  increment  of  Ax  =  0.045  were  used.  At  x<  0.0050  log¬ 
arithmic  extrapolation  was  employed. 

9.4  Integration  Results  —  Universal  Wind  Spiral 

To  facilitate  a  practical  application  of  the  wind  profile  computa¬ 
tion,  the  results  are  presented  not  only  as  sets  of  x,  F,  s,  s' ,  s", 
t,  t' ,  t",  y  values  (see  Table  1)  which  were  directly  written  out  in 
great  detail  by  the  IBM-704,  but  also  after  certain  transformations 
were  made.  We  recall  that  s'  and  t'  represent  the  two  rectangular 
components  of  the  dimensionless  geostrophic  departure.  Let 

T  =  ^  +  t'^  and  p  =  tan”*  (s' /t')  ,  (7) 

so  that  r  is  the  absolute  value  and  p  the  azimuth  of  the  dimension¬ 
less  geostrophic  departure.  These  functions  are  listed  in  Table  2. 

The  dependent  variable  is  quoted  in  two  forms,  x  =  (z  +  Zq)/X  and 
^  =  (z  +  Zq)/Z  with  4  =  2x  or  Z  =  X/Z  .  The  scale-height  of  the 
boundary  layer  (X)  was  defined  as  the  level  where  the  length-scale 
of  turbulence  has  the  maximum  value;  the  unit-height  Z  is  defined  as 
the  level  where  eddy  diffusivity  has  the  maximum  value  and  where,  in 
the  basic  solution,  the  cross -isobar  flow  has  a  maximum  value,  too. 

For  any  assumed  set  of  lower  boundary  values  (ro,  Po).  the  two 
rectangular  wind  components  parallel  with,  and  perpendicular  to  the 
geostrophic  wind,  are 


Table  1.  Basic  spiral  solution  as  a  function  of  relative  altitude,  x  s  z/X;  F  =  relative  scale  of  turbulence; 
s,  t  =  rectangular  components  of  relative  stress;  s',  s",  t',  t"  =  first  and  second-cx-der  derivatives  of  s  and  t 
with  respect  to  x  ;  y  =  FT  =  dimensionless  eddy  diffusivity  (where  T^  =  +  t* )." 
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u/U  =  1  -  {t/tq  )  •  cos  (P  -  Po  ) »  and  v/U  =  (t/tq  )  •  sin  (p  -  Po )  • 

When  these  two  terms  are  multiplied  by  the  speed  of  the  geostrophic 
motion  (Vg)»  the  dimensional  wind  components  are  obtained. 

In  addition  to  =  absolute  value  of  the  dimensionless  shearing 
stress  we  define  its  azimuth 


Y  =  tan'*(t/s)  .  (8)  * 

Values  of  T^,  T,  y  and  =  p  -  y  are  also  listed  in  Table  2.  The 
data  of  T^  and  y  serve  to  compute  the  stress  spiral,  considering  that 
the  relative  stress  components  parallel  with  and  perpendicular  to  the 
surface  stress  are 

(T*AJ)cos(y- Vo),  and  (TVT?)  sln(Y  -  Yo)  • 

The  angle  \\t  formed  by  the  vectors  of  geostrophic  departure  and  shear¬ 
ing  stress  shows  a  significant  height  dependency  at  x  <  0.5 .  It  will 
be  remembered  that  in  the  classical  Ekman-spiral  this  angle  is  con¬ 
stant  and  equal  to  45®. 

In  order  to  compute  an  actual  wind  spiral  it  is  necessary  to  spec¬ 
ify  the  external  conditions.  As  a  minimum  requirement,  this  means 
that  a  set  of  three  parameters  must  be  given  —  the  geostrophic  speed 
Vg ,  the  Coriolis  parameter  f ,  and  the  aerodynamic  surface  rough¬ 
ness  Zq  .  The  non-dimensional  combination  of  these  three  quantities 
is  the  surf  ace -Rossby  number  (see  Lettau  (1959)), 

Roo  =  Vg/zo  f  .  (9) 

Owing  to  the  boundary  condition  of  zero  wind  speed  at  the  surface, 
the  quantity  ro  (which  is  the  absolute  value  of  the  dimensionless 
geostrophic  departure  at  z  =  0 )  must  be  equal  to  the  dimensionless 
geostrophic  speed  (U).  Recalling  that  Xq  =  Zq/X  ,  it  follows  that 

ROo  =  U/xo  =  Tq/xq  .  (10) 

Consequently,  with  the  aid  of  the  listings  in  the  lower  part  of  Table  2, 
we  can  establish  the  functional  relationship  between  r/x ,  and 
r,  T^  ,  T,  or  the  angular  quantities.  Obviously,  at  z  =  0,  equals 
^0  ,  since  the  wind  vector  in  the  immediate  vicinity  of  the  surface 
will  be  parallel  with  the  surface  stress,  while  oq  is  the  angle  formed 
by  the  wind  nearest  to  the  surface,  and  the  geostrophic  motion.  It 
can  be  noted  that  r/x  varies  by  many  powers  of  ten.  Therefore,  the 


Table  2.  Universal  spiral  solution  for  geostrophic  departure  and  stress 
as  a  function  of  relative  height.  (  s  z/Z  =  2x :  r,p  s  relative  magnitude 
and  azimuth  of  geostrophic  departure ;  ,  y  -  relative  magnitude  and 

azimuth  of  horizontal  stress;  T  s  relative  magnitude  of  friction  velocity; 

r  p  -  Y  s  angle  formed  by  the  vectors  of  stress  and  geostrophic  de¬ 
parture. 


X 

( 

r 

(dea) 

T* 

Y 

(deal 

T 

(deo) 

4.  0 

8.0 

1.66 

377.8 

0.27 

323.6 

0.  52 

54.2 

3.8 

7.6 

2.74 

338.5 

0.  59 

286.  8 

0.  77 

51.7 

3.6 

7.2 

4.  45 

304.4 

1.  13 

255.0 

1.  06 

49.  4 

3.4 

6.8 

6.  74 

276.7 

1.96 

227.4 

1.  30 

49.  3 

3.2 

6.4 

9.68 

252.5 

3.  19 

203.3 

1.78 

49.2 

3.0 

6.0 

13.3 

231.  1 

4.  92 

181.  8 

2.  22 

49.3 

2.8 

5.6 

17.  1 

211.9 

7.25 

162.  5 

2.  69 

49.4 

2.  6 

5.2 

23.0 

194.5 

12.67 

145.  1 

3.  56 

49.4 

2.4 

4.8 

29.0 

178.6 

14.  27 

129.  1 

3.  78 

49.  5 

2.2 

4.4 

36.  1 

164.0 

19.21 

114.5 

4.  38 

49.  5 

2.  0 

4.0 

43.0 

150.5 

25.3 

101.  0 

5.  03 

49.  5 

1.  8 

3.6 

53.8 

137.9 

32.9 

88.4 

5.73 

49.  5 

1.  6 

3.2 

63.6 

126.2 

41.  4 

76.  6 

6.44 

49.6 

1.4 

2.8 

74.8 

115.2 

52.  1 

65.6 

7.  22 

49.  6 

1.2 

2.4 

86.  9 

104.  9 

64.4 

55.  3 

8.  02 

49.  6 

1.  (f 

2.0 

100.6 

95.  0 

78.7 

45.  5 

8.  87 

49.  5 

0.9 

1.8 

108.3 

90.3 

86.  6 

40.8 

9.  31 

49.  5 

0.8 

1.  6 

116.  1 

85.7 

96.7 

36.  2 

9.83 

49.  5 

0.7 

1.4 

124.2 

81.  1 

104.  3 

31.7 

10.  43 

49.  4 

0.6 

1.  2 

133.  0 

76.  6 

114.  0 

27.  3 

10.  68 

49.  3 

0.  5 

1.  0 

142.4 

72.  0 

124.  2 

22.  9 

11.  14 

49.  1 

0.  4 

0.8 

152.7 

67.4 

135.6 

18.  5 

11.  64 

48.  9 

0.  3 

0.6 

164.  5 

62.6 

147.2 

14.  3 

12.  15 

48.  3 

0.2 

0.4 

179.0 

57.  2 

160.  0 

10.  1 

12.  65 

47.  1 

0.  1 

0.2 

199.7 

50.  1 

173.  9 

5.6 

13.  16 

44.  5 

5.  10"^ 

1.  10*' 

219.4 

45.  6 

180.6 

3.2 

13.  45 

42.  4 

5.  lO'* 

1.  10** 

283.  1 

34.  1 

187.9 

0.4 

13.  70 

33.  7 

5.  10"* 

1.  10*' 

351.  1 

27.  0 

188.  5 

0.  1 

13.  74 

26.  9 

5.  10'* 

1.  10*'* 

413.7 

22.  1 

188.6 

0.  0 

13.75 

22.  1 

5.  10'‘ 

1.  10‘* 

497.4 

18.7 

188.6 

0.  0 

13.  75 

18.  7 

5.  10"’ 

1.  10*‘ 

572.5 

16.2 

188.6 

0.  0 

13.75 

16.  2 
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common  logarithm  of  the  surface -Rossby  number  was  utilized  for  the 
interpolations  which  resulted  in  the  data  summarized  in  Table  3.  With 
the  aid  of  Table  3,  the  boundary  values  (ro>  To*  Tqi  Po*  andyo) 
can  be  interpolated  for  a  given  set  (Vg»  Zo>  f)  or  the  subsequent 
surface -Rossby  numbers. 

At  this  point,  an  order-of-magnitude  discussion  will  be  illustra¬ 
tive.  While  both  Vg  and  f  will  normally  vary  by  no  more  than  a  fac¬ 
tor  of  4  (or  1/4)  about  their  average  values  (Vg  approximately 
1,000  cm/sec  and  f  approximately  0.00005  sec'M*  the  aerodynamic 
roughness  parameter  Zq  of  natural  ground  can  vary  from  more  than 
100  cm  (for  tall  vegetation)  to  less  than  0.  01  cm  (for  snow  surfaces, 
etc. ).  Consequently,  the  surface -Rossby  number  may  range  from 
approximately  lO’  to  10^.  This  is  considered  in  the  computations 
which  are  summarized  in  Table  3  and  the  following  tables,  where 
I  certain  variables  are  quoted  as  a  function  of  the  surface-Rossby 
number. 

Table  3.  Dependence  of  parameters  of  theoretical  wind  and  stress 
spirals  on  surface-Rossby  number  ^0  •  To  =  relative  magnitude  of 
the  surface  stress;  ao  =  angle  formed  by  surface  stress  and  geo- 
strophic  wind;  ro ,  Po  =  relative  magnitude  and  azimuth  of  the  surface 
geos  trophic  wind;  Yo  =  Po  *  <3^0  ^  To  =  relative  magnitude  of  the  sur¬ 
face  friction  velocity. 


4.  5 

187.  6 

35.  0 

269.  2 

36.  0 

1.  0 

13.  67 

5.  0 

188.  1 

31.  9 

299.3 

32.2 

0.  3 

13.  72 

5.  5 

188.  3 

28.  9 

330.  4 

29.  0 

0.  1 

13.  73 

6.  0 

188.  5 

26.  1 

362.  4 

26.  1 

0.  0 

13.  73 

6.  5 

188.  6 

23.  7 

395.  2 

23.7 

0.  0 

13.  73 

7.  0 

188.  6 

21.  7 

428.7 

21.  7 

0.  0 

13.  73 

7.  5 

188.  6 

20.  1 

462.  8 

20.  1 

0.  0 

13.  73 

8.  0 

188.  6 

18.  7 

497.4 

18.7 

0.  0 

13.  73 

8.  5 

188.  6 

17.  4 

532.  4 

17.  4 

0.  0 

13.  73 

9.0 

188.  6 

16.  2 

567.  8 

16.  2 

0.  0 

13.  73 

9.  5 

188.  6 

15.  1 

603.  5 

15.  1 

0.  0 

13.  73 
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For  the  conversion  of  dimensionless  distance  from  the  surface 
(x ,  or  I  =  2x )  into  absolute  or  metric  heights  we  consider  again  that 
r  :  U  s  Vg/Xf  s  Vg/2Zf .  Thus,  a  unit-height  coefficient,  iro  i 
can  be  defined  in  the  relation 

Z=  C-  V /f  i  or  X  =  2C„  V/f  .  (11) 

2  g  »  9 

It  follows  from  the  universal  wind-spiral  solution  in  Table  2  or  Table  3 
that  the  unit-height  coefficient  G2  is  a  unique -valued  function  of  the 
surface -Rossby  number.  Numerical  values  are  listed  in  Table  4.  With 
them  we  can  compute  the  unit -height  Z,  or  the  scale -height  of  the 
boundary  layer  X  =  2Z  ,  for  any  given  set  (Vg,  Zq,  f) . 


Table  4.  Non-dimensional  coefficients  of  theoretical  wind  spirals  as 
functions  of  the  surface-Rossby  number  Roo  .  C  =  geostrophic  drag 
coefficient;  0%  ~  unit  height  coefficient;  =  maximum-diffusivity 
coefficient;  C|y4  =  coefficient  of  mass  transport  across  the  isobars. 


log  10  Roq 

C 

10*  •  C2 

10’- Ck 

Cm 

4.5 

0.  0508 

187.4 

116.  9 

1.  13 

5.0 

.  0459 

167.  1 

92.9 

1.  06 

5.  5 

.  0416 

151.  3 

76.  2 

1.  00 

6.0 

.  0379 

138.  0 

63.  4 

.  94 

6.  5 

.  0347 

126.  5 

53.  2 

.  88 

7.0 

.  0320 

116.  7 

45.  2 

.  82 

7.  5 

.0297 

108.  1 

38.  8 

.77 

8.  0 

.0276 

100.5 

33.6 

.  72 

8.  5 

.  0258 

93.9 

29.3 

.  68 

9.0 

.  0242 

88.  1 

25.  8 

.  64 

9.  5 

.  0228 

82.  8 

22.  8 

.  60 

As  an  example,  let  f  =  10"*  sec"^  ,  Zq  =  5  cm ,  and  Vg  =  15.8 
m/sec.  Then  ROq  =  3.16  10^  or  logio  B9o  =6.50,  and  Vg/f  = 
1.57  10^  m.  For  this  surface-Rossby  number.  Table  4  gives  C2  = 
126.5  10"^  ,  whereupon  Z=  126.5  X  1.57  =  199  m.  The  unit -height 
Z  will  have  the  order  of  magnitude  of  10^  m  in  nearly  all  cases; 
exceptions  will  only  occur  for  extremely  low  geostrophic  speeds. 
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In  a  similar  manner  the  absolute  magnitude  of  the  surface  drag 
(tq  )  or  the  surface-friction  velocity  ('s/  tq  /p  )  can  be  computed  from 
To  /U  =  To  Ao  =  Vto  /p  /Vg  .  This  relationship  Is  the  basis  for  the 
definition  of  the  geostrophlc  drag  coefficient  C  =  To  Ao  in  the 
equation 


To  /p  =  CV  ;  or  To  =  p  C*  V*  .  (12) 

0  g 

It  follows  from  the  universal  wind-spiral  solution  that  the  geostrophlc 
drag  coefficient  C  is  a  unique  function  of  the  surface-Rossby  number. 
Numerical  values  are  listed  in  Table  4.  They  permit  us  to  compute  to 
(or  n/  To  /p  )  for  any  given  set  (Vg,  z©  ,  f). 


The  listings  in  Table  4  include  also  the  maximum  eddy- viscosity 
(or  diffusivity)  coefficient  Ck  which  is  defined  in  the  relation 


K  =  vVf 

max  K  g 


(13) 


where  ~  maximum  value  of  K  which  was  assumed  to  occur  at 

X  =  X*  =  0. 5  or  ^  =  1. 0 ,  1.  e. ,  z  =  Z .  Another  coefficient  is  that  of 
total  mass  transport  across  the  isobars,  which  is  defined  in  the 
relation 

00 

/  Vg(v/U)d2  =  C^ZVg  .  (14) 

where  v/U  =  (tAo  )  •  sin  (p  -  Po)  =  dimensionless  cross-isobar  flow, 
and  Z  =  unit-height,  as  defined  above.  This  “mass  transport  across 
the  isobars,  due  to  friction  and  the  subsequent  imbalance  between 
pressure  gradient  force  and  Coriolis  force,  is  well  known  for  any  solu¬ 
tion  of  atmospheric  boundary  layer  flow  —  for  example,  the  Ekmah- 
spiral  solution. 


Inspection  of  the  values  listed  in  Table  4  reveals  that  all  coeffi¬ 
cients  decrease  significantly  with  increasing  surface-Rossby  number. 
It  can  be  concluded  that  surface  drag,  boundary-layer  scale  height, 
maximum  eddy  viscosity,  and  mass  transport  across  the  isobars,  are 
significantly  dependent  on  the  aerodynamic  roughness  parameter  Zq  of 
the  terrain  under  consideration. 


Several  characteristic  levels  (relative  to  unit-height,  Z  )  of  the 
atmospheric  boundary  layer  are  listed  in  Table  5,  as  a  function  of  the 
surface-Rossby  number.  The  height  of  the  surface  layer,  h  ,  denotes 
that  level  at  which  the  wind  speed,  computed  from  the  spiral  solution, 
begins  to  deviate  by  more  than  5%  from  the  conventional  logarithmic 
law.  More  detail  will  be  discussed  in  connection  with  data  listed  in 
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Table  6.  The  geostrophic  displacement  thickness  6  is  defined  in 
analogy  to  the  displacement  thickness  of  flat-plate  boundary  layer 
flow.  The  wind  component  parallel  to  the  geostrophic  wind,  1.  e. . 
u/u  =  1  -  (r/ro ) '  cos  (p  -  Po )  is  considered,  and  6  is  defined  in  the 
relation 

'  00  00  X 

/  V  udz  i  /  V  Udz  ;  or  6  =  /  (1  -  u/tj)dz  .  (15) 

0  ®  -  6  ®  0 

6  can  be  interpreted  as  the  level  above  which  an  exactly  geostrophic 
motion  would  produce  the  same  total  mass  transport  (in  the  direction 
parallel  to  the  isobars)  as  is  produced  by  the  actual  wind  above  the 
level  z  =  0 .  In  the  unidimensional  shear  flow  above  a  flat  plate  in 
the  wind  tunnel,  the  air  between  0  and  6  can  be  thought  of  as  being 
brought  to  rest;  however,  atmospheric  boundary  layer  flow  is  two- 
dimensional  (spiral  flow),  and  it  is  necessary  to  consider  the  air  be¬ 
tween  the  surface  and  the  level  6  as  moving  parallel  to  the  direction 
of  the  pressure  gradient  from  high  to  low  pressure;  in  fact,  this  layer 
can  be  thought  of  as  absorbing  the  total  cross -isobar  flow,  which  was 
expressed  with  the  aid  of  the  coefficient  Cm  in  Table  4.  Obviously, 
horizontal  variations  in  this  total  cross -isobar  flow  will  result  in 
frictionally  induced  mean  vertical  motion  at  the  top  of  the  atmospheric 
boundary  layer. 


Table  5.  Relative  height  of  characteristic  levels  as  functions  of  the 
surf  ace -Rossby  number  Roo  (relative  to  unit  height  Z  =  Cz  )• 
h  =  height  of  the  surface  layer  (for  5%  deviation  from  log-law) ; 

6  =  geostrophic  displacement  thickness  ;  Zy  =  level  of  maximum  cross - 
isobar  flow ;  H  =  height  of  the  "geostrophic  wind  level." 


log  10  Roo 

h/Z 

6/Z 

Zy/Z 

h/z 

4.  5 

0.  106 

0.  265 

1.  59 

5.  69 

5.  0 

.  273 

.  242 

1.  42 

5.  61 

5.  5 

1  .  340 

.221 

1.  28 

5.  53 

6.  0 

.  403 

.  202 

1.  16 

5.  46 

6.  5 

.  465 

.  185 

1.  05 

5.  40 

7.  0 

.  526 

.  171 

.  96 

5.  35 

7.  5 

.  587 

.  159 

.  88 

5.  31 

8.0 

.  648 

.  148 

.  81 

5.  28 

8.  5 

.  709 

.  138 

.75 

5.  25 

9.0 

.  770 

.  129 

.  70 

5.  22 

9.  5 

.  830 

.  122 

.  65 

5.  20 

NORMALIZED  WIND  SPIRALS  FOR 
VARIOUS  SURFACE- ROSSBY  NUMBERS  (Ro  ) 


Fig.  2.  Normalized  theoretical  wind  spirals  for  extremely  low, 
moderate,  and  extremely  high  surface-Rossby  numbers. 
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The  level  of  maximum  arose -isobar  flow  is  also  listed  in  Table  5. 
It  will  be  remembered  that,  in  the  basic  spiral  solution,  this  level 
was  made  to  coincide  with  z  s  Z.  However,  Introducing  actual 
^  surface-reference  values  causes  this  level  to  become  a  function  of 

the  surface-Rossby  number. 

^  The  universal  spiral  solution  listed  in  Table  2  shows  that  the  ' 

angle  p  increases  continuously  with  height.  For  every  Po  there 
must  exist  special  levels  where  p  -  Po  equals  180*,  then  3.60*,  540*, 
etc.  At  any  of  these  special  levels  the  actual  wind  vector  is  parallel 
to  the  geostrophlc  wind  vector.  The  first  of  these  special  levels,  i.e., 
the  height  at  which  p  -  Po  =  180* ,  is  conventionally  referred  to  as 
the  "geostrophlc  wind  level,"  H,  even  though  only  the  direction,  but 
not  the  speed,  is  truly  geostrophlc.  Values  of  H  can  be  obtained 
from  the  ratio  H/Z  listed  in  Table  5.  It  is  Interesting  to  note  that  H , 
like  Zy  and  6 ,  increases  with  decreasing  surface-Rossby  number, 
while  the  height  of  the  surface  layer,  h,  behaves  in  the  opposite 
manner. 

A  direct  consequence  of  the  particular  choice  for  the  F-functlon 
in  the  boundary  layer  solution  is  that  the  wind  speed  profile  in  the 
lowest  air  layer  approaches  the  conventional  logarithmic  law, 

*  (V  /U)  u^  +  v^  a  k"*  n/  To/p  log  (z  +  zo  )/zo . 

g  e 

The  deviation  from  the  logarithmic  law  is  negligible  very  close  to  the 
'  surface,  but  increases  gradually  with  height,  due  to  three  causes: 

(1)  actual  veering  of  wind  direction  with  height,  (2)  existing  vertical 
gradient  of  shearing  stress,  and  (3)  non-linear  increase  of  the  i  term 
with  height.  Table  6  shows  at  which  levels  the  deviation  reaches 
values  of  2,  5,  and  10%.  All  these  h-levels  are  strongly  dependent 
on  the  surface-Rossby  number.  Previously,  estimates  of  the  height  of 
the  surface  layer  were  based  on  the  relative  momentum  divergence, 
i.  e. ,  on  a  certain  value  of  (t©  -  t)/to  ;  reference  can  be  made  to 
Lettau  (1957a).  In  this  case,  a  possible  non-linearity  in  the  height 
increase  of  the  length  scale  of  turbulence  and  also  possible  changes 
of  wind  direction  with  height  were  disregarded.  It  follows  from  the 
^  data  in  Table  6  that  the  effect  of  shear-stress  gradients  seems  to  be 
partially  compensated  by  the  effect  of  the  height  decrease  of  the 
gradient  of  length  scale  of  turbulence.  Evidence  of  this  is  the  fact 
that  the  listed  values  of  the  ratio  At/t©  are  nearly  always  signifi¬ 
cantly  larger  than  the  percentage  figures  denoted  by  "i"  in  Table  6. 

»  The  only  exception  occurs  at  lowest  surface-Rossby  number,  i.  e. , 

most  likely  for  extremely  rough  ground  surface.  Certain  consequences 
of  these  findings,  with  respect  to  wind  profile  curvature  in  the  adia- 
‘  batic  surface  layer,  will  be  discussed  in  a  forthcoming  report. 


tyilMCI  STRCtt  vimut  iurfacc  uostiiofhic  ifiid 


V-n  (M/SEC) 
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Fig.  3.  Surface  stress  as  a  function  of  surface  geostrophic 
speed.  (1)  From  wind  spiral  theory  and  geostrophic  drag 
coefficient  relationships;  (2)  after  Mintz  (1957);  and 
(3)  after  Elidssen  (1960). 
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Table  6.  Surface-layer  characteristics  as  functions  of  the  surface- 
Rossby  number  Roo .  h^  =  level  at  which  the  wind  speed  computed 
from  the  spiral  solution  deviates  by  i%  (for  i  =  2,  5,  and  10)  from  the 
logarithmic  law  (for  which  V  =  2. 5  To  Vg  ^  logQ  (x/xo ));  Z  s  unit  height 
“  Cz  Vg^(/f ;  At/to  =  AT*/(T*)o  =  (to  -*Th)/To  =  relative  stress  gradi¬ 
ent,  or  total  momentum-flux  divergence  between  the  levels  0  and  h^. 

log  10  Roo 

i  = 

2% 

i 

*  5% 

i 

s  10% 

h/Z 

At/to 

hj/Z 

At/to 

hi/Z 

AtAo 

4.5 

0.01 

0.  004 

0.  11 

0.  044 

0.37 

0.  141 

5.0 

.09 

.036 

.27 

.  087 

.59 

.214 

5.  5 

.  12 

.050 

.34 

.  122 

.73 

.258 

6.0 

.  15 

.062 

.40 

.  151 

^  3^' 

.  300 

6.  5 

.  18 

.073 

.  46 

.  175 

.99 

.341 

7.0 

.21 

.084 

.  53 

.  196 

1.  12 

.  381 

7.5 

.24 

.094 

•  59 

.  216 

1.  25 

.418 

8.0 

.26 

.  104 

.65 

.236 

1.  38 

.452 

8.  5 

.29 

.  114 

.  71 

.  256 

1.  51 

.481 

9.0 

.32 

.  124 

.  77 

.276 

1.64 

.  505 

9.5 

.35 

.  134 

.  83 

.  296 

1.77 

.  524 

9.5  Comparison  with  Observations  —  Applications 

With  the  aid  of  the  universal  spiral  solution,  i.  e. ,  the  r  and  p 
values  listed  in  Table  2,  together  with  the  boundary  values  ro ,  po  In 
in  Table  3,  a  wind  spiral  for  any  given  set  (Vg,  Zq,  f)  can  be  con¬ 
structed.  Figure  1  illustrates  the  general  procedure,  and  the  relation¬ 
ships  among  the  variables  of  the  spiral  solutions,  for  two  different 
^  boundary  values.  Figure  2  shows  normalized  wind  spirals  for  three 

different  surface-Rossby  numbers.  The  normalization  involves  the 
turning  of  the  geostrophic  wind  vector  into  a  horizontal  direction  and 
constructing  the  spirals  of  geostrophic  departure  in  the  relative  terms 
of  r/ro  and  p  -  po  •  These  theoretical  spirals  are  in  satisfactory 
»  agreement  with  observed  spirals  over  land  (for  low  Roo  )  and  over  the 

ocean  (high  Roo  );  reference  is  made  to  the  illustration  of  the  "Leipzig 
wind  profile"  and  the  "Scilly  wind  profile"  in  Lettau  (i^957b.  Figure  3). 
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Fig.  4.  Seasonal  contrast  of  dissipation  of  energy  In 
the  atmospheric  boundary  layer.  Computed  from  sea- 
level  geos  trophic  winds  over  the  North  American  con¬ 
tinent  in  January  and  July  1954,  with  the  aid  of  geo- 
s trophic  drag  coefficients. 
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Table  7.  A  quantitative  comparison  between  observed  and  calculated 
wind  spiral  data  for  both  land  and  ocean  boundaries. 

Observations;  Leipzig  (08  to  15  GOT  on  20  October  1931): 
f  =  1.14  10*^  sec**  ,  yg,o=  1'7*51  m/sec,  estimated  Zo  =  30  cm; 
consequently,  Roo  =  5.13  10^ —  Ocean  near  Scilly  Islands  (12  to 
15  GOT  on  4  January  1951):  £  =  1.11  10*^  sec**  ,  Vg,o  =  12.19  m/sec, 
estimated  Zo  =0.03  cm;  consequently,  Rgo=  3.66  10^ . 

Theory:  Data  calculated  from  coefficients  and  ratios  interpolated 
for  logic  Roo  values  of  5.71  (for  Leipzig)  and  8.56  (for  Scilly)  on 
Tables  3,  4,  and  5:  assumed  air  density  p  =  0.00115  g/cm*. 


Quantity 


Leipzig  (Land)  Scillv  (Ocean) 

Observed  Calculated  Observed  Calculated 


oio  (deg)  =  angle  be¬ 
tween  surface  wind 


and  isobars 

26.  1 

27.7 

13.  9 

17.  3 

To  (dynes/cm^)  = 
surface  stress 

5.  31 

5.  65  ■ 

0.  92 

1.  13 

Z(m)  =  unit  height 
(level  where  K  =  K^ax) 

235 

223 

100 

102 

Kjjjax(cmVsec)  =  maxi¬ 
mum  eddy  diffusivity 

144,700 

205,000 

28.000 

41,  000 

2y(m)  =  level  of  maxi¬ 
mum  cross -isobar 
wind  component 

250 

274 

70 

75  - 

H(m)  =  geostrophic  wind 
level  (first  level 
where  cross -isobar 
wind  component 
equals  zero) 

1070 

1200 

480 

534 

^  maxim)  =  (KVe)^ 
at  z  =  2Z 

32 

36 

13 

16 
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i 


Fig.  5.  Ocean-continent  contrast  of  dissipation  of  energy 
in  the  atmospheric  boundary  layer.  Computed  from  sea-level 
geostrophic  winds  over  the  American  continent  and  the  At¬ 
lantic  Ocean  in  January,  1953,  with  the  aid  of  geostrophic 
drag  coefficients. 
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D»88IPAT»0M  OF  ENEROY,  PER  UNIT  MEA  (E) 

NORTH  AMERICAN  CONTINENT  (ASSUMING  1^-100  CM) 
JANUARY  AVERAGES  FOR  INDICATED  YEARS 


ESTIMATE  BY  WHITE  f  SALTZMAN  (1956)  FROM  RATE  OF 
CONVERSION  BETWEEN  POTENTIAL  AND  KINETIC  ENERGY 
(5  WATTS/M^.  AVERAGE  35*  TO  70*.  2  TO  31  JANUARY,  1953) 


Fig.  6.  Year-to-ye^r  differences  in  the  dissipation  of 
energy  in  the  atmospheric  boundary  layer.  Computed 
from  sea-level  geostrophic  winds  over  the  North  Amer¬ 
ican  continent  for  the  months  of  January  in  1952,  1953, 
and  1954,  with  the  aid  of  geostrophic  drag  coefficients 
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wind  Spirals _ _ _ 

V  •  (8^82)  =  V*  pfi(V  -  Vg)  =  -piv-  IVg  =  pfVg  •  iV  =  pfVg  •  i(V  -Vg)  = 
Vg-(9T/8z)  .  (17) 

Consequently, 

E  =  -  J^Vg- (dT/dz)dz  =  -Vg-/^(8T/8z)dz  =  Vg*To  ,  or 
0  0 

E  =  Vg  To  cos  Qfo  =  pC^^q  cos  oro  •  (18) 

Results  of  sample  computations  of  the  total  energy  dissipation, 
based  on  daily  Vg*-estimates  from  surface  weather  charts  between  30'" 
and  70”  northern  latitude,  and  a  tentatively  assumed  Zo -value  of 
100  cm  for  North  America,  and  Zo  ^  o.l  cm  for  the  North  Atlantic 
Ocean,  are  illustrated  in  Figures  4,  5,  and  6.  They  show  the  seasonal 
contrast  (January-Iuly),  the  ocean-continent  contrast,  and  differences 
from  year  to  year,  respectively.  A  basis  for  comparison  is  given  for 
North  America  in-fanuary  1953  in  that  White  and  Saltzman  (1956)  de¬ 
rived,  from  a  laborious  investigation  of  the  conversion  between  poten¬ 
tial  and  kinetic  energy  in  the  troposphere,  a  mean  dissipation  value 
of  5  watts/m^  ,  which  compares  satisfactorily  with  an  average  of  the 
meridional  distribution  of  E  depicted  on  Figure  6. 
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